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Exact theory of self-expanding piston rings 
By C. Truesdell 


1. Formulation of the problem. A piston ring must be shaped, when free of stress, in such a 
way that when inserted within a given cylinder it shall fit snugly, This requirement is idealized by 
requiring that the ring shall close exactly when subject to uniform normal pressure by the cylinder 
wall. While a new cylinder is circular, a worn one is not, so that rings for replacement ought to be 
of a somewhat different shape than were those originally installed; it is possible also that particular 
conditions of loading may make a non-circular shape desirable. In the present investigation, the 
shape of the ring is supposed to have an axis of symmetry passing through the point where the 

_ gap has been closed, but otherwise it is taken as an arbitrary plane curve (Fig. 1). Special attention, 
of course, is given to nearly circular forms. — 


———— 


Fig. 1. 


In all other respects, we adopt the same idealizations of the problem as did Meldahl? in his 


_ solution for the case of the perfectly circularring. Effects of finite thickness of the ring are neglected, 


so that the normal pressure is supposed to be applied directly to the neutral line of a Bernoulli- 


_ Euler elastica. Moreover, the flexural rigidity E J is assumed to be constant: in effect, the cross- 
_ section of the elastica is taken as uniform. 


Thus, equivalently, our problem is as follows: 


What form must be given to a cantilever of uniform flexural rigidity in order 


- that it may be constrained by a uniform normal pressure to lie upon a given curve? 


This is a problem of determining the free shape of a beam. When the deflections are small, it 


3 may be replaced by the conventional one of determining the shape assumed bya given beam when 
_ subject to given load. For piston rings which are to be used in extreme conditions, large wall pres- 
sures are needed, and to obtain them, in turn, the ring must be fabricated with a large ratio of gap 
to mean diameter. For such rings, there will be large deflections near the gap, just where require- 
- ments of performance are most critical, and the linearized theory is inadequate for precise design. 


This was recognized by Meldahl in his analysis for circular rings. 

_ The above problem of free shape is solved here exactly, in the theory of finite deflection 
of aninitially curved elastica. 

At the end of Part II the exact solution, specialized to the case of elliptical cross-section, is made 


_ the basis of estimations of the error resulting from use of certain approximate theories. From this 


comparison it is inferred that while the theory of small deflections leads to results significantly in 


1 Meldahl, Teknisk Tidskrift 73 (1943) p. 99. 
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error, a method of perturbation about the exact solution for the circular ring yields formulae which 
are quite accurate for finite deflection of elliptical rings of small eccentricity. 
In Part IIL is developed a general theory of nearly circular rings. 


2. Historical origin of the problem of free shapes. In the theory of finite elastic deformations, 
the problem of determining the shape a body must have in order to assume a given shape when load- 
ed is different from that of determining the form taken on by a given body in consequence of load- 
ing. This was recognized by James Bernoulli in the very first paper on the mathematical theory 
of elasticity, published in 16941; in the list of ,,remaining“ problems at the end is included ,,what 
should be the curvature of the band in order that from an attached load or from its own weight or 
from both together it assume the form of a straight line (this would be useful in designing the arms 
of balances and axles, where it is required that the centers of the motion and of the suspended bodies 
be collinear). 

One problem of this kind had been solved, in fact, by James Bernoulli prior to his writing the 
paper, but his cryptic note stating the result was not published until 17442. In that same year 
appeared Euler’s great treatise on elastic curves ®, in which, among other problems, he took up those 
proposed by James Bernoulli and solved them. In particular, he obtained independently the 
solution of James Bernoulli for the form a band must have in order to be bent out straight by an 
assigned, normal, terminal load. The solution was expressed by each of its discoverers in terms of 
the integrals now named after Fresnel, and Euler remarked that the free shape is spiraliform. 

Apparently the non-linear problem of free shapes lay untouched thereafter until Meldahl’s 
solution for perfectly circular piston rings in 1943. In 1953 I outlined the general method 4 which is 
developed in more detail in the present study; in my earlier note, the explicit solution for a beam 
bent out straight by uniform normal pressure was given as an illustrative example. 


Part I. Theory of piston rings of arbitrary form 


3. Calculation of the bending moment resulting from uniform normal pressure. Consider a plane 
curve C to which a constant normal unit pressure P is applied (Fig. 2). We shall calculate the re- 
sultant moment, about an arbitrary point (X’, Y’), of the load acting on the arc 4B. For the mo- 


Fig. 2. 


ment dM exerted by the load on the element of arc dS at the point Q(X, Y), where the slope is: 
© = Arctan (dY/dX), we have 


dM = [—(P cos @) (X' — X) — (P sin @) (Y’ — Y)] dS 
= P[(X —X) dX +(Y—Y) dY] 
=+ Pd [(X—X")*+(Y—Y]. (1) 


1 Jak. Bernoulli, Curvatura laminae elasticae Acta erud. Leipzig J 1694. = 
» Cur icae .. es ‘ ig June 1694, p. 262—276 = Opera 1, 

Ne PiScLeipee 1 ie Fae work has been published in a German translation by Linsenbarth, Ostwalds Klassiker 
= 2 oe eas Ne: — of Varia pe Opera 2, p. 1084—1086 (1744); the editor, Nicholas I Bernoulli, 

0 was the author's nephew, was unable to see any c tion bet i 
eee ee ny connection between the problem stated and the solution 
' § L, Euler, ,,Additamentum I de Curvis Elasticis“ in Methodus inveniendi lineas c imi minini 
Pe eee Soa a Geneva, 1744 = Opera Omnia I 24, p. 231—297. Meo ean iiautanen - 

msenbarth, Ostwalds Klassiker No. 175, Leipzic 1910. Eneli i 7 

eats Pores eipzig nglish translation by Oldfather, Ellis & Brown, 

* C. Truesdell, Proc. First Midwest. Conf. Solid Mech., Univ. Illinois, p. 52—55 (1954). 
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The moment exerted by the load acting on AB is obtained by integrating this expression from 
A to B. But since P = const., the right-hand side is a perfect differential. Hence the integral of 
(1) from A to B is independent of the path connecting A und B. Hence follows the 

Equivalence principle. The resultant moments exerted by a constant normal 
pressure field on any two plane curves having the same endpoints are equal}. 

The only application we shall need is the simplest, when one of the two curves is a straight line. 
We may express the result as follows: 

Method of equivalent chords. To calculate the resultant moment exerted bya 
constant normal pressure field on any arc, about any point, replace the arc by its 
chord. 

The only application we shall need is the simplest, when the point about which the moment is 
_ taken is one end, B. We get 


1 1 


M= PL f=. — PT? (2) 


where L is the length of the chord AB. 


4, The curvature of the free ring. According to the Bernoulli-Euler theory, a bending moment M 
produces a proportional change of curvature: 


M 
Ej 7 K—x (3) 
where ~ is the curvature of the free ring, K is the curvature of the loaded ring, and E J is the 
flexural rigidity. 


In the present application, by (2) we thus obtain 


1 
oP 


(samcih Er oS 


L?. (4) 
Since the desired shape is prescribed in advance, all quantities on the right-hand side are known 
functions of arc length S. Since explicit knowledge of the curvature of a plane curve determines the 
curve itself, up to two arbitrary constants, our problem is now solved in principle. 


5. Calculation of the slope of the free ring. Let 6 be the slope of the free ring, @ the slope of 


the loaded ring. Then 
d@ dé 


__ 40 _ 4 : 5 
er ey we igi? (5). 


K 
where dS and ds are the respective elements of arc. In the Bernoulli-Euler theory, the deformation 


is inextensional: ds = dS. Hence (4) may be written 
1 
Ween 2h 6) 
dS dS EJ : 


Integration yields 
J 2 7 
0 = 6) EJ L?dsS, (7) 


where we have adjusted the arbitrary constant by taking the foot of the ring, where 0 = 0, as 
p= 0. 


Write 


il 

ae 
git) 8 
Beye (8) 


CS 


b L ici idification is given 1 i j ited. Professor Grammel 
1 f based on Stevin’s priciple of solidification is given in my earlier paper, just cited. : 
Metneine that the argument there presented was used by A. Féppl. Eq. (2) was obtained by James Ber 
noulli in No. 164 of his ,,Meditationes“*, written in 1691 but still unpublished. 
6* 


80 C. Truesdell: Exact theory of self-expanding piston rings Ingenieur-Archiy 


and write I = I(S) for the (geometrical) polar moment of inertia of the section of the loaded ring 
from S = 0 to the point of arc length S, taken about the free end (Fig. 3). Then (8) may be written 


6=0—cl. (9) 


Since J is an increasing function of S, we see from (9) mathematical proof of the evident fact 
that the maximum angular deflexion 0 — 6 occurs at the free end. 


Y 


Free end: xX 
toot 


Fig. 3. 


6. The exact shape of the free ring. Let a plane curve be given in complex form: z= x +7y 
= x(s) +iy(s). Then 
dz 


dx « Ye Se ety: 10 
Ba tig = 0089 +isind =e ‘ (10) 


Hence we may put (9) into the form 


de = OR egy > (6 a Os eo et —icr@Z ll 
Fee Fs Oe: ere == qs’ (11) 


By integration follows 


z=Z,+fe-‘!dZ, . (12) 


where we have adjusted the constant of integration so that the foot of the ring, whose affix is 
Zy = Xq +7 Yo, shall remain at the same point in both the free and loaded conditions. Separation 
of real and imaginary parts in (12) yields 


x Y Y x 
x = X, + f cos (c I) dX + fsin(c I)dY, y = Yo+ f cos (cI) dY —/fsin(e I) dX. (13) 
OG Y Yo X 


[Caution: these are line integrals. When the ring is convex and when co-ordinates are chosen 
‘asin Figure 3, dX is negative, while d Y changes sign. It is usually easier to proceed directly from 


(12). 


7. The maximum error of the linearized theory. If deviations of shape are negligible, we can 
get a sufficiently accurate result by taking the desired shape as the free shape, applying the load, 
and then reversing the displacement so obtained. Equivalently, we can reverse the load and inter- 
change the roles of the curvatures. By applying the method of the previous sections, we should 
then get in place of (4) 


Kreutz FP, (14) 


where 1 is chord-length in the free beam. Comparison with the exact result (4) shows that the 
linearized theory is valid if and only if chord lengths are not significantly altered 
by loading. 

The result of the small deflection theory can be obtained by taking the first term in the expan- 
sion of (12) or (13) in powers of c I, namely 


Z 
sw ky + f(L—icl)dZ, (15) 


or 


Y" be 
ameX+efIldY¥=x2%, ya Y—cfIldX=y. (16) 
¥} xe 
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To estimate the error made by using the small deflection theory, we subtract (15) from (12) 
obtaining 


Zz Zz 
ee pee) 02) {leone I) — 1] —i [sin (c I) —c I]} dZ. (17) 


Hence 


le — 20| x J I[cos (c I) — 1] —i [sin (c 1) —e I] |az| 


S 


Ss 
= { [cos (c I)| + |e IT—sin(c Ii] dS < [[den+4¢ 1s|dS. (18) 


0 


The error is worst at the free end. Denote by & the length of the half ring, and by & the value of 


I at the gap (i. e., § is the polar moment of inertia of the desired form of the half ring about the 
gap). Then (18) yields 


le—2| 1 


< F(e 9) 1+ 368), (19) 


When the load is small enough that the dimensionless quantity c 0 is small, then a sharp estimate 
of the worst percentage error made by using the small deflection theory is afforded by (19). If ¢ $ 
becomes appreciable with respect to 1, it is plain from (17) that the small deflection theory is no 
longer adequate. 


As an example, consider a circular ring of radius R. Since § = 22 R?, if we define a dimen- 
sionless load factor 


ap 
2 P R 
= en por eee 
| a=2eR2=2 EJ oe Ef? (20) 
we getc O=2 a, and (19) yields 
jz — | oes 3 


For a = 0.01, the worst error made in using the small deflection theory is less than 1/;%. For 
a = 0.05, however, our bound (21) for the worst error becomes 41/,%. Now our formula (21) does 
not prove that the error in the free gap as calculated by the small deflection theory will actually 
be as large as 41/,%, and the load ratio a = 0.05 is rather high. However, this result suggests 
that in some conditions now encountered the small deflection theory will not be trustworthy. 
Meldahl has given more extensive and precise comparisons which lead to the same conclusion. 


8. The quadratic theory and its maximum error. If we take one more term in the series expansion 
of the exponential in (12), we get 


Zz ; Zz 
xe of) =%+[(l-ier—ze P) dZ = 0 —+¢ { Paz. (22) 


2 


The method of estimating the error used in § 7 can be carried out again, and the result is 


lz con 20} ut CYY3 1. Ce 2 
a Seed (+48). (23) 
For the case of a circular ring we get then 
lz — 2D| 1 ee “ele 24 
R = Zuo (1 + 7 ua). (24) 


If a =0.05, the worst error is now less than 1/,%. For a =.10, the bound (24) becomes nearly 2%. 


Our numerical estimates concern only the circular ring, but since they are rather rough and 
unflattering to the theory, we can use them with confidence for any nearly circular shape. 


Ingenieur-Archiv 
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Part II. The elliptical piston ring 
9, The exact solution for the elliptical ring. Consider a ring which is to be elliptical when loaded: 
X= Reosy, Y=Ryl—@siny, (25) 


where e is the eccentricity and y is the running parameter. We shall suppose that the upper point 
of the gap occurs at y =a (Figure 4); in other words, the ring is to close perfectly when loaded, 


Y 


Fig. 4. 


and the gap occurs at the end of the major axis. (Other elliptical cases can be obtained from our 
result by simple changes of the variables.) 


For the chord-length L we have 
I? =(R-+ Reosy)? + (R yl —e* sin y) 


= ae (2 cos” 7?) + (1 — e?) (2 sin cos z¥)| =A i cos? > yp (1 — sin? . (26) 


[The case when e = 0 is equivalent to Meldahl’s Eq. (3).] Hence 
yy pak Se eae 
cI=c{ dS =cR{ 12/1 —ecosty dy =a V(y), (27) 
d 


where ais the dimensionless load ratio, again defined by (20), and where the dimensionless function 
Vv is defined by 


v(y) = ( (L-+cos &) (1 ee sin? } 8) yl —€ costé dé. (28) 
Furthermore, for the ellipse ee we have 
dZ = dX +idY =R [—siny +i/1—e cos p| dy =i R [1 — & cosy + isin y] dp 
= +i R[/1—2 (ev + ei) + e'v—e—t¥] dy 
=> iR[ + 1—e) v— (1 —1— 2) e—¥] dy. (29) 


Putting (27) and (28) into the general expression (12) yields 


‘ 5 y 
z= 1 Ps fic +1 — e) els—2V@) —(1 —y1—e) ein + eV} dy. (30) 
0 


By separating real and imaginary parts, we obtain the exact solution for the elliptical ring: 


r= ee +y1=e) [sin [7 —a V(y)] dy — (1 — ji —@) f sin [7 +4 V(n)] in 
: (cae 


J 


In these expressions, the function (7) is to be obtained from (28). Thus our result gives the para- 
metric equation of the free ring in terms of four double integrals. 


y 


RAZ At =A) J cos [In — a V(n)] dy — + (1— 1 —@) { cos In + @ V(y)] dy. 
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Mr. Warren de Witt has kindly reduced the function Vv (y) to an expression in terms of elliptic 
integrals, as follows: 


en ets 1 — e? esiny + yl —e cos? y 
Ah cpemesewerer cet 2 ies y + Vl —e? cos? y 
V (py) yl e cos? p e siny + 16 in y Cos v) = log fine 
sess aa gk 
+ “Z* siny Are sin (¢ cosy) + 2=* Ely sm 6] + (Ln, o} 


Gly dade e(ton] e 


The function (yp), being independent of the load, may be tabulated once and for all for each 
. eccentricity e. 

From (31) it is plain that double series expansions for x/R and y/R in powers of e and a may 
be obtained and will converge for all relevant values of the parameters, i. e., when e< 1] and ais 
arbitrary. 


10. The exact solution for the circular ring. When the ellipse reduces to a circle, we have e = 0, 
and (24) becomes simply 


V(y) =yp+siny. (33) 
Substitution in (24) yields 
” 
x aoe! — {sin [(1— a) yn — asin 7] dy =1—G(a, py), 
; (34) 


y 
z= { 8 1 —4) 4 —asiny] y= Ca, 9). 
0 


equivalent to Meldahl’s Kq. (8). The functions G(a, y) and C(a, y) we shall employ again in § 14 
and § 15. . 


11. The slightly deflected elliptical ring. When the load factor a is very small, or, equivalently, 
when deviations from the elliptical shape are insignificant, we may replace the exact solution (31) 
by the leading term in its expansion in powers of a; equivalently, we may put (27) and (29) into 
(15). The result is 


we 2 
~ = cosy +ay1—@ { ¥(n) cosy dy, x = /1—@siny +4 f V(q) sing dy. (35) 
f 6 


The function (7) is not simplified in any way. 


12. The elliptical ring of slight eccentricity. When the eccentricity e is small, the function 
V (y), defined by (28) in general, may be approximated as follows: 


vo fd + cos é) (1 —esint > &) (1 — get eos?) ds 


yp : 
~ fa + cos &) eae |i + e088) dé =(1— 7 4) (y +-siny) + fe sin? yp. (36) 
0 


Our exact solution (31) is now sufficiently approximated by 
~ 


x ce | (1 : e) { sin [(1 — m) 7 — msin 7 — nsin® 9] dy 
6 


y 


— ef sin [1 + a) + asin ny] dy, 
4 
: (37) 
v 
Z = (14a) f cos (1 —m) 1 —msiny —n int] dy 


0 
Boer 
—7e { cos [(l-+ a) + asin 4] dy, 


0 
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where 


1 
m=a(l—74), n= ad. (38) 


A series expansion for these integrals in terms of Bessel functions, similar to that obtained by 
Meldahl for the case of the circle, is evidently possible, but the result will be rather elaborate, and 
numerical evaluation of the quadratures seems to be more practical. 


13. Quadratic theory for the slightly elliptical ring. In the quadratic approximation for the 
ellipse of small eccentricity, we neglect not only terms of order in excess of e? but also terms of 
order a2 e2. The resulting approximate equations for the free ring are 
__ x l CREAR | ren nee roe 
=F =z 4/ 3 rd braey 


i a(1 e* + : a aoe yra)cos yp +a(1—e + a) py siny —— ap sin 2 y 


1 1 11 i 1 
tae 1 za + qe] cos 2p + 5z a? 608 Sy + 55 eens Ap 


(39) 


1 : 1 3 1 
he (1—z4)siny = 5 a(1 rejyta(l—zé +a 


ces a ae 
gre) sing 


—ay (1—7e +a) cosy— 5a (I yet ya)sin2y 


a e 


1 1 : : 
+ 7 @ycos2y + 5 a’ sin 3 y + jo9 sin 4. 


The accuracy of these formulae for the case e = 0 was established in § 8. The terms arising from 
the eccentricity are in most cases of opposite sign from those arising from the second power of the 
load factor, suggesting that in the case of slight eccentricity these formulae are even more accurate 
than for the circular case. 


The free gap fc is given by 


1 f 


Bake Ay 11 
2 RR R 


3 f 5 
= yan(l +5¢—H4)- (40) 


pan 


Since the maximum error occurs at the free gap, we get from (40) the following estimate of the 
accuracy of the linearized theory: 


; 5 a\i 
Maximum percent error ~ (z eae ?) Fe (41) 


[Caution: while this estimate is considerably sharper than (21), which was derived from more 
general considerations, it has been obtained by regarding the quadratic theory as correct. Thus 
it is not valid for large a, while (21) is valid (for the circular ring) without any restriction whatever. | 
In general e? will be considerably less than a; then increase in eccentricity for an oblong ring with 
fixed load ratio increases the accuracy of the linearized theory. For example, if a = 0.03 ande = 0, 
ny maximum error is 21/,%, while if a = 0.03 and e = 0.10, the maximum error is reduced to about 
. 2%: ; : 


In a practical situation a is not known, but f¢/R is. From (40) we have, with considerable 


accuracy, 
2 fe 12s 
cele ly aos 0 NIR Ret 
3aR een te te ©: : a) (42) 
= (sac 34 ° 


As an example, we select a ring for which + f,/R = 0.15, e = 0.10. Then the load ratio accor- 
ding to the linearized theory is a = (10 2) = 0.03183, while the more accurate formula (42) 
yields a value about 2% smaller: a = 0.3116. 
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14. A nearly elliptical ring for which there is a simple but very accurate solution, valid for all loads. 
We return now to the results of § 12, which give the exact solution for an ellipse of small eccen- 


tricity. In (36) the eccentricity correction is proportional to y+siny Loge sin? y. The percent 
error made by neglecting the last term is ; 


Leahy 
3 ant y 
EQ Soe (43) 
y + sin y — ea sin’ y 


This function vanishes both at the foot of the ring and at the free gap, and its value nowhere 
exceeds 20%. When we recall that this is the percent error in a very small correction term, we 
see that no very large error, when e is small, will be made by using 


ry) = (1 x #) (y + sin y) (44) 


in place of (36). Equivalently, we shall consider a nearly elliptical ring for which (44) yields the 


solution. 


To get the result of using (44) in (31), we need only put n = 0 in (37). Comparison with (34) 
yields 


x 


ot (1 4) S(m.y) 4 &S(—a, W)s * =(1 =e) C(m, yw) = e?&(—a,y). (45) 


The result (45) expresses the free shape approximately in terms of the functions © and € occur- 
ring in the exact theory of the circular ring. It shows us how to adjust and combine the results 
for circular rings in order to get those for nearly elliptical rings. The load ratios are altered, and 
negative values must be employed. Now negative load ratios occur in the problem of a ring sub- 
jected to internal rather than external pressure. Thus the result (45) indicates that the nearly 
elliptical ring may be obtained by properly superposing two circular ones, one subject to external 
and the other to internal pressure.. 


The advantage of the present solution is its simplicity. Let a table of solutions for the circular 
ring, for various load ratios, both positive and negative, be constructed once and for all. Then (45) 
enables us to use these same tables to get the free shape of elliptical rings by simple numerical 
substitutions. 


15. Numerical specimen. So as to compare the accuracies of the various formulae, Mr. Harrison 
Hancock has calculated values for the elliptical ring in which e = 0.085 and a = 0.03. The exact 
free gap is given by + f,/R = 0.144251. The free shape of this ring was calculated according to 
the following four formulae: 


1. The exact solution (31). 

2. The approximation (37) for small eccentricity. 
3. The further approximation (39). 

4. The semi-conjectured approximation (45). 


In Table 1, the entries are the values of the rectangular co-ordinate to radius ratios, x/R and 
y/R, for the free ring, calculated at 15° intervals of the parameter y. The columns headed ,, Error“ 
were calculated from results to seven decimal places, not shown in the table. 


The quadratic approximation, given by (39), is seen to be not sufficiently accurate. For example, 
it gives a ratio fo/R in error by nearly 13/100%. The approximation (37) is accurate to within 8 
_ millionths of a radius; this is more accurate than is necessary. The semi-conjectured approximation 
(45) is accurate to within 40 millionths of a radius and yields a ratio f¢/R in error by less than 
3/100%. 

Thus, if the numerical example is a typical one, the semi-conjectured approximation (45) is 
adequate for current needs of design. It is to be emphasized that this formula yields a speedy 
solution for any elliptical ring of slight eccentricity in terms of a sufficiently extensive set of values 
of the two functions occurring in the exact solution for the circular ring. 
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Part III. Theory of nearly circular rings 
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Table 1. 
Ba Sy) ee 
Exac xima Error | ximate Error Approximate Error 
" (degrees) sii (31) eto (37) s 108 | eoitecns(9) x 10° goletiant (45) x 108 
a  —— 
x Ri 0 1 if 0 1. 0.0 if 0 
15 0.967 929 0.967 929 0.1 0.967 929 0.1 0.967 929 0 
30 | 0.873 587 0.873 587 0.1 0.873 588 0.8 0.873 586 — 0.5 
45 0.722 518 0.722 518 0.1 0.722 521 3.5 0.722 516 — 2.3 
60 0.523 708 0.523 709 Osby 0:523°717 8.8 0.523 703 — 5.2 
715 0.289 194 0.289 194 0.2 0.289 210 16.0 0.289 185 — 8.4 
90 0.033 472 0.033 472 0.1 0.033 495 22.6 0.033 462 —10.2 
105 | —0.227 276 —0.227 276 0.3 —0.227 251 25.5 | —0.227 286 — 9.6 
120 | —0.476 162 —0.476 162 0.2 —0.476 141 21.2 | —0.476 169 — 7.5 
135.) —0.696 702 —0.696 702 0.2 —0.696 694 8.2 | —0.696 707 — 5.0 
150 | —0.874 008 —0.874 008 0.2 | —0.874 022 —13.6 | —0.874 011 — 3.6 
165 | —0.995927 | —0.995 927 0.2 —0.995 970 —A2.8 | —0.995 930 — 3.0 
180 | —1.054027 | —1.054 027 0.2 —1].054 105 —77.0 | —1.054 030 — 2.9 
| 
y/R = 0 0 0 0 | 0 0 0 0 
15 0.258 226 0.258228 | 1.6— 0.258 227 1.5 0.258 228 1.6 
30 0.500 860 0.500 863. | 3.3 | 0.500 862 2.7 0.500 863 3.0 
AS 0.713 126 0.713 131 4.8 | 0.213) 1305 a 4.1 0.713 129 3.2 
60 0.881 878 0.881 884 5.9 | 0.881 886 | 1.7 0.881 879 a 
715 0.996 379 0.996 386 6.6 | 0.996 395 | 16.4 0.996 374 — 5.1 
90 1.049 031 1.049 039 7.9 1.049 063 32.7 1.049 018 —12.9 
105 1.035 988 1.035 996 8.0 | 1.036 044 56.4 1.035 966 —21.7 
120 0.957 589 0.957 596 1.5 0.957 674 85.6 0.957 559 —29.4 
135 0.818 532 0.818 538 6.6 0.818 649 VTA 0.818 497 —34.3 
150 0.627 728 0.627 734 D3 0.627 875 146.8 0.627 691 —37.1 
165 0.397 803 0.397 807 3.8 0.397 974 170.9 0.397 764 —38.9 
180 0.144 251 0.144 253 2.1 0.144 438 187.3 0.144 210 —40.6 


16. General solution for nearly circular rings. In the general theory developed in § 6, the ring 


has been assumed to have one axis of symmetry but not necessarily two. The fact that the crank is 
inclined to the axis of the cylinder, or special conditions of installation, may make a somewhat 
ovoid form more desirable than the elliptical form, having two axes of symmetry, which was 
subjected to study in Part II. 

If the numerical results given in § 15 may be taken as typical, we are not justified, in general, 
in adopting even a quadratic approximation in powers of the load ratio, but we are justified, for 
nearly circular shapes, in using a solution slightly perturbed from the one for the circle. Such a 
solution, for a general oval form, we now proceed to derive from the exact formula (12). 


Let an equation for the upper half of the loaded ring, in polar co-ordinates r, ®, be 
r=R[l+ef(®)], 0x @G<xz, 
where ¢ is a small parameter defining the eccentricity of the ring. The function f(®) we shall regard 


as given and fixed; there is no loss in generality in supposing that f(0) = f(a) = 0. We shall expand 
the exact solution (12) for this shape in powers of ¢ and discard all squares and higher powers. 


(46) 


For the chord length L we have 
L? = R? + r?— 2 Rrcos (x — ®) 
= R4+R(1+2¢ef)+2R(1+ef)cos® =4 R(1 + ef)eots ®, 
Then, since dS = R(1 + ef) dD, we have 


(47) 
Ss ® © 

cI=c [ 1*ds =e TR(L +ef)db =2a f [1 + 2ef(@)] co + Edé =al¥ (0) +e Q(®)], (48) 

where = ; 


®D 
To(P)=OH+sind, — MO) =A | f(E) cow® +E dé. (49) 
0 
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First putting the exact solution (12) into the form 
Zz 
a= Zeit tic [ etZ1ds (50) 
Zo 
and noting that Z = re'® = R(1 + ef) e'®, by use of (47) and (48) we obtain 
oD 
2 . . es . . 
Arse (1 + ef) e(—4 Vo e—iee2 1 Dig f pet ye —1082 (Le ef) cos? D d@ 


0 
PD 


= (@—2 Vd) [1+ ¢(f—-iaQ)] +2ia J ei(P—2 Vo) [1 +6 (3 f—iaQ)] cost > Dd. (51) 
i 


We now express this result in terms of deviation from the exact solution 2° = 2°(a, ®) = x® + iy 
for the circle of radius R, given by (34). Setting 


Ge is bx by (52) 
we write (51) in the form 
Se = oil V0 (f—iaQ) +2 a { eer (3if +42) cos? — Dd. (53) 
Hence 
ot = f(®) cos (P— a V(P)) + a QP) sin (b — a ¥()) 


| 
+2a f [a 22(€) cos (§ — a V,(€)) — 3 f(€) sin (§ —a V,(€))] cos? Edé, | 
and ; L (64) 
Y= f(®) sin (® — a V,(D)) — a QO) cos (6 —a V (9) 
+ 2f(®) J [3 f(&) cos (é = V o(é)) + a Q(€) sin (é =a V o(¢))] cos? > Edé. 


These results express the deviation from the exact solution for the circle of radius R. They are 
similar in type to that obtained in § 12 for the elliptical ring 1 in that only the eccentricity parameter | 
is supposed small, while the load ratio may be arbitrarily large. Thus, probably, the formulae (54) 
are more accurate than is necessary and may be simplified by approximations such as those leading 
to the semi-conjectured result (45) for the elliptical ring. 


(Eingegangen am 17, Marz 1960.) 
Anschrift des Verfassers: Professor C. Truesdell, Bloomington, Indiana (USA), 801 North College Avenue. 


1 1 ; ; ; if ee 
1 To compare the results, put ¢ = ag 2, O= y — a esin2y, f= sin? yp = sin? + 7 e? sin® 2 @. 
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Berechnung der Strémung des Tandemgitters 
mit bewegter zweiter Schaufelreihe 
(2. Hauptaufgabe) 


Von E.-G. Feindt 


1. Einleitung. In letzter Zeit sind mehrere Arbeiten erschienen, in denen die Durchstromung 
des aus zwei unbeweglichen Schaufelgittern bestehenden stationéren Tandemgitters behandelt 
wird! 3% 4, Fiir die Praxisist jedoch das instationare Tandemgitter besonders wichtig, bei dem die 
erste Schaufelreihe feststeht und die zweite mit konstanter Geschwindigkeit U in Gitterfrontrich- 
tung bewegt wird entsprechend der Leitrad-Laufrad-Anordnungen der Turbomaschinen (Abb. 1). 


Leitrad Laufrad 
Abb. 1. Das Tandemgitter (¢-Ebene). 


Eine genaherte Berechnung der inkompressiblen Durchstrémung des instationaren Tandemgitters 
wurde kiirzlich von O. S. Miiftiioglu' durchgefihrt. In der vorliegenden Arbeit wird ein Verfahren 
angegeben, nach dem die exakte Lésung fiir diese Strémungsaufgabe berechnet werden kann. Das 
Verfahren wird am Beispiel des Tandem-Plattengitters erprobt. 


2. Bezeichnungen. Es sollen bedeuten 

1 die Profiltiefe, 

t die Teilung, 

s _ die Spaltweite, 

€ =&§+i%y komplexe Koordinaten der physikalischen Strémungsebene, 
2=x-+ty komplexe Koordinaten der doppelperiodischen Strémungsebene, 


w =u—ziv die komplexe Geschwindigkeit, 
U die Geschwindigkeit, mit der sich das zweite Gitter bewegt, 
WwW, die Normalkomponente der Geschwindigkeit, _ 
W, die Tangentialkomponente der Geschwindigkeit 
(w, und w, sind reelle GréBen), 
Q die Quellstarke, 
yt die Wirbelstarke, 
q die Quellverteilung. 


1 W. E. Spranglin, Flow through cascades in tandem. NACA TN 2393 (1951). 

2 K. Nickel, Ing.-Arch, 23 (1955), S. 102. 

8 P. F. Byrd and M. F. Huggins, Ing.-Arch., 21 (1953), S. 191. 

4 H. Ohashi, Ing.-Arch. 27 (1959), S. 201. 

5 S. Miiftiioglu, Stromungsuntersuchungen an Tandem-Schaufelgittern. Jahrbuch 1957 der WGL, S. 208. 
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Ferner sind 3, 03, 0, Thetafunktionen und 


mr gr 1. og _ dg 


dz’ aly Pea? en da: 


ihre Ableitungen. 


Endlich sind K(k) und K’(k) = K(y1 — k2) die vollstandigen elliptischen Normalintegrale erster 
Gattung mit dem Modul k. 


Die Indizes sollen folgende Bedeutung haben: 


a ein Profil des ersten Gitters, 

b ein Profil des zweiten Gitters, 

Q von Quellen Q induziert, - 
I’ von Wirbeln I induziert, 

q von Quellverteilung q induziert, 

1 weit vor dem Gitter, 

2 weit hinter dem Gitter. 


3. Das Rechenverfahren. In inkompressibler Strémung pflanzen sich die Stérungen mit unend- 
lich groBer Geschwindigkeit im Strémungsfeld fort. Daher wird bei der instationaren inkompressib- 
len Strémung das Geschwindigkeitsfeld in jedem Augenblick von den Randbedingungen bestimmt, 
die in dem Augenblick vorhanden sind. Die Vorgeschichte geht nur dann ein, wenn freie Wirbel 
in der Strémung sind, deren Lage im Strémungsfeld von der Vorgeschichte abhangt. Im folgenden 
wird angenommen, da sich in der Strémung keine freien Wirbel befinden. Die Strémung durch das 
Tandemgitter kann dann als Potentialstr6mung berechnet werden, die in jedem Punkt der Profil- 
kontur des ruhenden Gitters die Randbedingung w,, = 0 und in jedem Punkt der Profilkontur des 
bewegten Gitters die Randbedingung w, = U cos f erfillt (Abb. 1). Zur Lésung dieses Randwert- 
problems wird das Doppelgitter zunachst mit Hilfe der Funktion 


22 
¥ ze 
Cts 6 (1) 


auf eine €*-Ebene konform abgebildet. Unter der Voraussetzung, daB die Teilung t fiir beide Gitter 
gleich ist, erscheint in der ¢*-Ebene nur je eine Schaufel des ersten und des zweiten Gitters. Das 
AuBere dieser beiden Schaufeln wird nach bekannten Methoden! weiter abgebildet auf das Innere 
eines Rechtecks 4 BADCD (Abb. 2). Dabei gehen die Schaufelkonturen in die vertikalen Rechteck- 
seiten AA bzw. DD iiber. Durch die nacheinander ausgefiihrten konformen Abbildungen werde 
der unendlich Ferne Punkt € = —oo auf den Punkt z, = x, +1 y, der z-Ebene abgebildet. Ent- 
sprechend sei 2, = x, +1, das Bild des Punktes € = + co. Wie beim Einzelgitter kann man 
schlieBen, daB in z, eine Wirbelquelle 


Q,—itl, =(u,—in,)t 
und in z, eine Wirbelsenke 


Q,—il, =(— uy +i r)t 


liegt. Durch Spiegeln an den Rechteckseiten AA und DD sowie durch periodisches Fortsetzen 
iiber die beiden Seiten AD wird das Strémungsfeld des Rechtecks auf die volle z-Ebene erganzt. Es 
entsteht so eine doppelperiodische Anordnung von Wirbelquellen, die die Geschwindigkeit induziern 


gf) eo eee) 
ve 0. — ibe 4K Cpa e eek ae 
CE eeta eG 1K) Qn z+(4—iK’) 
( 4K Os AK 
a | | 
ip Omner co i Ona ot 4K : (2) 
2g (2— Gat K) amy (2+ @—iK) 
«( 4K | 4K 


1 4, Betz, Konforme Abbildung, Berlin/Géttingen/Heidelberg 1948. 
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Fiir das Geschwindigkeitsfeld (2) sind die Rechteckseiten 4A und DD Stromlinien, d. h. das 
Geschwindigkeitsfeld (2) entspricht dem stationéren Doppelgitter (U = 0). Die Strémung fiir das 
Doppelgitter mit bewegter zweiter Schaufelreihe ergibt sich, wenn gemaB Abb. 2 auf den Vertikalen 


x=K-+4nK (n=1, 2, 3,.%.) 
Quellverteilungen 
g = 2 w,(2) 
angebracht werden, wobei 


w,(z) = Ucosf | 


= ii 


unten 


Gril; | Yul; 


Abb. 2. Anordnung der Singularititen in der z-Ebene. Die spezielle Lage von 2, und 2, gilt fiir das Zahlenbeispiel. 


die auf die 2-Ebene umgerechnete Geschwindigkeitskomponente normal zur Schaufeloberflache des 
zweiten Gitters ist. Die Geschwindigkeiten, die die Quellverteilungen auf den Schaufeloberflachen 
induzieren, werden fiir die Schaufel DD aus der Gleichung 


ee Z jibe 4 4K 
walle finn | a) tent ae 
q y) Tas 9 (SD dy (3) 
4K 
und fiir die Schaufel 4A aus der Gleichung 
: 5 fe ORES Sip eee 
a ete Ooi TK = 
Ks 4 4K 


berechnet. Das Integral (3), dessen Integrand fiir y = y’ singular wird, ist als Cauchyscher Haupt- 
wert zu verstehen. 
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Zu den Geschwindigkeiten (2) und (3) bzw. (4) der z-Ebene kann man noch eine beliebige verti- 
kale Parallelstrémung w = —ic hinzufiigen. Dieser Parallelstrémung entspricht in der ¢-Ebene 
eine Zirkulationsstrémung mit entgegengesetzt gleich groBen Zirkulationen 


te 2K! biws b= 2 ¢ K’ 


der Schaufeln des ersten bzw. zweiten Gitters*. Hiermit erhalt man fiir die resultierende Geschwin- 
digkeit in der z-Ebene 


w(z) = We,r + wy—ic. (5) 
Im allgemeinen wird w, = u,—iv, und damit auch Q, = u,, Q, = — ug = — u, und J =v, 
gegeben sein. J, = — v, und c folgen dann aus der Joukowsky-Bedingung, nach der in den Punk- 


ten B und D der z-Ebene die Geschwindigkeit w verschwinden mu8. Wenn die Geschwindigkeits- 
verteilung w(z) auf den Rechteckseiten 44 und DD berechnet worden ist, folgt die gesuchte Ge- 
schwindigkeitsverteilung des Tandemgitters (¢-Ebene) aus der Gleichung 


= De ee | 
w(C) = w(2) Fe (6) 
und die Schaufelzirkulationen ergeben sich aus den Gleichungen 
K’ 
PS hh Im {w(— K +i y'} dy’, 
eed 
K’ (7) 


I,= f Im{w(K +iy’)} dy’. 
, BKC 

Offenbar ist 
i hg — Uae rs 


4. Zahlenbeispiel. Als Zahlenbeispiel wurde die Strémung fiir das in Abb. 3 dargestellte Doppel- 
plattengitter berechnet, dessen Geometrie durch folgende Daten festgelegt ist: Teilung t = 1 fiir 
beide Gitter, Tiefe der Platten / = 0,8 fiir beide Gitter, Spaltweite s = 0,2. 


10; 
7Y, U, 
Y 
Dgben 
A a, B 
Uy unten Uz 


Abb. 3. Das Tandemplattengitter des Zahlenbeispiels. 


3 Durch die konforme Abbildung (1) wird das Tandemplattengitter auf zwei einzelne Platten ab- 
- gebildet, die mit den Strecken 


0,00187 < €*< 0,28461 und 1< &* < 152,40954 
der £*-Achse zusammenfallen. Durch die lineare Abbildung** 


; __ 1,533.48 2* — 0,818 09 (8) 
: T = 0,46051C* + 0,24887 


- werden die beiden Platten auf die Strecken _ 
3 SST PAT ns eerie eee und l<t< 3,2643 


* Vel. S. 245 des in Fufnote 1 von S. 89 genannten Buches von Betz. 
** Vol. S. 231 des Buches von Betz. 
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der reellen Achse der t-Ebene abgebildet. Nachdem die beiden Platten in diese Normallage ge- 
bracht worden sind, ergibt sich fiir die anschlieBende konforme Abbildung auf ein Rechteck ein 
elliptisches Normalintegral erster Gattung 


T 


dt (9) 


mit dem Modul k = 1/3,2643 = 0,30635. Da die beiden Platten sowie die Bilder t = + 3,2872 der 
Punkte ¢ = ++ co in der t-Ebene auf der reellen Achse liegen, braucht das Integral (9) nur fiir 
reelle t ausgewertet zu werden, was mit Hilfe der bekannter Tafeln! geschah. 


-1 


ei 


Nee 


-08 —O4 0 O4 08 
Abb, 4. Auftragungen zum Zahlenbeispiel. 


In der Umgebung der Plattenendpunkte t = + 1 und t = + I/k reicht die Genauigkeit der 
Tafeln fiir die vorliegenden Zwecke nicht aus. Daher wurde der Integrand von (9) in diesen Be- 
reichen durch Naherungsfunktionen ersetzt, die geschlossen integrierbar sind. Die Naherungs- 
funktionen ergeben sich, wenn man in der Umgebung von t = 1 die Naherung 


Lac: Seoreeh ke ; 
jase fis ya —k23 =) 


und in der Umgebung von t = 1/k die Naherung 


pbscncens ee ena 1 1 
ren yet apt 
Vii ~@) 
in (9) einsetzt. 
Die halben Rechteckseiten K und K’ sind gegeben durch die vollstandigen elliptischen Normal- 
integrale K(k) = 1,60973 und K’(k) = 2,60800. Fiir die Bildpunkte z, und z, der unendlich fernen 
Punkte ¢ =— oo baw. € = + oo findet man z, = — 1,48624 + i K’ bzw. 2, = 1,48624 +7 K’. 
Die Quellverteilung gq wird aus der Gleichung 

= dt dt dt* dr 

= 2 w,(z) = 2 eMees 

PN ae WNC eee areca ag 


aioe mit Hilfe von (1), (8) und (9) leicht ausgewertet werden kann. In Abb. 4 ist die 
uellverteilung q zusammen mit den von ihr induzierten Geschwindigkeiten (3) und (4) aufeet 
Die Integrale (3) und (4) wurden numerisch ausgewertet. : eee 


(10) 


1 K. Hayashi, Fiinfstellige Funktionentafeln, Berlin 1930. 
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Die 3-Funktionen wurden aus den beiden ersten Gliedern ihrer Reihenentwicklung berechnet 
gema} den Formeln 


i Oy Zz 
(Fe) wah sin 1 ae —2 hs sin 3 ae ; 


s(x) * 1+2hcos2z ce 7 


ile ww 1—2hceos2x TK? 


POOR 2 
ond) mo Wes 


Cl ipa 
mit h=e ** =—0,078478. 


Wie beim ungestaffelten einfachen Plattengitter * lat sich auch in dem vorliegenden Beispiel 
des Doppel-Plattengitters eine wesentliche Vereinfachung dadurch erzielen, da man zunachst 
die Rechnung fir die gitternormalen Geschwindigkeitskomponenten u, = u, = 0 durchfiihrt und 
erst am Schluf der Rechnung eine Parallelstr6mung w = u, = uy tiberlagert. Es ist dann in (2) 
Q, = Q, = 0 zu setzen, und das Geschwindigkeitsfeld der z-Ebene ergibt sich als Uberlagerung der 
von den Wirbeln J), und J induzierten Geschwindigkeitsfelder Abb. 5a bzw. 5b mit demin Abb. 5¢ 


D 0 
4 
g 
“ C 
D D 
2k- 2k 

i ilstré — den Wirbeln I, induzierte Strémung; b) von den Wirbeln I’, induzierte Stré- 
Abb. 5. Skizzen der Teilstrémungen der z-Ebene. a) Von di eln I, in TEC SS ake ed oa alae 
: mung; c) von der Quellverteilung induzierte Strémung. — Zu ele ioe See itt als vi g 


* Vgl. das Buch von Betz. 
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skizzierten Geschwindigkeitsfeld, das von den Quellverteilungen induziert wird, und der Parallel- 
stromung w = —Uc. 

SchlieBlich sei noch erwahnt, daB sich eine einfache Rechenkontrolle daraus ergibt, daB die 
Geschwindigkeitsverteilungen fiir die folgenden beiden Str6mungen gleich sein missen: 


a) u=u,=—0, 71 =%=1, ,=f,=0, 
beide Schaufelreihen feststehend 
b)u,=u=0, y1=%=0, W=l,=0, 
beide Schaufelreihen mit der Geschwindigkeit U = — 1 bewegt. 


U4 


Ser eee 


10 y 


~O6 06 40 


Abb. 6. Geschwindigkeitsverteilung fiir die Schaufeloberseiten des Tandemplattengitters. 
Die Geschwindigkeiten der Schaufelunterseiten sind entgegengesetzt gleich den Geschwin- 
digkeiten der Schaufeloberseiten. — 1 Stationéres Tandemgitter (U = 0) mit gleicher 
Zirkulation 'g = Ip = —1 fiir beide Plattenreihen. Strémung weit vor bzw. weit hinter 
dem Gitter: u, =u, = 0, v, =—v, = 1.— 2 Stationares Tandemgitter (U=0) mit 
gleich groBer aber entgegengesetzter Zirkulation beider Plattenreihen, [g = 1, [p = —1. 
Strémung weit vor bzw. weit hinter dem Gitter: u, = u, = 0, v, = v, = 0. — 3 Instatio- 
nares Tandemgitter, U = 1, mit der Zirkulation Null fir beide Plattenreihen(['g = Ip = 0). 
Strémung weit vor bzw. weit hinter dem Gitter: u, = u, = v, = v, = 0. — 4 Strémungs- 
bedingungen wie bei Kurve 3 (also U = 1, u, = u, = v, = 0), jedoch mit glattem AbfluB 
an den Plattenhinterkanten, wodurch sich v, = 0,929 ergibt. 


Im Fall b) ergibt sich die Geschwindigkeitsverteilung in einfacher Weise durch sinngemaBe 
Addition der beiden Geschwindigkeitsverteilungen 3 von Abb. 6. Dies folgt aus der Symmetrie der 
beiden Schaufelreihen. 


5. Rechenergebnisse. Die Rechenergebnisse sind in Abb. 6 und 7 dargestellt. Abb. 6 zeigt die 
Geschwindigkeitsverteilung der Schaufeloberseiten fiir die in der Bildunterschrift angegebenen vier 
Falle. Besonders auffallend ist die geringe gegenseitige Beeinflussung der beiden Schaufelreihen 
bei der instationaren Strémung ohne Schaufelzirkulation, Kurve 3. Besonders interessant ist auch 
die starke Symmetrie der Geschwindigkeitsverteilungen J und 2. 

Unter der Voraussetzung, das die Joukowsky-Bedingung erfiillt ist, wird das stromende Medium 
von der bewegten Schaufelreihe mitgenommen, und man kann einen Mitfiihrungskoeffizienten 
(v, — v,)/U definieren. Im Fall des Beispiels 4 von Abb. 6 hat dieser den Wert (v,— v,)/U = v,/U 
= 0,929. Dieser Wert gilt offenbar auch noch nach Uberlagerung einer beliebig groBen gitternorma- 
len Geschwindigkeitskomponente u, = Ug. 
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Als Letztes ist in Abb. 7 der Einflu® der Umfangsgeschwindigkeit U auf die Geschwindigkeits- 


verteilung dargestellt. Fiir U = 0 stimmt die Geschwindigkeitsverteilung mit der nach Nickel* 
berechneten iiberein. 
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16 a \ =) T 1 
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Abb. 7. Einflu§ der Umfangsgeschwindigkeit U auf die Geschwindigkeitsverteilung des 
Tandemplattengitters mit Zustrémwinkel «1, = 30°. Die Rechnung ergibt die Abstrém- 
winkel: a, = 0,21° fiir U/|w,| = 0; 0, = 29,8° fiir U/|w,| = 0,5; x. = 42,1° fir U/|w,| = 1. 


Zum SchluB sei noch angemerkt, das das beschriebene Verfahren nach Fortlassen der kon- 
formen Abbildung (1) auch geeignet ist fiir die Berechnung der ebenen Str6émung um zwei be- 
liebige Kérper, die sich gegeneinander bewegen und sich dabei evtl. auch noch in be- 
liebiger Weise deformieren. In der Literatur scheint diese Aufgabe noch nicht gelést zu sein. 
Als Beispiel ware etwa der Klappenfliigel mit schnell ausschlagender Klappe zu nennen. 


(Eingegangen am 11. Marz 1960.) 


Anschrift des Verfassers: Dr.-Ing. Ernst-Giinther Feindt, Hamburg-Garstedt, Breslauer Str. 13 


* Siehe FuGnote 2 von Seite 88. 
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Turbulente Ausbreitung eines runden Heifluftstrahls in ruhender Au8enluft 
Erster Teil: Kernbereich 


Von W. Szablewski 


1. Einleitung. In einer kiirzlich erschienenen Arbeit! habe ich die turbulente Ausbreitung 
runder Heiftluftstrahlen in bewegter Luft theoretisch beschrieben. Der Fall ruhender AuSenluft 
war dabei in die genannte Untersuchung nicht mit einbezogen worden, da er gegeniiber dem Fall 
bewegter Aufenluft eine wesentlich gréBere numerische Genauigkeit in der rechnerischen Be- 
handlung erforderte. Aber auch sonst weist dieser Fall Besonderheiten auf, auf die wir noch zu 
sprechen kommen werden. 

Die Methode der Untersuchung besteht, wie in Arbeit III, in der Anwendung eines Fortsetzungs- 
verfahrens, bei dem wir, ausgehend von den Anfangsprofilen, das Strémungsfeld approximativ in 
Differenzenschritten aufbauen. Mit Hinweis auf Arbeit III werden wir uns im Folgenden in der 
Schilderung der Einzelheiten der Methode kurz fassen kénnen. 

Die vorliegende Untersuchung beschrankt sich, wie schon Arbeit III, auf den sog. Kern- 
bereich des Mischungsfeldes, der durch den von der Mischung noch nicht beriihrten Strahlkern 
gekennzeichnet ist. Nach den in Arbeit III gemachten. Erfahrungen begniigen wir uns bei der 
Erfassung des Kernbereichs mit zwei Fortsetzungsschritten, die als praktisch ausreichend anzu- 
sehen sind. 

Unter Vernachlassigung von Reibungseinfliissen der Diisenwandung werden die Profile der 
Geschwindigkeit und Temperatur in der Diisenmiindung als rechteckig angenommen. 


2. Gleichungssystem des Mischungsfeldes. Unter der Voraussetzung konstanten Drucks lauten 
die Gleichungen? des rotationssymmetrischen Mischungsfeldes (Abb. 1), das als im zeitlichen 
Mittel stationér angenommen wird: 


Impuls: 
érom , arpud _ ou ~ 8 
aes = = ola) (ro + Bra), (1) 
Energie: 
or @ cp Tu ér@ cp Tv 00 cp T 
MBM a icin ey or (2) 
Masse: 
orou , Or@v _ o [. a 
a ee = Bele) (Fr): (3) 


(Die Querstriche, die zeitliche Mittelwerthildung bezeichnen, werden im Folgenden fortgelassen.) |} 
Es bedeuten dabei x, r Zylinderkoordinaten und u, v die entsprechenden Geschwindigkeitskompo- | 


Abb. 1. Geschwindigkeitsfeld. 


nenten, Q die Dichte, T die absolute Temperatur, c, die spezifische Warme bei konstantem Druck. 
E ist eine empitische Zahl und ein MaB dafiir, da® den substantiellen Kigenschaften eines Fliissig- 
keitselements eine andere Lange des Mischungsweges fiir den turbulenten Austausch zuzuschreiben 


1 W. Szablewski, Ing.-Arch. 26 (1958) S. 358; zitiert als Arbeit III 
® W. Szablewski, Ing.-Arch. 25 (1957) S. 10; zitiert als Arbeit IT. 
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ist als der Druckschwankungen wahrend des Transports ausgesetzten Geschwindigkeit; wie in den 
Arbeiten II und III setzen wir hier 


Fi 2 (4) 
Die scheinbare kinematische Zahigkeit der Strémung wird gleich 


E(x) = Hy by(x) Uy (5) 
gesetzt, wo u, die Geschwindigkeit des Diisenstrahls in der Miindung, b,(x) die Mischbreite des 


Geschwindigkeitsfeldes und x, einen empirischen Koeffizienten bedeuten. Mit T, bzw. T, werden 


des weiteren die Temperatur des Diisenstrahls in der Miindung bzw. die Temperatur der Auf®enluft 
bezeichnet. 


Die GréBe c, wird als Konstante behandelt. Mittels der Zustandsgleichung fiir ideale Gase 
bei konstantem Druck 


o T = konst. (6) 
ersetzen wir in den Gleichungen g durch 3 (= T — T,): 


konst. 1 
SRS Cece 7) 


Es werden weiter die dimensionslosen Funktionen g = u/uy und ¥ = 8/0), wo 0) = Ty) — T, 
ist, eingefiihrt und die Langen mit dem Diisenradius r, dimensionslos gemacht. Transformation 
der Gleichungen (1) bis (3) auf die Strahlenkoordinate 


L pia § 


pS (8) 


x 


ergibt dann, wenn wir noch mit dem Streckungsfaktor 


o = (4 =) (9) 


ca | x 


_ die Variablen & = (1/c) (x/rp) und y = o 7 einfiihren und y fiir o v/uy schreiben, schlieBlich folgendes 


Gleichungssystem (vgl. Arbeit ITI): 
Impuls: 


ap | oy (Fo/T) Ox/éy E ap 


Energie: 


Oe On. 1.2 (o/T1) Ox/ay g 2 OF 
yt oy [EY —W— 214 @ytyy IF Ey|— EPPA” ” 


Masse: 


op oy g at 
Soe ty yt Pitty te 


Im Vergleich mit dem Feld konstanter Dichte treten die unterstrichenen Glieder zusatzlich auf. 


Dabei haben wir (vgl. Arbeit IIT) vorausgesetzt, daB fiir das einzelne Modell, gekennzeichnet 


durch den Geschwindigkeitsparameter (uy — u,)/Uy (u, Geschwindigkeit der AuBenluft) und den 


Temperaturparameter ?,/T;, 
o(x) = konst. (13) 


ist. Man kann diese Voraussetzung auch so lesen, dafi fiir den Kernbereich x, und 6,/x als Kon- 
stanten angesehen werden kénnen. 


Die Randbedingungen lauten im Kernbereich, wenn man beachtet, daB die Strahlachse die 
Koordinate y = — 1/€ hat, 


Ll ¢% — Ue 14 
P46 fiir y>| ne (14) 


Weiter ist Verschwinden der Querkomponente auf der Strahlachse zu fordern: 


y >0 fir y—>—l/é. 7 (15) 


9) il 


98 W. Szablewski: Turbulente Ausbreitung eines runden Heifluftstrahls in AuBenluft —— Ingenieur-Archiv 


Als Anfangsbedingung ergeben sich fiir > 0 aus den Gleichungen (1) bis (3) Profile der Geschwin- 
digkeit q, y) und Temperatur 7, die wir bereits in Arbeit II berechnet haben und die, wie wir 
anmerken wollen, die turbulente Auflésung des Randes eines ebenen Heifluftstrahls beschreiben: 


P(E. ¥) > Pol) 
pl, ¥) > poly) fir >0. (16) 


x(&.¥) > Xo) 


Gema®B der Annahme (13) o(x) = konst. kénnen wir o(x) der Arbeit II [o = o(0)] entnehmen. 
Dort hatten wir die Breite o c, der Geschwindigkeitsverteilung im y-Mafstab zwischen den Gren- 
zen 0,95 und 0,05 der berechneten Verteilung abgelesen. Der Vergleich mit den Messungen von 
O. Pabst! ergab dann fiir den empirischen Koeffizienten x, den Zahlenwert 0,0082. 


3. Fortsetzungsverfahren. Wir berechnen im Folgenden die Geschwindigkeits- und Temperatur- 
felder bei Ausbreitung des Diisenstrahls in ruhender Aufenluft fir die Temperaturparameter 


Oo/ Ty — 0, Ue De 
a) Erster Fortsetzungsschritt. Differentiation der Gleichungen (10) bis (12) ergibt (vgl. 
Arbeit III), fiir die bendtigten Ableitungen (@p/0£):~o baw. (07/0&)z=0 baw. (Ay/0E)e—0, be- 

zeichnet mit ®, bzw. X, baw. WY, die Gleichungen 
®, + P(E | +96 


7 
0 


8/1) Xo \’ 
2(Pyy—%) —E+D (TE GiT ym) t1]—2%%=9 0) 


mit den aus (14) folgenden Randbedingungen 


d, +0 fiir y> ae ; (18) 
weiter 
xX” tis Xe , {2 ey yo (Io T) Xo : bee 2 ais 
0 + Xo x 4 Xo E (Dy Y) 2 (GIT) i care Ev X0 = (19) 
mit den (14) gemaBen Randbedingungen 
xX, 0 fiir yal, (20) 
CO 
schlieBlich, eingedenk (15), 
y y 
By =— S ody +yPy—2 JDy dy. (21) 


Die in den Gleichungen auftretenden Funktionen qo, yo und 7%, 7) Sowie Yy wurden, wie erwahnt, 
bereits in Arbeit II berechnet. 


Das Gleichungssystem (17) bis (21) haben wir mittels des in Arbeit IIT beschriebenen Iterations- 
formalismus gelést. Die nach zahlreichen Iterationen erhaltenen Verteilungen fiir ©), X, sind in 
Abb. 2 dargestellt. Die Verteilungen zeigen, wie in Arbeit III, eine mit dem Pama Ogle 
gehende Linksversetzung der Kurven, die dem analogen Verhalten der Anfangsprofile (Atbert II) 
entspricht. Auch die mit 3)/T, wachsende Amplitude entsprechend der mit gréBer werdenden 
Temperaturdifferenz zu erwartenden starkeren Kriimmung der Isotachen und Isothermen zur 
Achse hin tritt wieder in Erscheinung. Als Besonderheit gegeniiber den sich im Fall bewegter 
AuBenluft ergebenden Verteilungen ist jedoch zu konstatieren, daB hier unsere Verteilungen fir 
groBe positive y vor dem Abklingen gegen Null zu positiven Werten iibergreifen, was einer Ab- 
kriimmung der 4uSeren Isotachen und Isothermen von der Achse weg bzw. einem konkaven Ver- 
lauf dieser Kurven entspricht. Diese dem Fall ruhender AuBenluft zakommende Erscheinung ist 
eine Higentiimlichkeit der zugrundegelegten Theorie, wie ein Vergleich mit dem Voskalten det 
analogen Kurven fiir den Fall (uy —u,)/uy = 0,95, also fast ruhender AuBenluft, zeigt (Abb. 3). 


<a sieht, daf} auch hier diese Eigentiimlichkeit, wenn auch in stark reduzierter Form, festzu- 
stellen ist. 


1 0. Pabst, U. u. M. 8004 (1944). 
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Ausgehend von den Anfangsprofilen qo, y) ist dann der erste Fortsetzungsschritt durchfiihrbar 
und liefert mit der Schrittweite a approximativ die Profile 


%1 =P) +aQ,, Wi =%o +aXy. (22) 


In der vorliegenden Rechnung wurde die Schrittweite a — 0,2 genommen. 


y 


3. 


ate es See te) 
8 20 a 


Abb. 2. 


Abb. 3. 


b) Zweiter Fortsetzungsschritt. Nachdem wir oben (7, und y, ermittelt haben, berechnen 
wir des weiteren zundchst das Profily, der Querkomponente der Geschwindigkeit an der Stelle 


& = a gem4S der Formel (vgl. Arbeit ITT) 
y 


(l+ay) ln + PY |—@+ 1) 
1 Pi 


(Po/ Ts) x4 w I 
— d 
esphaek aeO ic F Fi 


1 + (9/ Tr) 44 1 +ay 


(23) 


Damit ist erst einmal der erste Schritt abgeschlossen. 
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Die fiir den zweiten Schritt benétigten Funktionen (0p/0£):=_ baw. (04/0&):—., bezeichnet 
mit ©, bzw. X,, berechnen wir dann (vgl. Arbeit IIT) aus den expliziten Formeln 


eeee Cat (Ole a | 
= Py “is Py 2 (Py V1) , fe | » 1 + (9o/Ty) Xy 5 l oe ay 24 
D, = 7S 3 i is ee ) ( ) 
aQy 
und 
f Flee (9o/T 1) x3 a 
’ : pA = = Sr a 
x ss ane x1 p (ay 1) it ae (9o/T) % ‘i ale ay (25) 
1 (2/E) a9, 
Hierbei ist nach (22) 
P= +taH, . ee ee (26) 
M1 =o + 4X; M1 =X +4Xo- 


Man bemerkt, da in den Formeln (23) bis (25) fir y—>oo nicht nur die Zahler, sondern auch, 
im Gegensatz zu dem Fall bewegter AuBenluft, die Nenner Null werden. Das erklart schon, warum 
der Fall ruhender AuBenluft eine besonders groBe Rechengenauigkeit erfordert. 

Die berechneten Verteilungen ®,, X, haben wir in Abb. 4 dargestellt. Im Zusammenhang 
mit der obigen Bemerkung muf dabei auf Folgendes hingewiesen werden: Die Rechnungen ergaben 
nur fiir den in Abb. 4 wiedergegebenen y-Bereich einen normalen Verlauf der Verteilungen; fiir 
groBere positive y-Werte — untere Grenze derselben in Abb. 4 gestrichelt markiert —- ergab die 
Rechnung einen sehr unregelmaBigen (und hier nicht dargestellten) Verlauf der Kurven. Es zeigte 
sich damit, dafs fiir diese gréBeren Werte von y die von uns mit Tischrechenmaschinen erzielte 
Rechengenauigkeit nicht ausreichend war, um die in den Formeln (23) bis (25) fiir y +0oo auf- 
tretende Stelle der Unbestimmtheit zu bewaltigen. Doch betrifft dieser von uns derart nicht 
erfaBte Bereich praktisch nur die 4uBeren Grenzen des Mischungsfeldes. 


Mit den so erhaltenen Funktionen ®, und X, kénnen wir dann den zweiten Fortsetzungsschritt 
durchfiihren und erbalien mit der Schrittweite b approximativ die Profile 


(ox Soe een 
was bX, =% +a4X,+5X,. 


Die Schrittweite 6 wurde zur Ermittlung des Geschwindigkeitsfeldes bis zum Kernende, das die 
Grenze der Anwendbarkeit von (5) und damit unseres Gleichungssystems bezeichnet, in den Gren- 
zen 0,15 bis 0,50 variiert (Tabelle 1). Wir gelangen dabei z. T. etwas iiber das Kernende hinaus 
doch ist fiir die hier vorgenommene geringe Uberschreitung des Kernendes (Isotache 0,95) ie 
Fehler als gering anzusehen. 

Der Abfall der Geschwindigkeit innerhalb der angezeigten Grenze ergibt sich dabei in folgender 
Weise: Tragt man die Kurven 7,(y) fiir die jeweilige Schrittweite b iiber y auf, so wird ihr Existenz- 
bereich zu negativen y hin durch die Koordinate y = — 1/§ der Strahlachse begrenzt. Hat man 
mit einem Schritt b das Kernende iiberschritten, so hat P2(y) an der Stelle y = — 1/& einen Wert 
kleiner als Eins. AuBerdem ist aus Griinden der Symmetrie zu verlangen, daB Y(y) an der Stelle 


(27) 
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shes 1s die Tangente Null hat. Daraus ergibt sich die Konstruktion des Abfalls der Geschwin- 
digkeit langs der Strahlachse (Abb. 5). Entsprechend ergibt sich innerhalb der angegebenen 
Grenze (y = 0,95) der Abfall der Temperatur langs der Strahlachse. 


4 Pstrarlachse 
ak 4 
a 499 
085|- a 
|- a 
E 1 
Lt ! 
-4 =o 


Abb. 5. Konstruktion des Abfalls der Geschwindigkeit langs Strahlachse. 


4. Feld der Isotachen und Isothermen im Kernbereich!. Wir beschreiben das Mischungsfeld 
im Kernbereich mittels des Verlaufs der Isotachen und Isothermen: r/r, iiber x/r, fiir y = konst. 
bzw. vy = konst. Die effektiv vorhandenen Geschwindigkeiten und Temperaturen ergeben sich zu 


T =(T) — T;) x + T;. 


waa y 3 


Tabelle 1. Isotachen. 


§ = (1/0) (x/r) 


1 1,55 1 1 1,54 | 1,81 2,09 
1 1,41 1,72 1,93 il 1 
1 1,19 1,34 1,43 1 1 
1 1,08 1,13 1,15 1 1 
1 0,99 0,95 0,92 1 1 
1 0,85 0,72 0,62 1 1 
aE 0,80 0,62 0,49 il 1 
1 0,69 0,43 0,20 1 1 


: 05 1 1 ied TE Pe Sal 3,01 1 2,29 ~- — 
1 il 1 1,72 | 2,34 2,81 i} 1,70} 2,29 2,89 
3 1 1,31 1,56 1,76 1 1,24 | 1,38 1,48 1 1202 29 1A 
Sy Ye 1 1,12 ate Lg 1 1,05 | 1,04 1,01 dh 1,02 | 0,99 0,94 
TT) 1 0,97 0,90 0,84 i 0,91 | 0,81 0,71 1 0,88 | 0,78 0,67 
9 1 0,78 0,577 | 0,42 1 0,74 | 0,53 0,34 1 O71) 0555 0,36 
95 1 0,70 0,43 0,22 1 0,66 | 0,40 0,16 1 0,64 | 0,44 0,21 
99 1 0,55 0,15 — 1 0,51 | 0,12 — 1 0,51 | 0,24 — 
95 


1 Die Durchfiihrung der langwierigen und schwierigen Rechnungen verdanke ich meiner Mitarbeiterin 
Frau Lilli Schulze. 
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Nach der in Abschn. 2 gegebenen Definition der Koordinaten YH Cl Diese Voy ar 
haben wir 


= oe; (28) 

T) we 1 aa Ey ? ro a€ ? 
gema8 der in Abschn. 2 getroffenen Annahme (13) o(x) = konst. fiir das einzelne Modell ent- 
nehmen wir o der Arbeit II (o = 0(0)): 


AP 
% 
=Q05 
a 
He tg ee 
7, GS Ug ? i = 
“Al = Q7 
G9 Z 
| ~ 285 595 = : % 
G99 y/ 
= 
LtLy ——— 
oe 03 
G2 
ot 
q 
zed 
% 
See 
B2 Ug-u. vB, 
(i hess } 
VitsAl = e Up ih Foe 
7, — 
7 
09 @ 
d ee 70 


Abb. 6. Feld der Isotachen und Isothermen, 


Tabelle 1 enthalt das Ergebnis unserer Rechnung. 
und Isothermenfelder fiir den Kernbereich dargestellt, 
peratureffekt eine longitudinale und laterale Verkiirzu 
entspricht der aus Arbeit IT folgenden Drehung des W. 


In Abb. 6 sind die erhaltenen Isotachen- 
Es ergibt sich, wie in Arbeit III, als Tem- 
ng des Mischungsfeldes; diese Verkiirzung 
inkelraumes der Mischung, mit Scheitel in 
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der Diisenkante, zur Strahlachse hin. Weiterhin ist, wie in Arbeit III, ersichtlich, daB das Tem- 
peraturfeld, bedingt durch den sich aus dem Experiment ergebenden Wert E = 2 fiir das Uber- 
tragungsverhaltnis [vgl. (4)], eine wesentlich andere Struktur als das Geschwindigkeitsfeld auf- 
weist, wie sich in seiner gréferen Breite und 
geringeren Tiefe zeigt. 


Tabelle 2 enthalt die Werte der Querkompo- 
nente y =o v/u, die in Abb. 7 dargestellt ist. 
Fiir die effektive Querkomponente der Geschwin- 


digkeit folgt v = Uy plo. 


Tabelle 2. 


0,4 1,54 2,29 2,75 
0 1,78 2,99 3,00 
0,4 1,43 2,19 2.63 
0,8 +0,56 1,40 1,88 
1,2 —0,55 +0,46 1,02 
1,6 1,62 —0,41 40,21 
2,0 —2,49 —1,14 —0,47 
2,4 —3,11 —1,69 —1,00 
2,8 —3,50 | —2,07 —1,39 
3,2 —3,71 —2,32 —1,67 
3,6 —3,79 == 9 AF Ei 85 
4,0 —3,79 —2,53 —1,97 
4,4 —3,75 —2,55 —2,03 
4,8 —3,70 —2,54 —2,05 
252 —3,66 —2,50 —2,04 
5,6 —3,63 —2,46 For Abb. 7. 


5. Vergleich mit Messungen. Zum Vergleich mit der Theorie haben wir einige neuere ameri- 
kanische Messungen herangezogen. 


S. Corrsin! hat das Mischungsfeld eines Diisenstrahls, der die kleine positive Temperatur- 
differenz von ungefahr 13°C (%)/T, ~ 0) gegeniiber der umgebenden Luft aufwies, gemessen. 
In Abb. 8 haben wir Experiment und Theorie fiir das Geschwindigkeits- und Temperaturfeld des 
Kernbereichs dargestellt. Fiir das Geschwindigkeitsfeld ist die Ubereinstimmung als gut zu be- 
zeichnen. Auch fiir das Temperaturfeld ist die Ubereinstimmung als befriedigend anzusehen bis 
auf den 4uBeren Teil, wo die duBeren Isothermen weiter nach auSen greifen. Dennoch ist die 
von der Theorie angezeigte verschiedene Struktur der Geschwindigkeits- und Temperaturfelder 
— gréBere Breite und geringere Tiefe des Temperaturfeldes — gut erkennbar. 


Auch die von S. Corrsin und M. Uberoi? durchgefiihrten Messungen an HeiSluftstrahlen mit 
Ubertemperaturen von 15° und 296° C (A/T, ~ 0 und 1) zeigen die geringere Tiefe des Tem- 
peraturfeldes gegeniiber dem Geschwindigkeitsfeld; insbesondere erweisen die Messungen auch 
den von der Theorie angesagten Temperatureffekt der longitudinalen Verkiirzung des Mischungs- 
feldes (sowohl fiir das Geschwindigkeits- wie das Temperaturfeld) mit wachsender Temperatur- 
differenz. Das zeigen in guter quantitativer Ubereinstimmung mit der Theorie die Messungen 
langs der Strahlachse (Abb. 9). Uber das Geschwindigkeits- und Temperaturfeld in seiner Gesamt- 
heit wurden fiir den Kernbereich in der Messung leider keine Angaben gemacht. 


1 §. Corrsin, NACA-Wartime Rep. W-94 (1946). 
2 §. Corrsin and M. S. Uberoi, NACA Rep. 998 (1950). 
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Abb, 9. Geschwindigkeit und Temperatur langs Strahlachse. 


6. Zusammenfassung. Es werden die Geschwindigkeits- und Temperaturfelder der turbulenten 
Ausbreitung von Heifluftstrahlen in ruhender Aufenluft fiir den Kernbereich des Diisenstrahls 
berechnet. Grundlage der Rechnung bildet das vom Verf. in einer vorangehenden Arbeit (Ing.- 
Arch. 20) aufgestellte Gleichungssystem, das mittels eines Fortsetzungsverfahrens gelést wird. 
Der Vergleich mit der Messung ergibt gute Ubereinstimmung. 


(Eingegangen am 21. Marz 1960.) 


Anschrift des Verfassers: Dr. W. Szablewski, Berlin W 8, Mohrenstr. 39. 
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Der schlanke Stab unter konservativer Belastung als Variationsproblem 


Von H. Leipholz 


1. Einleitung. In friiheren Arbeiten!? ist die Knickung sowohl des spannungsfrei verwundenen 
Stabes als auch der tordierten Welle, beide mit Einzelkraft und kontinuierlicher Langskraft be- 
lastet, behandelt worden. Im ersten Fall, fiir den Stab, waren zwei verschiedence, im zweiten Fall, 
fiir die Welle, zwei gleiche Haupttragheitsmomente vorausgesetzt worden. In der vorliegenden 
Arbeit sollen alle Einschrankungen fortfallen, so daB fiir Stab oder Welle, auBer der schon ge- 
nannten Druckbelastung, sowohl spannungsfreie Verwindung als auch Torsion und noch zwei 
verschiedene Haupttragheitsmomente zugelassen werden. 

Dieses verallgemeinerte Problem lat sich als Variationsproblem lésen, wenn man die tibrigen 
Voraussetzungen beibehalt, namlich 


1. die Stabachse sei urspriinglich gerade, 


2. die fiir die Stabquerschnitte beziiglich ihrer Schwerpunkte ermittelten Haupttragheits- 
momente seien iiber die Lange des Stabes konstant, 


3. die unter der Belastung auftretenden elastischen Verformungen seien klein, 


A, die Lagerung der Stabenden und die auf den Stab einwirkende Belastung seien derart, daB 
das ganze Problem konservativ wird. 


2. Die Differentialgleichung. Es sei ein Stab zu Grunde gelegt, dessen Beanspruchung im 
allgemeinsten Fall aus einem axialen Torsionsmoment T, aus einer Einzelkraft P und aus einer 
kontinuierlichen Langskraft q besteht. An den Stabenden greifen T und P an; q ist nach einem 
vorgeschriebenen Gesetz q(s) iiber die Stabachse verteilt, wobei s die Lange der Stabachse vom 
unteren Stabende an gezahlt ist, und P und q sind richtungstreue Druckkrafte, die die Richtung 
der unverformten Stabachse beibehalten (Abb. 1). 


VE 


| 
| 
i 
| 
| 


Abb. 1. Abb. 2. 


Ferner sei der Stab gema Fall 1 von Abb. 2 gelagert; nur wenn im besonderen das Torsions- 
moment T entfallt und lediglich die Druckkrafte P und q wirken, kénnen auch die Lagerfalle 2 
und 3 beriicksichtigt werden. 


1 H. Leipholz, Ing.-Arch. 29 (1960) S. 262 (kiinftig zitiert mit I.). 
2 H. Leipholz, Ing.-Arch. 30 (1961) S. 42 (kiinftig zitiert mit I1.). 
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Bei diesen Bedingungen fiir Belastung und Stablagerung ist das Problem nicht nur konser- 
vativ, sondern auch homogen, weil an den Stabenden keine Querkrafte auftreten oder diese.zu- 
mindest von héherer Ordnung klein sind. 

Wegen des konservativen Charakters kénnen zur Ableitung der Differentialgleichungen die 
Gleichgewichtsbedingungen der Elastostatik auf ein herausgeschnittenes Stabelement angewendet 
werden. 

Als Bezugssystem sei das Hauptachsensystem! (&,7,¢) des Stabes verwendet, das die Dreh- 
geschwindigkeit 1p hat, wenn man mit der Geschwindigkeit Eins an der Stabachse entlang geht. 
Dann fiihren die Gleichungen von Kirchhoff und Clebsch nach einigen Umformungen, wie sie ganz 
entsprechend in der Arbeit I durchgefiihrt worden sind, auf die Differentialgleichungen des ver- 
allzgemeinerten Problems 

a é, + [(a@ 4-6) (€ + t%) — ce] e+ [—b(e +a)? +e(e+m)¢+ P+ Q(s)] é, ae (1) 
b ¢, —[(a + b) (¢ + %) —ce]e, +[—a(e+%m)? +e(e+H)e+ P+Qs)]e=0. 
Hier sind eg und e, die Komponenten eines Einheitsvektors ¢ in der urspriinglichen Stabachse; 
die iibrigen Bezeichnungen in (1) werden alsbald noch genauer erklart werden. 

Man muf bei jenen Umformungen nur beachten, da$ man jetzt nicht nur den Vektor der 
elastischen Deformation u zu beriicksichtigen hat, sondern daB auch noch die Vorverwindung 
des Stabes mit p = Tt) ¢; zugelassen ist. Dabei ist t) die Verwindung der Langeneinheit der un- 
belasteten Stabachse und e, der Einheitsvektor in der ¢-Achse. Deshalb besteht zwischen den 
Komponenten der Vektoren tv, u und 0, weil die Auslenkung der Achse des belasteten Stabes 
von der Lotrechten klein sein soll, die Beziehung 


To i 0.) Cues OAL (2) 
Te E To 


Es sind, wie friiher, wg, w,, Tt; die Komponenten von {v, ferner wz und w, die Kriimmungen der 
Achse des deformierten Stabes und ¢ die zusatzliche Torsion der Langeneinheit des Stabes in- 
folge T. 

AufBerdem war fiir die Umformungen das Gleichungssystem 


M; = a,, M, =bw,, i; =e (3) 


zu benutzen. Hierbei sind M;, M, die Komponenten des Biegemomentes, T; das Torsionsmoment, 
a und b die konstanten Biegesteifigkeiten und c die konstante Torsionssteifigkeit des Stabes. 


Zu den Gleichungen (1) ist noch zu bemerken: 
1. Q(s) ist eine Abkiirzung fir fae ds , 
2. die unabhangige Variable, Wack der differentiiert wird, ist s. 
Ferner ist wieder, wie schon gezeigt? worden ist, 
ct; = T,; =konst. also T,=T 
und daher, wegen 1; == € + t) und c = konst., auch ¢ + 1) = konst. 


In dem Differentialgleichungssystem (1) sind einige bereits bekannte Probleme als Sonderfalle 
enthalten: 


LP. = 0, Q(s) =0, t) =0 ist das Problem der tordierten Welle, das sowohl fiir den Fall 
a = b (gleiche Haupttrigheitsmomente), als auch a b (verschiedene Haupttragheitsmomente) 
schon gelést worden ist’. 


2. Q(s) = 0 ist das Problem der tordierten und gleichzeitig gedriickten Welle (durch Einzel- 
kraft P an den Stabenden), das eingehend von H. Ziegler*, A. Trésch®, M. Beck® entweder fiir 


1 H. Leipholz, I., S. 263. 

2 H. Leipholz, I1., S. 45 

R. Grammel, Z. angew. Math. Mech. 3 (1923) S. 262. 

é are Ziegler, Z. angew. Math. Physik 2 (1951) S. 265 und 3 (1952) S. 96; Schweizer Bauzeitung 66 (1948) 
° A. Trésch, Ing.-Arch. 20 (1952) S. 258. 
6 M. Beck, Ing.-Arch. 23 (1955) S. 231. 
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a=b (gleiche Haupttrigheitsmomente) oder auch fiir a b (verschiedene Haupttragheits- 
-momente) dargestellt worden ist. Dabei haben die Autoren noch von den beiden Méglichkeiten 
T) = 0 oder t) ~ 0 Gebrauch gemacht. 


3. Q(s) =0, ¢ = 0 ist das Problem der »»Knickung verwundener Stabe unter Druck“, wie es 
durch die Arbeiten von H. Ziegler! und E. Hui? bekannt geworden ist. 


4. ¢€ = 0 ist das Problem der Knickung verwundener Stabe unter dem Druck einer konser- 
vativen, kontinuierlich und gleichmaBig verteilten Belastung, das ich schon friiher? behandelt 


habe. 


5. @ = b, t%) = 0 ist das Problem der tordierten und gleichzeitig gedriickten Welle mit gleichen 
Haupttragheitsmomenten, wobei die auf die Welle wirkenden Druckkrafte aus der Einzelkraft P 
und aus der richtungstreuen Langskraft q(s) bestehen; auch dieses Problem ist von mir? behandelt 
worden. 


Man ersieht daraus, da (1) das Problem der tordierten und gedriickten Welle (unter den 
bereits genannten Voraussetzungen) in sehr allgemeiner Form darstellt. Eine Lésung von (1) 
diirfte daher wertvoll sein, da sie die Lésungen aller Sonderfalle einschlieBen wiirde. Lésungs- 
méglichkeiten sind durch Reihenansatz oder durch Verwendung der Differenzenmethode gegeben. 
Das wurde von mir bereits in einer anderen Arbeit® angedeutet. In der jetzigen Arbeit soll 
gezeigt werden, da man die Lésung auch iiber ein Variationsproblem erhalten 
kann. 


3. Das Variationsproblem. Es werden die folgenden Abkirzungen eingefiihrt: 


a=—[a+N(e+m)—ce], Pa—rl-ble+n)P+ee(e+%)], 
(4) 
y =F lla +8) (6 +r) eel], 6 =fl-ale+n)§tee(e +m), 


a o 
old © 


«2 =rV2 (7 


und hieraus, wenn man fiir die linke Seite von (7) die Abkiirzung (5) verwendet, 


fae Oe (8) 


Jetzt wird noch eine Variablentransformation durchgefihrt: 


Man erkennt leicht, daB 


ist. Aus (6) ergibt sich 


b 
Mit (4) und (9) 14Bt sich (1) umschreiben zu 
ates (on P+0(s)\]/a — 
é + poe (a p 2400)., =0, Vee—r a, +(8 BeOS) Says 0: (10) 


Aus (10) erhalt man mit (5), (8) und durch Erganzung mit + u*e, in der oberen Zeile und 
mit -++ u? é; in der unteren Zeile 


é, + 2 ue; — ure, = o(s) &,, eg —2 we, — u? ee = y(s) ee, (11) 


1 Siehe FuBnote 4 von S. 106. 

2 KE. Hui, Osterr. Ing.-Arch. 9 (1955) S. 288. 
3 H, Leipholz, I. 

4 H. Leipholz, Il. 

5 H. Leipholz, 1. 
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wobei noch die neuen Abkiirzungen 


o(s) = —(u# +B+ ze: =) F y(s) = —(¥ L § + z ad (12) 


a 


benutzt worden sind. 

Das Differentialgleichungssystem ist in der Form (11) besonders geeignet, den weiteren Gang 
der Rechnung erkennen zu lassen. Zu diesem Zweck wird eine Analogiebetrachtung durchgefiihrt. 
Die Gleichungen (11) lassen sich namlich auffassen als die Bewegungsgleichungen eines Punktes, 
dargestellt in einem rotierenden Bezugssystem, wenn man nur die unabhangige Variable s, nach 
der differentiert wird, als ,,Zeit“‘ deutet. 

Ist és, e, das gleichformig mit der Drehgeschwindigkeit u rotierende orthogonale Achsenkreuz, 
so laBt sich auf Grund dieser Auffassung die Beziehung mit einem raumfesten orthogonalen Achsen- 
kreuz (x, y) vermége 

x= e,coSus—e, sinus, y = egsinus + @ cosus (13) 
und 

@: =xcosus+ysinus, e, =ycosus—xsinus . (14) 
herstellen. : 

Das System (11) lieBe sich formal, unter Verwendung der Transformationsformeln (14), in ein 
System fiir die neuen Variablen x und y umschreiben. Viel schneller und eleganter kommt man 
aber dazu, wenn man den Gedanken der kinetischen Analogie weiter verfolgt. 

Im Sinne der Kinetik miissen die Bewegungsgleichungen im raumfesten (x, y)-System lauten: 


%=K,, ¥=K,. (15) 


x 


Man hat jetzt nur die Aufgabe, die ,,Kraftkomponenten“ K, und K, aus den bekannten ,,Kraft- 
komponenten“ 


Kz = y(s) eg, K, =o(s) ¢,; (16) 
die man aus (11) abliest, zu errechnen. Das geschieht bekanntlich durch 
ge ee de Ses. de 
Ke A ik, reel Ky = Kye t as, - (17) 


Fihrt man die Rechnung gemaB (17) unter Verwendung von (14) durch, so erhalt man 
K,, = y(s) (x cos* us + y sin us cos us) —o(s) (y Sin us cos us — x sin® us) , (18) 
K, = y(s) (« sin us cos us + y sin? us) + o(s) (y cos? us — x sin us cos us). 


Durch Verbindung von (15) und (18) kann man das Ziel erreichen, das Differentialgleichungs- 
system (11) in den neuen Variablen x und y zu schreiben. Fiihrt man noch trigonometrische Um- 
formungen durch, so lauten sie 


zaxtt* tx" cos2us +y"“*sin2us, 
jay Pte y ">" cos2us +x "sin2us. a 

Ebenfalls durch die kinetische Analogie wird man zu den Begriffen 
T= (#@ +5), U =—[y(s) (x cos ws + y sin u s)? + o(s) (y cos us — x sin us)?] (20) 


gefiihrt, die fiir den bewegten Punkt die kinetische und die potentielle Energie bedeuten. Man 
sieht, da tatsachlich 


ont ii d aT 
K, ox? "ds dx ° 
und 21 
Pohet aor (21) 
ef ay’ ds ay 


ist, so daB man fiir (15) schreiben kann 


s\n, gen) 


‘ds\ dx Ay ds (22) 


| 
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Da aber T von x und y und U von x und y frei ist, so kann man nach Einfihrung der Funktion 
lige : = 
La = T—U oder La = ~ (** + y? + y(s) eZ + o(s) e?) (23) 
(Lagrangesche Funktion im Sinne der Kinetik) das Differentialgleichungssystem (15) in 


d /éLa La 0 d /éLa oLa 
ds\ ax ane ae) ay 


=0 (24) 


umwandeln. Dabei sind die GréBen ee und e,, die in (23) vorkommen, als Abkiir- 
zungen fiir Beziehungen zwischen x, y und u gem4B (14) aufzufassen. 


Mit (24) ist man aber zum Ziel gekommen. Denn es sind dies ja die Eulerschen Differential- 
gleichungen eines Variationsproblems: 


l 1 
aes == as) Extr, oderalso —f (x* + y? + y(s) ez + o(s) e?) ds = Extr.. (25) 
b 


Jetzt ist man zur Behandlung des Problems der gedriickten und tordierten Welle, das durch 
das System (1) in der allgemeinsten Form gegeben war, nicht mehr auf die sicher recht schwierige 
Auflésung des Differentialgleichungssystems angewiesen, sondern kann das Variationsproblem mit 
einer der direkten Methoden, z. B. mit dem Ritzschen Verfahren, lésen. 


4. Die Randbedingungen. Zu beachten ist noch, da8 man bei der Durchfiihrung des Ritzschen 
Verfahrens fiir die Funktionen x(s) und y(s) Naherungsfunktionen verwenden muB, die die Rand- 
bedingungen des Problems erfiillen. Man mufB also auch noch diese Randbedingungen kennen. Sie 
sind fiir die Variablen e; und e, gegeben und lauten bei den angenommenen Lagerfallen J, 2 und 3 
von Abb. 2 

G07, e, = 0 (26) 
fiir eingespannte Stabenden und 


e:—e,(€ + 7%) =90, e, = @ (€ + %) =0 (27) 
fur Stabenden, die von Biegemomenten frei sind. 


Durch die Beziehungen (9) und (13) findet man aus (26) und (27) leicht die entsprechenden 
Randbedingungen, die fiir die Funktionen x(s) und y(s) gelten miissen. 


5. Ein Sonderfall. Einer besonderen Betrachtung soll der Sonderfall 


unterzogen werden. Es ist dies die tordierte und gedriickte Welle ohne Vorverwindung, etwa mit 
kreisrundem oder quadratischem Querschnitt, wie sie von mir! behandelt worden ist. Fiir diesen 


Fall wird wegen @ = b | : 
rae coma (es) a ae (29) 


oder weil ce = T ist, 
ws) =o) = Vs) =—|(z,) +. (30) 


Zanichst einmal gibt das, wie man mit y = o = V aus (18) abliest, 


K, = V(s)x, K,=V(s)y- (31) 
Fihrt man nun die komplexen GréBSen 
rea+iy, o=¢+1e, (32) 


ein, so wird aus den Differentialgleichungen (15) 


7 = V(s)r. (33) 


In diesem Sonderfall ist also alles auf die Lésung dieser einen Differentialgleichung (33) zuriick- 
gefiihrt, die deshalb Basisgleichung heiBen soll. 


1 H. Leipholz, Il. 
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Kennt man die Lésungen r(s) von (33), so kann man die Lésungen eg und e, sofort angeben. 
Denn setzt man in 9 = eg + ve, die Beziehungen (14) ein, so 14Bt sich schreiben 


o(s) = r(s) e—™*, (34) 
und damit kann man auch die allgemeine Lésung des Problems und seine Kigenwerte ermitteln. 


In diesem Sonderfall ist also eine Integration der Differentialgleichungen méglich. 
Will man fiir diesen Sonderfall das Variationsproblem lésen, so erhalt man aus (20), unter Be- 


achtung von (30), zunachst 


U =— Vs) (@ +9) (35) 


und damit aus (23) 
La = [+2 + Ms) (+9). (36) © 


Das Variationsproblem lautet dann fiir die spezielle Annahme 


I 
q(s) = konst. und also Q(s) = f q(s) ds = q (/—s) 


l 
: ¢ LPN TE l 
[fe +54 [2 —(e{)— 4] @ +9] = Ear. (37) 
0 
Zur Vereinfachung der Rechnung seien die Abkiirzungen 
_ ql Bail Va ©. 
ot aad (z) aa 


eingefiihrt, womit 
I 


[ret hy-@—a 


(x? + x) ds = Extr. (38) 


folgt. 
Will man insbesondere den Lagerfall 1 von Abb. 2 zugrunde legen, so hat man wegen zweier 
eingespannter Stabenden die Randbedingungen 


= 0,5.6, = 0 - fir s S0ehusds =e (39) 
Das gibt wegen (9) und (13) fiir x und y die Randbedingungen 
a= OL 4 Sy Sia sc — Oita (40) 


Diese Randbedingungen werden von den Naherungsfunktionen 


<a ade 
x =4a,sins, y =a, sins (41) 


erfiillt. In diesen Naherungsfunktionen sind a, und a, unbestimmte Koeffizienten. 


Mit (41) wird aus (38) 


I l I 
2 
(aj + a3) F | ot ss —(K +4) | sint [sds +77 | s sig’ 7s ds] = Exte., (42) 
6 j 


0 


Die Ausrechnung der Integrale fiihrt auf die Bedingung 


; me KET eae 
Sage) (rays — Soe z| = Ese (43) 
Das Ritzsche Verfahren schreibt aber vor: 
as OF 
Di NG foe 
was die Eigenwertbedingung 
a (KA Ad 2al/a\? Eye 
21 US Weg ea ieee, a. 


zur Folge hat. 
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Hierfiir ein Zahlenbeispiel: In meiner zitierten Arbeit! waren folgende Zahlen gegeben: 
T = 104 kpm, Pao kp. q=60kp/m, a=10'kpn?. 


- Geht man mit diesen Zahlen in (44) ein, so erhalt man durch den zweiten Ausdruck eine Iterations- 
vorschrift fiir 1. Damit Konvergenz eintritt, mu8 man allerdings ein Verfahren zur Konvergenz- 
verbesserung anwenden, z. B. ein solches, wie es von R. Zurmiihi2 beschrieben worden ist. 

Beginnt man die Rechnung mit |, = 22 m, so erhalt man als I, = 22,6 m und dieser Wert 
stimmt mit dem von mir® gefundenen genauen Wert / = 22,5 m recht gut iiberein. Man hat also 
durch das Variationsproblem auf einfachste Weise in (44) eine sehr brauchbare Vorschrift zur Be- 
rechnung von Kigenwerten fiir den Sonderfall (28) und den Lagerfall 1 erhalten. 


6. Ein Versuch zur Transformation auf ,,zyklische“ Variable. Wendet man die Analogiebetrach- 
tung bereits auf das System (11) an, so hat man entsprechend mit 


Tle + & +2 u (G6, —tge,) + w+ ol, 
U = —- (y(s) & + o(s) 2), (45) 


1B ze as ‘ = Ss). (ez Ss 2 
La =T—U= 5) f+ +2 u (Ge, —ee,) —(0 +" EOS) ae —(p +2 £20) 0 | 


zu rechnen, wobei T = T, + T, + T, mit 


is g Lees : 
T.=5V(+4), T=u(%e,—ee), => (E+ %) (0) 
ist. 
In Polarkoordinaten sei 
€g = 1r(s) cos y(s) , e, = 1(s) sin y(s) . (47) 
Dann ergibt eine Umrechnung 
T=S[P+P(u+ OP], U=—+[y(s) 7? cos? + o(s) sin? g], 
(48) 


La = 2 [r? + r? (uw + q)? + y(s) r? cos? y + o(s) r? sin? gy] , 
was fiir (11) wieder die Darstellung 


d /éLa éLa _ 0, d Whe <s éLa ne (49) 
ds \ ér or ds\ dp oy 


erlaubt [diesmal aber mit La gemaf (48)] und damit zu dem Variationsproblem 


i [r2 + r2 (wu + ~)? + y(s) r? cos? p + o(s) r? sin? yp] ds = Extr., (50) 
6 


mit den Funktionen r(s) und ¢(s), fiihrt. ; 
Im Sonderfall a = b hat man o(s) = p(s) = V(s) und damit 


La = [# +2 (u + 9) + Vs) 1. (51) 


e , ; : ee : : ; 
Da La nun von der Variablen ¢ frei ist, so ist g eine »»zyklische“ Variable, und man kann eine 


Integration des Systems (49) durchfiihren. Wegen 0La/dp = 0 folgt 


d /aLa\ dy i= 0 DG i, OTS 6: (52) 
eee l= ly a, ee 
Setzt man insbesondere die Konstante C auf der rechten Seite gleich Null, so hat man, da r = 0 
trivial ware, paras (53) 


1 H. Leipholz, Il. 3M 
2 zn Teseinitbh Praktische Mathematik, S. 25—27, Berlin/Géttingen/Heidelberg 1953. 
3 H. Leipholz, II. i 
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Verwendet man (52) mit der Annahme C = 0, um T gemaf (48) zu schreiben, so gibt das 


r 
T=". (54) 


In bekannter Weise wird nun mit Hilfe der Routhschen Funktion 


pg ee Bee ed (55) 


die ,, Bewegungsgleichung“ fiir die ,nichtzyklische“ Variable r als 
d OR oR 6U (56) 


dsé; -ér .- or 


angesetzt. R ist durch (55) gegeben. Fiir U erhalt man mit o = yp = V aus (48) 
U =—5 Vs), (57) 
und so wird aus (56) 
= Visine. (58) 
was bereits als Basisgleichung (33) bekannt ist. 
Setzt man wieder 
(Li ee a i ens 
o=r er? , (59) 


und da nach (53) p bereits bekannt ist, wird wieder die Lésung fiir die eigentlichen Variablen e; 
und e, des Problems in komplexer Schreibweise durch 


@ =r(s) e~**s (60) 
gegeben, wofern man die Lésungen r(s) der Basisgleichung (58) kennt. Das ist genau wie schon 
durch (34) beschrieben. 


Damit ist gezeigt worden, daB fiir den Sonderfall a=bdie Auflésung der Dif- 
ferentialgleichungen (11) besonders einfach ist, weil dann eine der Variablen 
zyklisch wird. 


so gibt das mittels (47) 


Im allgemeinen Fall seien die ,,verallgemeinerten Koordinaten und Impulse“ 


1 = > da =e: 
aT | 8T — = RE OT). a (61) 
=p ST = Cte — UE! — => = 
Pi oan ee é n? Po = 8%, Be, e, + U eg 


eingefuhrt, und es sei in bekannter Weise die ,,Hamiltonsche Funktion“ 


H=T,—T,+U 62 
aufgestellt. Das gibt 2 : oe 
: 
H = — (pi + 4 gz)? + (py — u qi)? — (uv? + y(s)) gi — (uv? + o(s)) 3], (63) 
und hiermit erhalt man fiir das System (11) die aquivalente ,,kanonische Darstellung“ 
dy 2H dy _ Hd, dy 
ds ep,’ ds ap,” sae egS ds 3 ag 
bzw. bei Verwendung von (63) 
dq, ass dq, dp d 
ds Pit Uae» Oe io are a ae =U Pat ¥(s) Gi» Fe = — 4 Pi +9(s) qa (64) 


Die in der Kinetik bekannte Jacobi-Hamiltonsche Methode faB®t dieses System erster Ordnung 
(64) als das charakteristische System einer partiellen Differentialgleichung 


ow 2 ow ow \2 ow | | 
Ee +2u iq, 2 + ae —2u ee w(s) gi — o(s) g3 = konst. (65) 
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auf, deren vollstandiges Integral dann den Schliissel zur Lésung des ganzen Problems geben 
- wurde; denn die Riickfiihrung auf die partielle Differentialgleichung (65) stellt den systematischen 

Versuch dar, auch im allgemeinen Fall auf zyklische Variable zu transformieren. Dieser Versuch 
ist gelungen, sobald man ein vollstandiges Integral zu (65) angeben kann. 


7. Die elastomechanische Bedeutung des Variationsproblems. Die Gleichungen (11) und die 


Beziehungen (45) kénnen dazu dienen, den elastomechanischen Inhalt des Variationsproblems zu 
erklaren. Das soll jetzt geschehen. 


Unter Verwendung von La nach (45) hat man als Variationsproblem 
I 


I 
| ta ds =y|[e+4 +24@e,— ee) —(0 aed: += FES) ef] = Exe. 
0 0 


Das laBt sich mit T = c ¢, sowie mittels (4), (5) und (9) umschreiben in 
1 


1 
[ Peas ay [+08 +e +H) +r) — Tes, eee) + alot ate 
0 0 


+b (e+ 19)" 6! —[T (e + 1) + P+ O(s)] (8 + )} ds = Extr.. (66) 

Da jetzt T = 0 sein soll und das Variationsproblem von den homogenen Differentialgleichungen 
(11) abgeleitet worden ist, gelten fiir seine Anwendung folgende Voraussetzungen: 

1. eingespannte Stabenden (damit das axiale Torsionsmoment konservativ bleibt), 

2. querkraftfreie Stabenden (weil die Ausgangsgleichungen (11) homogen waren). 

Folglich ist insbesondere Lagerfall 1 von Abb. 2 zugrunde zu legen, fiir den diese Voraus- 
setzungen erfillt sind. 

Wendet man auf die gedriickte und tordierte Welle das Prinzip der virtuellen Arbeit an, so 
lautet dies, wenn A? die virtuelle Arbeit der inneren und A? die virtuelle Arbeit der 4uBeren Krafte 
ist, 

A? + AZj=0. (67) 

Die virtuelle Arbeit der auBeren Krafte setzt sich dabei aus einem Anteil Af von den Druckkraften 
und einem Anteil A vom Torsionsmoment herriihrend zusammen: 

Ai = Ap+ Az. (68) 

Wenn es nun méglich ist, diese virtuellen Arbeiten als negative erste Variationen eines Potentials 
darzustellen: 


A? + Ab + AY =—d6(l1;+1Ip + L1p) , 


so folgt aus (67) 
6 (i; + ITp + Iz) = 0 
und damit die Extremalforderung 


W=1,+p +p = Extr., (69) 


die als Prinzip vom stationaren Wert des elastischen Potentials bekannt ist. : 
Man wird vermuten, daB die Forderungen (66) und (69) in engem Zusammenhang stehen. LaSt 
1 


sich das zeigen, so ergibt sich als elastomechanischer Sinn von { La ds eine enge Verwandtschaft 
0 


zum elastischen Potential J/. : r . 
Die virtuelle Arbeit der inneren Krifte ist, bei Beschrankung auf die von den Momenten her- 


riihrenden Anteile, 
A? = — f (Mz dw, + M,, dw, + T de) ds . (70) 
6 


Jetzt sei beachtét, daB (3) gilt (mit T; = T = konst.). Setzt man dies in (70) ein, so erhalt man 


I 
A? = — f (amg dug + bw, dw, + T de) ds , 
0 


I 
A= —6 | (pack +p buh + Teds, 
0 
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woraus man fiir das Potential der inneren Krafte sofort 


t 
T= { (owt + bo} +2 T e) ds 
0 


abliest, was wegen der Bedeutung von w, und w, 
l 


T= | la ey + %e( +20)? + b(— ae + oe +)? +2 T e| ds (71) 


ergibt. 
Es ist die virtuelle Arbeit A}, von den Druckkraften herrihrend, 
I l I 
AY, a [p + fq as) (e¢ deg + e, de,) ds =f (P + Q(s)) (¢: deg + @, de,) ds , 


oder 
! 


seal 
A =6 | 5 (P +09) (8 +64) ds, (72) 
i) 
so daB man fiir das Potential 


Tp fe + Q(s)) (2 + e2) ds 


erhalt. 
Die virtuelle Arbeit A’, vom Torsionsmoment T herriihrend, ist 


l 
Ay = T f dr, ds. 
0 
Unter Verwendung von 
OT, = & OWs + &, 0, + € 
und. damit 


OT, = &- [de,, ab (e + To) deg] += en [= deg a (e sn To) de, A de 
gibt das: 


1 
Ap = ch [es de,, — @, dee + (€ + To) (€¢ de: + @, de,) + de] ds , 
was weiter auszurechnen ist. Sogleich hat man 
! 
: 1 : 
Ay 5 T{ [—2eSe, 4 2 o, See Dee) le Be, 1-6, de:) 2 del ds wc ae) 
: 0 


Ferner erhalt man durch partielle Integration 


1 ; Pi i, 
ae e: Oe, ds = — J 2252, ds Fiat ez Oe, ds — [ee de,], ; 


ol ee. (1) 
se e, Oe ds = H e, Oeg ds — f @, deg ds + [e, dee], : 
0 
Setzt man (74) in (73) ein, so kommt 
l 
; 1 ie ; ; 
Ay = ics T | [e, des + eg de,— e, deg — eg 6, — 2 (& + To) (eg Seg + e,, 0@,) — 2 de] ds 
: 
1 

aes T [—e, deg + ey den] ; (75) 


Aus (75) folgt aber 
l 


irae 1 
Ap =—5 | Tle ty —e 6, —(0 + T ) (e2 + e) — 2] ds +5 fa; de, —e, Sel. 
0 
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Da die Voraussetzung gemacht worden ist, daB das axiale Torsionsmoment konservatiy sein soll 
damit ein Potential existiert, so mu8 man Randbedingungen fordern, fiir die eg und e, an den 
Stabenden zu Null werden, namlich eingespannte Stabenden. Dann wird auch der torende Rand- 
term zu Null, und man erhalt 


l 
: 1 : 
Ap =—8 [5 Tyee — eb, —(@ + %) (e2 +e; ) —2e]ds, 
6 


so daB das Potential 
1 


1 : ‘ 
My = ae 1h | le, €¢ — eg &, — (€ + T) (e2 + e) —2 6] ds (76) 


0 


wird. GemaB (69) bekommt man jetzt mit (71), (72) und (76) fiir das Gesamtpotential 
l 


=> | (el, tele +P +? [— e, + e, (¢ + 1)? — [P + Q(s)] (e2 + #2) 
0 + T [e, @g — e e, —(e + t) (e2 + e)]} ds. 


Durch einfaches Umordnen ergibt sich hieraus 
I 


1 A : : 
Maz | (b+ ad —T (qs, — ee) + ale +r E+ d(e +r a 
@ [7 (e + 1) + P+ O(s)] (ef + ef) +20 (e +1) 7 6,—2b (6 +7) O a} ds. (77) 


Durch partielle Integration findet man 


; I I l 
2a(e at Fond Open ds — 3 (2 + %) J @n ee ds + a(e + %,) lien eelo— fe e, ash 
0 
(78) 
I I I 
—2b(e + to) J ey &¢ ds =— bd(e atta) eee toa. © (e + To) [en ely — f én oh 

0 


Setzt man (78) in (77) ein und beachtet noch, da wegen der Einspannung der Stabenden die 
Randterme [e, es]; wieder verschwinden, so hat man endlich 
l 
1 : é : : 
=z | ba +o, + [(@ + b) (e + t%) — T] (eg &, — eg &,) + a (e + tT)? ef +b (e + 1)? e, 
; —[T (e + t%) + P + Q(s)] (e2 + @; )} ds. (79) 
Der Vergleich von (79) mit (66) zeigt aber, daB 


I 
IT =a f Lads (80) 
0 


ist. Damit ist nachgewiesen, daB die Extremalbedingung (66) mit der bekannten Forderung nach 
dem stationaren Wert des eleastischen Potentials tibereinstimmt. 


Zusammenfassend kann man sagen: Das kinetische Analogen, bei dem die wirkenden Kriafte 
l 

aus einem Potential ableitbar sind, fihrt auf eine Extremalbedingung { Lads = Extr.; 
0 


diese ist einer elastomechanischen Extremalbedingung // = Extr. aquivalent, weil das hier be- 
trachtete elastomechanische Problem konservativ und deshalb fiir die virtuelle Arbeit der 
auBeren Krafte des elastomechanischen Problems ebenfalls ein Potential vor- 
handen ist. 

Als Letztes laBt sich noch zeigen, das fiir das hier behandelte Problem der Energiesatz fir die 
wirkliche Formanderungsarbeit A; und die von den virtuellen Arbeiten der duferen Krafte her- 
riihrenden Potentiale [7 und [/7 die Beziehung 


fi Oil, (81) 
liefert. 
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Da auBerdem bekannt! ist, da8 das Potential der virtuellen Arbeit der inneren Krafte gleich 
der negativen Forminderungsarbeit ist, so daB noch [/; = — A, gilt, folgt aus (81) 


Tee a, (82) 


Der Vergleich von (82) mit (69) und (80) zeigt aber, das der Wert des Extremums von //7 eben 
gerade Null sein muB. 

Wegen der Tatsache, da® bei den Extremalforderungen (25), (37), (50), (66) der im Integral 
stehende Ausdruck homogen quadratisch ist, muf ein eingliedriger Ritz-Ansatz, mit gleichem 
Beiwert fiir beide im Integral auftretende Funktionen, als Bedingungsgleichung fiir den Kigenwert 


(mit La gebildet fiir zulassige Vergleichsfunktionen) auf 
t= 
f Lads = 0 (83) 
0 


fiihren. Das bedeutet, da® hier das Raleighsche Verfahren giiltig ist. Da auBerdem, wegen der 
Anwendbarkeit des Energiesatzes, noch 


I 
f Lads = Extr. = 0 (84) 
0 


gilt, so folgt daraus die Méglichkeit einer Abschatzung zwischen dem wahren kritischen Wert des 
Problems und seiner Naherung. Denn man kann schlieBen, daB 


ae I 
f Lads = f Lads =0 
0 0 
sein mu, was eben die Abschatzung des kritischen Wertes erméglicht. 


8. Zusammenfassung. Die Anwendung der kinetischen Analogie auf das mit bestimmten 
Voraussetzungen formulierte Problem der tordierten und gedriickten Welle erméglicht es, die dem 
Prinzip vom stationadren Wert des elastischen Potentials entsprechende Extremalbedingung in der 
besonders handlichen Form 


l 

f (2? + y? + y(s) eZ + o(s) ef] ds = Extr. 

0 
darzustellen. Das dadurch erhaltene Variationsproblem erlaubt aber die Bestimmung von Eigen- 
werten fiir ganz beliebige Gesetze q(s) der Langskraftverteilung auch bei zwei verschiedenen Haupt- 
tragheitsmomenten von Stab oder Welle, eine Aufgabe, die auf dem Wege iiber die Differential- 
gleichungen zu erheblichen Schwierigkeiten fiihren wiirde. 

(Eingegangen am 23. Marz 1960.) 


Anschrift des Verfassers: Horst Leipholz, Stuttgart-Untertiirkheim, Fiechtnerstr. 51. 


1 A. Pfliger, Stabilitatsprobleme der Elastostatik, S. 37, Berlin/Gottingen/Heidelberg 1950. 
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Berechnung stationarer Temperaturverteilungen 


Von A. Kneschke 


1. Einleitung. Beider Anwendung des Ritzschen Verfahrens zur Lisung von Randwertproblemen 
der Elastizitatstheorie, die vom Prinzip des Minimums der potentiellen Energie eines elastischen 
KG6rpers ausgehen, sind Ansatzfunktionen erforderlich, die die auf der Oberflache oder Teilflachen 
davon vorgeschriebenen Werte der Verschiebungen annehmen. D. Riidiger! hat das Prinzip vom 
Minimum der potentiellen Energie durch Zusatzglieder dahingehend erweitert, daB auBer den 
Bedingungen fiir die Oberflichenkrafte auch die Randbedingungen fiir die Verschiebungen als 
natirliche Randbedingungen des Variationsprinzips auftreten. Damit entfallt hinsichtlich der 
Randbedingungen die erwahnte Beschrankung in der Wahl der Ansatzfunktionen. Fiir das er- 
weiterte Minimalprinzip von der potentiellen Energie sind die Ritzschen Ansatzfunktionen somit 
an keinerlei Randbedingungen gebunden. 


Das von E. Trefftz* zur genaherten Lésung von Randwertaufgaben als Gegenstiick zum Ritz- 
schen Verfahren angegebene Verfahren der partikularen Lésungen hat durch die Erweiterung des 
Prinzips vom Minimum der potentiellen Energie eine beachtenswerte Anwendungsmiglichkeit er- 
fahren. Die als partikulare Lésungen der elastischen Grundgleichungen herangezogenen Ansatz- 
funktionen, die die Giltigkeit der Gleichgewichts- und Kompatibilitatsbedingungen im Kérper von 
vornherein gewahrleisten, brauchen keine Randbedingungen zu erfiillen, wenn vom erweiterten 
Minimalprinzip ausgegangen wird. Daf im Gegensatz zum Ritzschen Verfahren hierbei keine Ge- 
bietsintegrale, sondern nur Randintegrale auftreten, ist obendrein ein nicht zu unterschatzender 
Vorteil des Trefftzschen Verfahrens. Natiirlich wird man Randbedingungen, die sich ohne weiteres 
durch die Ansatzfunktionen erfiillen lassen, zur Vereinfachung der Rechnung bereits im Ansatz 
beriicksichtigen. 


Im folgenden soll ein fir stationére Warmeleitvorgange giiltiges Minimalprinzip angegeben und 
zur Berechnung der stationdaren Temperaturverteilungen das Trefftzsche Verfahren herangezogen 
werden. An zwei Beispielen wird der Vorzug dieses Integrationsverfahrens gegeniiber dem von 
Ritz besonders deutlich hervortreten. 


2. Variationsproblem und Trefftzsche Gleichungen. Im Innern eines festen Kérpers K mit der 
konstanten Warmeleitzahl J besteht eine vom Ort abhangige Warmequellung mit der Ergiebigkeit W 
je Volumen- und Zeiteinheit. An den Teilflachen F,; ({ = 1,2,..., n) der Oberfliche F von K 
sollen die Temperaturen T;, an den Teilflichen Fy; (j =1,2,...,m) die Warmestréme 


oT 
on Fo; 


—A =—Anegrad T| =} 


2j 


vorgegeben sein, wabrend an den restlichen Teilflichen Fy), (k = 1, 2, ..., /) ein natiirlicher Warme- 
iibergang mit den Warmeiibergangszahlen «, bei den AuSfentemperaturen T7,;, stattfindet. Fiir die 
im Kérper entstehende stationdére Temperaturverteilung T besteht dann die Randwertaufgabe 


AAT =A div grad T=— W, 
Peer pt: Gee een) 


oT : (1) 
ee Ua = ® (ji Vim tf) 5 
oT 
lA oe +o Cat (k =1,2,...,]). 


1 PD. Riidiger, Ing.-Arch. 27 (1960) S. 421. 
2 EB. Treffiz, Gee eet, Kone. f techn, Moch., Zirich 1926, 5.131. 


118 A. Kneschke: Berechnung stationarer Temperaturverteilungen Ingenieur-Archiv 


Sie ist, wie sogleich gezeigt wird, der Variationsaufgabe 


HT} = | 


K 


1 
Sy Lagi os? y [ o(¢ 7 = Extr. (2 
a [Ata ig lo + fo T do - [ (3 eee T) do Ear) 
fie ee 27 ~~ F345 


4 (grad T)?— W r| dr — 


gleichwertig. Die Extremalbedingung 


OJ = f(Agrad T grad dT — WOT) dr — 
K 


Se a ce tT) do— > [ 25n Te 


Fy; coh as 
+> [ seta +> | ax? — Tax) 9T do = 0 
J Fj k F 3}, 


geht hei Beriicksichtigung der Greenschen Forme] 


[ era T grad 0T dr =—[ ATOT dr + [ Fe OT do 
K K F 


und von 


in 


sJ=—<fUAT + W)dTde— 
K 


=> Ct, [ (5p + 95) 97 do sf 
: Fy; 


ee 


+> [Aa toe (7— T,)|6T do =0 
FF 


uber. Die Warmeleitgleichung 1 AT + W =0 ist demnach Eulersche Differentialgleichung und 
die Oberflachenbedingungen sind natiirliche Randbedingungen des Variationsproblems (2). Das 
Treffizsche Verfahren der partikularen Lésungen ist auf (2) anwendbar. 


An die approximierenden Funktionen v, v, des Ansatzes 


me P 
T = v% + x Vy 
werden die Bedingungen 


AAvy+W=0, dv, =0 (=, 255 Pp) 


gestellt. Da dann —AAT = W gilt, besteht nach dem Gaufschen Integralsatz die Beziehung 
, a oT oT 7 
ree a [aea> [ae=-t [4Te =| War. 
an TFS; * F3, oS oS 


Hierin kommt der Energiesatz zum Ausdruck, da der Uberschuf der ausstrémenden Warmemenge 
uber die einstrémende gleich der im Kérper in der Zeiteinheit freiwerdenden Warme ist. 
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Mit dem Ansatz fiir T ergibt sich aus (2) die Extremalbedingung 


a a grad T grad v,— Wv,) dr — 


= [a@— mo — >» [Agr edo + 
Fi; Fy 


Wir wenden auf den ersten Teil des Raumintegrals wieder die Greensche Formel 
= oi eT 
| era T grad v, de=— | AT», de + [Be 
K K ra 


an und erhalten bei Beriicksichtigung von 1 AT + W = 0 mit w=1,2,...,p die Trefftzschen 
eee 


l aL ve 
> far 11) 28 dp — [Qt +a)», o> ae (3a) 
i= 2 Fi; j=1 Fo; k=1 F3), 
oder 


=> Piss wae] eae o+ > | fas 


k=1 F3), 


Ds 7 fans st do — . See && eydo— >) J 03 a ont) 84 do 
vy=1 


Ste oa (te wt) do, (3b) 


Fir die allgemeinen @rthosonalen Koordinaten u,, u,, vu; mit den metrischen Fundamentalgréfen 
bzw. Laméschen Koeffizienten 


ee a ee 
Boa = Le = (ena Aa 1 (ine EEE 


gilt bekanntlich 


und damit ist in (3a) und (3b) allgemein 


P : (ns) 2 
= cos (n, Ug 
an pile 2 fe dUg 


o=1 


zu setzen. Hierbei sind die cos (n, u,) die Richtungskosinus der nach auBen weisenden Flachen- 
normalen und die ¢, die Basisvektoren des allgemeinen orthogonalen Koordinatensystems. 


Fir die erste Randwertaufgabe mit der Randbedingung 
fees IS 
lauten die Treffizschen Gleichungen 


P 

0 
> leg 2 aa ie ae (4 =1,2,...,p), (4a) 
v=1 F Ee 
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fiir die zweite Randwertaufgabe mit 


ov Ovy 
| (v, TAP Vn Av,) dt = | (°»n es, si do 
F 


wegen Av, = Av, = 0 


folgt, 
Cc le do ria A gale do (i = 2 Pp) 4b 
> on on } H ans oie i ( ) 


a 
Fiir die dritte Randwertaufgabe mit 


peer 
F 


a On 


bestehen die Gleichungen 


Pp 
Die | (Age + 2%) do = — | [ASP +o (09 — .)| v,, do (i= 1, 2h. pe 
TY F 


3. Beispiele. a) An der Grundflache eines Kreiszylinders vom Halbmesser R und der Héhe a 
herrscht die konstante Temperatur T,, an den iibrigen Flachen die Temperatur Null. Die stationare 
Temperaturverteilung T(r, z) ist Lésung der Randwertaufgabe 


o*T t oT eT 
AT= oe + > or 4a 2 9 (5) 
T\,—0 = faces YW = 0 ) én ee = 0. 


Fiir die Koordinatenfunktionen 1,,, v,, des Ansatzes 


bol P i 
T= 2 (4, Vy se Coy Vpy) 


mit 
Bae =t de (3, R) cosh 6, ‘e - 
- a AC) =0 
Voy = Jol>s 3) sinh o, R 9 
gilt 


Av,, = Ary, =0, %,(R, 2) = %,(R, 2) = 0. 


Die aus (3a) hierfiir abgeleiteten Trefftzschen Gleichungen lauten 
R 


R a 
oe 7 Op pw mm OV, 
[re | eae ewe! Pe | = 0(0 =1,2; 7 =1,...,p), 
0 0 


woraus wegen 


Ov. 
i) 1» 
2=0 : 0z 


z=0 
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nach Einsetzen von T 


R 


Pp es R 
v1 
> few [ ee | r dr + cay | 
a! Gee 
0 0 
R 


R 
. ov ov 
au 2M 
> eo 01,5, | ..rdr—Y |v,,—*| rdr c 
v=1 ; {| "oa Ie {| 1” 8s The ao 
0 0 0 


folgt. Bei Beachtung der Integralrelationen 


a 


R 0 fir ee a ae 
ae (5.7) J (3 R)® =) R 
i DLR Le ae = J%(5,) fir w=», 


R 


r R 
: r Jo(O RE) = 5, (0,) 
0 
gehen diese Gleichungen in 


c,, cosh 6, z + ¢», sinh 6, z == 0% 
Cis sinh}, > +c, cosh bz SS eee! 
dy Jy(5y) sinh 6, x 
uber. Daraus ergibt sich 
eT, saree aa a 
pee hey 1 eggs) Oo 
und die genaiherte Temperaturverteilung 
Soe liareaey 2 ae 
T(r, z) =2T, 2 BIA) [cosh Oy Fast coth 6, R sinhd, x) ; (6) 


Fiir p oo geht diese in die strenge Lésung! des stationdren Warmevorganges tiber, die man 
durch unmittelbare Integration des Randwertproblems (5) erhalten kann. 


b) In einem unendlich langen Stab mit einem quadratischen Querschnitt von der Seitenlange a 
findet eine Warmequellung mit der konstanten Ergiebigkeit W je Volumen- und Zeiteinheit statt. 
Es soll die stationire Temperaturverteilung T(x, y) ermittelt werden, wenn an der Staboberflache 
die Temperatur Null herrscht. 


Dem Warmeleitvorgang liegt die Randwertaufgabe 
OL OT W 
oa yt Ah (7) 


Thee =0 p) Dees =0 p) De af2 = 0 ? Tap =0 


zugrunde. Fiir den Ansatz 


P P 

— W : 

T =v + > Cy Vy = 59 (ax — 2%) + » c, sin 6, — cosh 6, =, 
ve v=1 


in dem 
0, = V1 


gesetzt werden soll, lauten die aus (3b) ableitbaren Treffizschen Gleichungen 
P 
ov 1) 
Dei aa | 25! « (ia Dp) > 
ae a (8 (C) 


1 4, Kneschke, Differentialgleichungen und Randwertprobleme, Bd. 2, S. 300, Berlin 1960. 
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in denen die Integrale itber den Querschnittrand (C) zu erstrecken sind. Man erhalt 


a 0. tir fp =H 


Opn Oy 
0G ds = 2 | [oF dx = 5) 6 

on CY |y als inh <4 OB. ss SS 
® ; y=a/2 6, sinh 5 cosh 5 fir p= pi 3 


a 


(mov Ie Gin 
tojn d= 2 | |e a os 3 sinh Of 
(6) 6 bs To te fir pp = 1,3, 02. 


und damit bei ungeradem p 
40 fur 2,4,...,p—l, 
4a? W 1 

A 


Oy 
6: cosh > 
Die stationire Temperaturverteilung wird also durch 


1 : i ene a: 
W a ax —x? [5 | ain Oh oad x cosh Gy })a, . 
a a 1 (2 y — 1)? 2? cosh (2¥7—)a 


v 


genihert dargestellt, wobei |? a die groBte ganze Zahl bedeutet, die kleiner oder gleich © = : 


ist. Fiir p oo geht die Naherungslésung in die strenge Lésung ier. 


(Eingegangen am 1. April 1960.) 


Anschrift des Verfassers: Professor Dr.-Ing. A. Kneschke, Freiberg (Sachsen), Richard-Wagner-Str. 13. 
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Der Strémungseinflu8 auf den Wellenwiderstand von Schiffen* 
Von F. Kolberg 


1. Einleitung. Schiffsmodellversuche zur Klarung des EinfluBes einer Stroémung auf den 
Widerstand von Schiffent? haben ergeben, dafs inshesondere der Wellenwiderstand durch die 
Strémung stark beeinfluBt wird. Es zeigte sich dabei, daB bei gleicher Maximalgeschwindigkeit 
in der Strommitte am Ort des Schiffes der Wellenwiderstand bei Fahrt gegen den Strom kleiner 
ist als bei Fahrt auf stehendem Wasser, wogegen bei Fahrt mit dem Strom der Wellenwiderstand 
gréBer ist als bei Fahrt auf stehendem Wasser. 

Die vorliegende theoretische Untersuchung soll zeigen, daB diese auf Grund von Versuchen 
erhaltenen Ergebnisse im Rahmen einer linearen Theorie in befriedigender Weise erklart werden 
kénnen. Wir werden zunachst den Fall eines exponentiellen Geschwindigkeitsprofils behandeln, 
da hierbei das auftretende Randwertproblem exakt lésbar ist. Spdter untersuchen wir dann, wie 


weit die so erhaltene Lésung eine Approximation darstellt fiir den Fall anderer Geschwindigkeits- 
profile. 


2. Differentialgleichung des Problems und Lésungsmethode. Das Schiff bewege sich geradlinig 
gleichférmig mit der Geschwindigkeit U, an der Wasseroberflache eines seitlich unbegrenzten 
Stromes der Tiefe h (Abb. 1). Wir beziehen uns auf ein schiffsfestes kartesisches Koordinaten- 
system; die x-Achse weise entgegen der Fahrtrichtung, die y-Achse liege in der Wasseroberflache 
und sei senkrecht zur x-Achse und schlieBlich weise die positive z-Achse senkrecht nach oben. 


Abb. 1. Koordinatensystem und Bezeichnungen. 


Sind dann 9 die Dichte, g die Erdbeschleunigung, p’ = p —o gz der Druck, U(z) + u,v, w 
die Geschwindigkeitskomponenten in x-, y-, z-Richtung, so lauten die linearisierten Bewegungs- 
gleichungen fiir unsere stationére Bewegung 


Gee ie el Op: 1 
ee OE iti 0 ex? (1) 
dv ta OP 2 
Boe Big He (2) 
ow RCP a (3) 
“6x Mig ne On OR 


* Veréffentlichung der Arbeitsgemeinschaft: Institut fir Mathematik und GroBrechenanlagen, TH Aachen, 
Versuchsanstalt f. Binnenschiffbau, Duisburg, Inst. a.d. TH Aachen. — ; 
1 W. Graff, Der Strémungseinflu8 auf den Form- und SES oR on aaa von Binnenschiffen. For- 
bericht Nr. 333 des Wirtschafts- und Verkehrsministeriums ; ; 
ey, rai Wieeeauchube der in stehendem und strémendem Wasser festgestellten Anderungen des Schiffs- 
widerstandes durch Druckmessungen. Forschungsbericht Nr. 618 des Wirtschafts- und Verkehrsministeriums 
NRW. 
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Dabei ist U(z) das Strémungsprofil, in welches die Schiffsgeschwindigkeit U, mit einbezogen ist. 
Hinzu kommt die Kontinuitatsbedingung, welche fir die Fundamentallésung, eine Quelle der 


Starke Q im Punkte (é, 0,2), lautet! 
a 4 8 4 OO 8(x — 8) (y) Be —0) (4) 


Hierin ist 0(t) die Diracsche Delta-Funktion. Die rechte Seite der Gleichung (4) ist also Null in 
allen Punkten mit Ausnahme des Punktes (é, 0,¢), aber das Volumenintegral der rechten Seite 
von (4), erstreckt iiber irgendeinen den Punkt (é, 0,¢) enthaltenden Raumbereich, ist Q. 


Differentiiert man die Gleichungen (1), (2), (3) nach x baw. y baw. 2 und addiert dann diese 
drei Gleichungen‘, so erhalt man unter Beriicksichtigung der Kontinuitatsbedingung (4) 
dU éw 


see +E + 20S FE =— 00 Ul) Oe —8) ly) He —2)- 


Wegen (3) kénnen wir hierfiir schreiben 


dU 
p-., Op ep di: Opa 35 Aye. 23 Dr 
fo 1 Oey Oe 2 & P_. 99 UGE)8(@—8) dy) e—0), (5) 


womit die Differentialgleichung fur den Druck p(x, y, 2) abgeleitet ist. 


Durch Integration der Gleichungen (1), (2), (3) ergeben sich die Geschwindigkeitskompo- 


nenten? zu 
x 


dU 
1 dz 
= —_——* ae 6 
3 o U(z) P aah Na es (6) 
ae 1 ; ap(A, y> 2) 
2 UG | en 2 : 
7 1 = aptA, y> 2) | 
w =e iG) | = dd. (8) 


3. Die Randbedingungen an der freien Oberflache und die Gleichung fiir die freie Oberflache. 
Betrachtet werde ein auf den Bereich B der Wasseroberflache wirkendes Drucksystem p(x, y). 
AuBerhalb des Bereiches B ist der Druck an der Wasseroberflaiche z =¢(x, y) stets gleich dem 
konstanten Atmospharendruck py, den wir ohne Einschrankung der Allgemeinheit als py) = 0 an- 
nehmen kénnen (Abb. 2). Bezeichnen wir nach Stokes mit D/Dzt die substantielle Ableitung, so 
gelten an der Wasseroberflache z = C(x, y) die Randbedingungen 


Dp’ _ ap’ op’ or 
eee ge © ee Us (9) 
: fiir z == C(x, ¥) 
eS ac 
Dem ee eT Gy” =w, (10) 
In (9) ist 
P ~Pp—ogz.- (11) 


Werden C(x, y), 0¢/dx , 0C/dy ebenfalls als von erster Ordnung kleine GréBen vorausgesetzt, so 
erhalt man durch Linearisierung der Gleichungen (9) und (10) wegen (11) 


op a ap 
og aot re aes ae 
at (auBerhalb B) . 


1 N. J. Lighthill, J. fluid mech. 3 (1957/58) S. 113. 
® Th. v. Kaérmén u. H. S. Tsien, Quarterly appl. math. 3 (1945) S. 1. 
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Differentiieren wir (12) nach x und beriicksichtigen wir noch (3), so erhalten wir als Randbedingung 
an der freien Oberflache 


ep a : 
dat | UG) Oz xk? aes us 
Eliminiert man w aus (12) und (13), so erhalt man 
ac lee Ls 
ax OB ox, LP t Pli=o » 


und durch Integration nach « zwischen den Grenzen —oo und x ergibt sich hieraus die Gleichung 
der freien Wasseroberflache zu 


Clsy) == Lp + Plno- (15) 


Hierbei wurde beachtet, da8 C(x, y) >0 und p(x, y, 2) +0 fiir x ~—0oo gelten mul. 


Abb. 2. Koordinatensystem und Bezeichnungen. 


Mit den Gleichungen (14) und (15) haben wir die Randbedingung an der freien Oberflache und 
die Gleichung der freien Oberflache aufgestellt fiir den Fall, da®B auf den Bereich B der freien 
Oberflache ein Drucksystem wirkt. Wird jedoch ein Michellsches Schiff betrachtet, bei welchem 
die Schiffsoberflache erzeugt wird durch eine in der Lingsschiffsebene angebrachte Quell-Senken- 
belegung, so hat man in (14) und (15) p(x, y) = 0 zu setzen und erhalt als Randbedingung an 
der freien Oberflache 


ep pa ae e3 Nek 
aa = U%@ oe 0:, fin-7 =0" (16) 


Die Gleichung der freien Oberflache ergibt sich hier aus (15) zu 
Mh 
C(x, ¥) =—— [p].-o - (17) 


4, Bestimmung der Elementarliésung bei unbegrenzter Strémung. Betrachtet werde eine Quelle 
der Stirke Q. Die Quelle liege im Punkt x =&, y=0,z=€. Das Geschwindigkeitsprofil der 
Anstrémung sei U = U(z). Nach (5) haben wir dann fiir den Druck die Differentialgleichung 


dU 


ap 0*pi, (02p (ISS | er oa ite a 
Peay ee be 0 Q U(z) 6'(x —é) d(y) d(z—¢) . (18) 
Wir nehmen an, daB der Druck p(x, y,z) eine Fouriertransformierte p(s, t, z) besitzt: 


Bls,t,2) = ff plxy,2) exp [i(sx +ty)] de dy. (19) 


—00 —0o 
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Wenden wir auf die Differentialgleichung (18) die Fouriertransformation an, so erhalten wir fiir 
die Fouriertransformierte p(s, t, z) die gewohnliche Differentialgleichung 


dU 
BP EO +P 
=—9 QU) d@—2) [ { exp[i(s +ty)] 8% —€) Hy) de dy. (20) 


Beachtet man die Relationen 


fe) le—9) =F) 9 —9) 
SFr) 00) ty = £00) 


F fue) 0 — 8) de = — J fuls) (x —8) de =f), 
so erhalt man aus (20) 
cereal ee 
oF _ 2 & BP _ (2+ 2) p =i 9 QU(C) de —C) exp [is]. (21) 


Beschranken wir uns jetzt auf den Fall eines exponentiellen Geschwindigkeitsprofils 
U(z) = V exp [kz], (22) 
so geht (21) tiber in die gewéhnliche Differentialgleichung mit konstanten Koeffizienten 
oP 2k (2 + 2) p ig Q U(l) sz —L) exp [is]. (23) 
Die allgemeine Lésung der zu (23) gehérigen homogenen Differentialgleichung ist 
Pr(s, t 2) = exp [k(z —)] - {4 exp [(@ —2) fi? + * +] + Bexp[—@—O fP+% +e}. (24) 
Gesucht wird hier eine Lésung von (23) mit den folgenden Eigenschaften: 


1) p(s, t, z) soll fir z—€ >0 der Differentialgleichung (23) geniigen und fiir z > + co den 
Grenzwert Null haben. Nach (24) folgt daraus 


p(s, t, z) = Bexp [k(z —C)] exp [— (ze —O VR +3?+ 8] , fur 2—CS 0. (25) 


2) Fir z—¢ < 0 soll p(s, t, z) der Differentialgleichung (23) geniigen und fiir z ->—oo den 
Grenzwert Null haben. Daraus folgt 


P(s; t, 2) = A exp [k(z —Q)] exp [(z —0) Vl? + 2? + 2] , fir 2=-C <= 0. (26) 
3) Fiir z =¢ sei p(s, t, z) stetig. Nach (25) und (26) muB also 
A=B (27) 


sein. 
4) Die erste Ableitung von p(s, t, z) habe fiir z >¢ eine Unstetigkeit. Und zwar gelte 
lim p’'(s, t, 2) — lim p(s, t,z) =io Q U(¢) sexp [is¢]. 
2 ae 0 Gees 0 


Unter Beriicksichtigung von (25)—(27) ergibt sich daraus 


A=B=—5 Q U0) Feige lied 


und somit 
P(s, t, 2) a —£ Q U(¢) 


a exp [k (z —C)] exp [—|z —¢| /? + s? + 2] exp [is ¢] 
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oder wegen (22) 


n — Q is = 2) Fk Lo Bo 
P(s, t, 2) = —5 Q Ul) (| ae [—|z—¢| /i? + 8? + 2] exp [is é]. (28) 


Das Inversionstheorem der Fouriertransformation ergibt dann 


Ml is : 
PHI 2) =— $0 UG) ga [ a [—i {s(x —£) + ty} 


—|z—C| ke + 8? +] ds dt. | (29) 


Aus einer Tafel der Integraltransformationen! entnimmt man 


Se ers ts ei) e+e s+ 2] Pr yes exp [—k + y+ 2°] 
je+ete 2% jetyfre ~ 


—0o —00 


Hieraus folgt durch Differentiation nach x 


V+ s 


= | a gage: al 


—00 


=== 1 exp[—k x2 + y? See lle 
25 Coals ene a 
Hiermit erhalten wir aus (29) 
a x k V(x ty? + (2 —C)? 
OS ae (> oe Vee ae ai Rima ) ne 


_ und die zugehérigen Geschwindigkeitskomponenten u, v, w ergeben sich aus (6) bis (8) in Uberein- 
stimmung mit Lighthill? zu 


_ @(_ Q expl—kV¥@—F¥ + ¥ +E 
oA) 5. ee V@—F2 + + EO) : 
Oke , y) fewl-Ki@ er tt e—OF 
vrs ie mn ) | Va—s+y en 2 see 
ole ys) <2. {0 op -H@ =P TFT EH (32) 
ee y\ ae eRe t ey!’ 
wie Q. exp [—k Je — 2 $y? + E—O 3 
w(x, Ys 2) =(a + ')( An Ve= 2) + 97+ op ) ss 


Damit ist die Grundlésung einer Quelle in einem seitlich und der Tiefe nach nicht begrenzten 
Strom mit exponentiellem Geschwindigkeitsprofil abgeleitet. 


5. Die Elementarlésung bei flachem Wasser und freier Wasseroberfliche. Wir nehmen nun- 


mehr an, daB das Strémungsfeld begrenzt sei durch den Tankboden z = —h und durch die freie 
Oberflache bei z = 0. Gesucht ist also eine Lésung der Differentialgleichung 
a , 
oe ee Be 2 ee = 00 UL) 8 8) A(y) Bly) Hl) (34) 
die der Randbedingung am Tankboden w(x, y, —h) =0, also nach (10) 
&p = 0 (35) 
Oz |, 4 


1 A. Erdélyi, Tables of integral transforms, New York 1954. 
2 Siehe FuBnote 1 von S. 124. 


g* 
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und der pecan an der freien Oberflache 


F) bg / 
- Po i = 0, firz=0 [U(0)=V] (36) 


genugt. 
Mit (30) haben wir nun eine Lésung der inhomogenen Gleichung (34) bestimmt. Beachten 
wir noch die Integraldarstellung (29), so kénnen wir mit 


s=Kcos@0, t= Ksin O 
die Lésung (30) auch in der Form 


p(x ¥,2%) =— wee Je os 0 dO i ~ exp [—i K {(x —&) cos 9+ y sin 0} —|z —C| K*] dK 
: (37) 
schreiben. Dabei wurde zur Abkiirzung 


‘(RR (38) 


gesetzt. Zur Anpassung an die Randbedingung (35) addieren wir zu (37) den Druck 


Pil, ¥> 2) = ee [wou [x ; exp [—i K {(x — &) cos O + y sin O}] 
x F(O, K) exp [— K* z] dK, (39) 
der eine Lésung der zu (34) gehérigen homogenen Differentialgleichung ist. Die Randbedingung 


am Tankboden z = —h 
a(p + Ps) 


dz z=—h : 
liefert 
K* +k 
F(0, K) =— pe, op [-K* 24+). 
Damit wird 
U 
Pi(% ¥> 2) = — oa) os 0 dO [ Ze arse em exp [—i K {(x —€) cos 0 + y sin O}] 
x exp [— K* ees dK. (40) 


Zur Erfiillung der Oberflachenbedingung (36) addieren wir zu p(x, y,z) + p,(x, y,z) einen 
weiteren Druck p,(x, y, z), der bereits der Randbedingung 


Pa 


Oz 9 


z=—h 


genigt. 


Fiir p(x, y, z) machen wir den Ansatz 


e L UG r r 
P(X ¥.2) = ANS Pai 0 dO | ae exp [—1 K {(x —&) cos 9 + y sin OV] 
j 7 


meet 4 


| x {6,(0, K) cosh [K*(z + h)] + 6,(0, K) sinh [K*(z + h)]} dK . 
Die Randbedingung am Tankboden 


£Pa 


OB Se 


42=—h 
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liefert dann 
h 
C,(0, K) = —=; G,(0, K). 


Setzen wir G,(0, K) = G(0, K), so ist also 


> U hi4 co 2 
PAX, ¥s2) = Se cos 0 dO | = exp [—i K {(x —&) cos © + ysin 0}] 
1 0 
G(O, K) [eosh [K* (« + h)] —e sinh [K* (2 + Di dK (41) 
eine der Randbedingung 
ép. 
“a z=—h =f 


geniigende Lésung der zu (34) gehérigen homogenen Differentialgleichung. Die Summe p(x, y, z) 
+ py(x, y, 2) + p(x, y, 2) muB jetzt noch der Randbedingung (36) an der freien Oberflaiche 
gentigen. Mit (37), (40) und (41) ergibt sich daraus G(O, K) zu 


__ 2exp[— K* h] K* {cosh [K* (¢ + h)] — (k/K*) sinh [K* (¢ + h)} 
ee) (K* —k) K cosh [K* h] ; 


* 
K4 Ke (5 —x) 0 


Kok 
* “K* —(K) + k cos? ©) sek? © teh [K* hl] * 


Hierbei wurde K, = g/V? gesetzt. 


Damit erhalten wir 


K* 
K+ Ko( —— ——] sek? O 
(x, y,2) = wore weil (x 2 
Pa% Y> Ao K* — (K, + k cos? @) sek? © tgh [K* h] 


yz {cosh [K*(e++h)] —(h/K*) sinh [K*(e-+h)}} {cosh [K*(C + h)] —(h/K*) sinh [K*C + h)} KK* 
(K* — k) cosh [K* h] 


Xx exp [— K* h] exp [—i K {(x — &) cos 0 + y sin O}] dK, (42) 


oder, wenn dieser Ausdruck noch etwas umgeformt wird, 


a 


Pt, ¥.2) = oe [es 0 dO 


ee 


[o.e) 


ak feos ney + h)] = sinh [K*(z + ny feos [K*(¢ + h)] Ke sinh [K*(¢ + h)] 
x 


cosh [K* h] 
0 


x exp [—i K {(x —) cos @ + y sin 0} — K* h] 


K* —(K, + k cos* @) sek? © 43 
x K* — (K, + k cos? @) sek? © tgh [K* h] us o 


Dieses Integral ist im Sinne des Cauchyschen Hauptwertes zu verstehen. 
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Mit (30), (40) und (43) haben wir damit fiir den Gesamtdruck Pg einer Quelle in exponentieller 
Grundstrémung, wenn jeweils die Realteile der Ausdriicke (40) und (43) genommen werden, . 
> fre Q exp[—kV@—SF +¥+E—OF 
Pg(% ¥> 2) = —O Ares res W@—+y + @—o) 


co 


n/2 
0 Q Ulz) Ree 
a cos 0 d@ K* Kk 
0 6 


< exp [— K*(z + 2 h + 2)] cos [y K sin Q] sin [(x — &) K cos 0] dK 
n/2 fee) 


aoe R feos [K*(2 + h)] — = sinh [K*(e + i} 
wea | CO OOO) | rela ae 
0 0 


ok K* {K* -- (Ky + k cos? @) sek? 0} 
x feosh [K* (¢ + h)] — xx sinh [K* (¢ + h)]} ae (K, 2k a @) sek? © tgh [K* h] 


Xx exp [— K* h] cos [y K sin 0] sin [(x — ¢) K cos O] dK. (44) 


Dieser Druck p,(x, y, z) erfiillt alle bisher geforderten Randbedingungen. Nun ist aber klar, 
daB weit vor der Quelle (« -~—.o) keine Druckstérungen bzw. Wellen vorhanden sein diirfen, 
d. h. es muB 

lim p, = 0 


gelten. Um nun den Grenzwert lim p,(x, y,z) zu untersuchen, verfahren wir wie folgt. Der 
x => co 


Grenzwert lim p,(x, y,z) des ersten Summanden in (44) ist selbstverstandlich Null. Von dem 
Integral *>+° 


So pe eip [— K* (2 + 2h-+)]- cos [y K sin 0] sin [(x —€) K cos 0] dK 


1aBt sich durch Aufspaltung des Integrationsintervalles (0,00) in die Teilintervalle (x, A) und 
(A, co) unter Benutzung der Dirichletschen Formel nachweisen, daB es fiir x > + 00 gleichmafig 
in 0 (0 < O <a/2) gegen Null konvergiert. Damit ist auch der Grenzwert des zweiten Sum- 
manden von (44) fiir x >-+ oo gleich Null. 


Zu untersuchen bleibt also lediglich der dritte Summand von (44). Der Integrand dieses Sum- 
manden hat Polstellen in den Wurzeln von 


K* —(K, + k cos? @) sek? © teh [K* h] = 0. (45) 


Bei festem 0 (0 S O <z/2) ist von diesen Wurzeln genau eine reell und gré8er als Null, 
weshalb bei der Integration beziiglich K der Cauchysche Hauptwert zu nehmen ist. Zur Unter- 
suchung des Grenzwertes fiir « +-++ oo schreiben wir den dritten Summanden von (44) in der 
Form 


: m/2 
ala, ys2) = Re 20-00) i 68 6 
0 
| feos [K*(e +] — eq sinh [Ke + ny feos [K*C + hy] — 5 sinh [K*E + ny 
2s be cosh [K* h] 


0 
K* {K* + (Ky + k cos? @) sek? O} 
K* — (Ky +k cos? @) sek? © tgh [K* h] 


xX {exp ee K {(x —) cos O + y sin O}] 
+ exp [—i K {(x —&) cos 9 — y sin O}]\ dK . (46) 


exp [— K* h] + 
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Bevor wir dazu tibergehen, in dieser Integraldarstellung den Integrationsweg der Integration 
_ beziiglich K umzuformen, wollen wir die Lésungen der Gleichung (45) diskutieren. Wir schreiben 
(45) in der Form 

K*h 


Ky, h ( + — cos? 6) sek? © 
0 


= tgh [K* h]. (47) 


' Fir die Existenz reeller dieser Gleichung geniigender K* h > 0 ist die Bedingung 
I 
Ky (1 +e cost 6) sek? @ > 1 
0 


notwendig und hinreichend. Da uns hauptsichlich die Singularitat auf dem Integrationsweg 
K = 0 interessiert, haben fiir uns nur die reellen Wurzeln Bedeutung, fiir die h K* > hk ist, 
was auf a 


k 
Kyh ( + cos? 6) sek? O > 6 (48) 
0 
fiihrt, wo 6 eine gewisse Zahl gréGer als 1 ist und sich zu 
ae: 
— teh [kh] (49) 


bestimmt. Mit (49) kénnen wir (48) umschreiben zu 


K,h +k hcos O> waren o* 0, 
woraus 
peo ae Neh : (50) 
ea ilies kh FP ee kh 
tgh [kh] 2 (a eh. < ] 
folet. Hier wurde zur Abkiirzung die Froudesche Tiefenzahl 
1 Ve 
Fy (51) 


Relics gh 
eingefiihrt. Setzen wir 


@, = are e]/ wets a ee : 225) (52) 
AAR ganas kh 
teh [el ls Vote —* 


so bedeutet die Ungleichung (50), daB reelle Lésungen K* => k baw. K = 0 der Gleichung (47) 
nur vorhanden sind, wenn 


O=0,. (53) 


ist. Hieraus und aus (52) entnimmt man, daBf fiir 


kh 
F, Vien — Fh < 1 


die Gleichung (47) fiir jedes O mit 0 < O < m/2 genau eine Wurzel mit K = 0 besitzt, wogegen fiir 


eh 
F.Vaeng—e Bok 


die Gleichung (47) genau eine reelle Wurzel mit K = 0 besitzt, wenn 0 im Intervall 


1 1 
are a er = Ops Ors F (54) 


tgh [k h] ages 


liegt. : ba 
Wir untersuchen noch das Auftreten anderer Wurzeln unserer Gleichung. Setzen wir K* =i K, 


so geht (47) tiber in a 
uae : — te [Kh]. (55) 
Ky, h sek? © (2 + K cos” 6) 


0 
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Zu festem O mit O < O <an/2 existiert also eine unendliche Folge rein imagindérer Wurzeln 

K* =iK, mit lim K, =oo. Die zugehérigen Werte K, bestimmen sich aus der Gleichung . 
/Ki1 #7 K, 

zu 


KS L344 Ga ee (56) 


Da K, = 0 eine Lésung von (55) ist, haben wir als rein imaginare Wurzeln von (47) mit kleinstem 
Absolutbetrag 

Ky=ttk. 
Mit den rein imaginaéren Wurzeln (56) und den vorhin bestimmten reellen Wurzeln sind alle in 
der Halbebene Re K => 0 gelegenen Wurzeln der Gleichung (47) bestimmt. 


imag. 
imag. 9 


reel! 


reel] 
Abb. 3. a) Integrationsweg fiir (x — &) cos @ + y-sin © > 0; b) Integrationsweg fiir («x — £) cos 0 + y-sin © < 0. 


In (46) betrachten wir nunmehr die Integration beziiglich K. Zunachst nehmen wir an, da 
der Integrationsweg beziiglich K in (46) nicht die positive reelle Achse sei, sondern wahlen als 
Integrationsweg den in Abb. 3 dargestellten, in der komplexen K-Ebene liegenden Weg. Die 
Polstellen des Integranden werden auf kleinen Halbkreisen umgangen, so daB der Integrand im 
Inneren der in Abb. 3 dargestellten Wege eine regulare Funktion darstellt und somit dieses Integral 
nach dem Cauchyschen Integralsatz verschwindet. Die gewahlten Integrationswege liegen hierbei 
in der unteren bzw. oberen Halbebene, je nachdem 


(x—é)cssO+ysnO>0 baw. (x—&)cosO+ ysinO < 0 
ist. 

Nun streben die Integrale langs des Viertelkreisbogens mit R—> oo gegen Null, und die Inte- 
grale langs der Halbkreisbégen um die singularen Stellen lassen sich mittels Residuenrechnung 
auswerten. Der Cauchysche Hauptwert des in (46) auftretenden Integrals beziiglich K langs der 
positiven reellen Achse ist damit gleich + zi mal Residuum des Integranden in der auf der po- © 
sitiven reellen Achse gelegenen Polstelle plus | 2i mal Summe des Residuen der auf der nega- 
tiven bzw. positiven imaginaren Achse gelegenen Polstellen des Integranden plus dem Cauchyschen 
Hauptwert (Polstellen werden durch kleine Intervalle ausgespart) des Integrals iiber die negative 
bzw. positive imagindre Achse, je nachdem (x —&) cos@+ ysin@ > 0 bzw. (x —&) cos O 
+ ysin O < 0 ist. Von diesen beiden zuletzt aufgefiihrten Bestandteilen 1aBt sich nachweisen 
daB sie fiir |x| oo gegen Null streben. Damit haben wir 


70/2 
P2l%, Ys 2) = pen | 2 O exp [— K* h] 
6 


y. Ke {cosh [K*(e + h)] — (&/K*) sinh [K*(-+h)}} {cosh [K*(¢ + h)] —(k/K*) sinh [K*(C-+h)} 
K cosh [K* h] 
K* + (Ko + k cos? @) sek? © 


x 1— (Ky + k cos? @) h sek? © sekh? [K* h] cos [y Ksin 0] cos [(«—&) K cos 0] dO +o(1) 
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fiir x—> + oo. Hierbei ist K als Funktion von O implizit durch die Gleichung 
/K? + l®? — (K, + k cos? @) sek? @ tgh [/K? + B- h] —0 
bestimmt. 


Um nun zu erreichen, daB der Gesamtdruck Pg weit vor der Quelle («-+—oo) gegen Null 
strebt, miissen wir also zu Pz nach (44) noch den Druck 


n/2 
P3(X> ¥; 2) eet oo [ es O exp [— K* h] 
6 


5 K* {cosh [K*(z-+ h)] — (k/K*) sinh [K*(z + h)}} {cosh [K*(C-+ h)] —(k/K*) sinh iKec chp 
K cosh [K* h] ; tl a ao 


K* + (Ky + k cos? @) sek? © 


* T— (Ky +k cos? ©) h sek® © sekh? [Kh] 008 Ly K sin @] cos [(x —£) K cos @] dO (57) 


addieren und erhalten damit fiir den Gesamtdruck 


Pa(%s ¥, 2) =— oe U(z) 2 ( Q expl—kV@—fP+¥+¢ =) 


An Vee + + EOF 


m/2 for) 
UG K2 K* 
— Fa | 08 6.40 | K* a exp [. K*(72+2h+)] cos [y K sin (G)} 
0 0 


n/2 
x sin [(x — £) K cos 0] dK + 2200) ic 0 dO 
0 


5 Ke {cosh .[K*(z-+-h)] — (k/K*) sinh [K*(z+h)]}} {cosh [K*(C+-h)] —(k/K*) sinh [K*(C+-h)}} 
cosh [K* h] By Saray 


0 
K* + (Ky +k cos? @) sek? © 
K—(K, + k cos? @) sek? © tgh [K* h] 


x exp [— K* h] cos [y K sin 0] 


n/2 
xX sin [(x —¢é) K cos 0] dK eu [ es O exp [ — K* h] 
0, 
yz. Ki? {eosh [K* (2-+h)] — (b/K*) sinh [K* (2 + )}} {cosh [K*(C-+-h)] — (b/K*) sinh [K* C-+h)} 


K cosh [K* hl] 
cos [y K sin @] cos [(x—&) K cos 9] dO. (58) 


: K* + (Ky + k cos? @) sek? © 
X 1 —(K, £ k cos®@) h sek? © sckh? [K* h] 


Mit (58) haben wir die Elementarlésung unseres Problems bestimmt. 


6. Die Lésung fiir das Michellsche Schiff und der Wellenwiderstand. Nach Michell}? wird ein 
Verdrangungsschiff in erster Naherung ersetzt durch eine Quell-Senken-Belegung der Mittschiffs- 
ebene. Ist S die Mittschiffsebene, sind € und ¢ die Koordinaten in der Mittschiffsebene, und ist 
n =n(&,C) die Gleichung der Schiffsoberflache,so lautet die kinematische Randbedingung an der 
Schiffsoberflache 


7) 7) te 
£(U+wgto(-ltwy=0,  firn=+n(60)- 
Durch Linearisierung ergibt sich hieraus 
0 
v(é, 0,6) =+ U (59) 


Zum anderen haben wir nach (32) fiir die Geschwindigkeits-Komponente v(x, y, z) einer Quell- 
Senken-Belegung der Mittschiffsebene S von der Dichte o(&, ¢) 


U(X» Y> 2) =a] [oes iy ie es dé de 
Ss 


\G= 2) oe Gt) 


1 J. H. Michell, Phil. Mag. 5 (1898) S. 106. 
2 JK. Lunde, Trans. SNAME 59 (1951) S. 24. 
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92) =ZaI | [ o(é,0) exp [— hk V(x—EP? +¥ +e —F] 


1 kk | 
; ae dé dC. 59’) | 
a fe —f+y¥+e—08 aa Cea aC 4. fd CP 


Der Grenzwert lim v(x, y, z) 1aBt sich mit den Hilfsmitteln aus der Theorie der Potentiale von |} 
y>0 k 
Flachenbelegungen! bestimmen, und es ergibt sich 


lim’ v(x, y, 2) = = a(x, z) . (60) | 
a 


Sind £,¢ die Koordinaten in der Langsschiffsebene, so gilt nach (60) 


(640,02) =+ Fo, 


und mit (59) ergibt sich die Dichte der unser Schiff ersetzenden Quell-Senken-Belegung zu 
on 


a(é, 6) =2 U(0) (61) 

Betrachten wir jetzt anstelle der Einzelquelle in (58) eine Quell-Senken-Belegung der Mitt- 
schiffsebene S mit der durch (61) gegebenen Dichte in einer exponentiellen Grundstrémung, so 
haben wir analog zu (58) fir den Gesamtdruck 


on “( 1 exp[—kV@—s ty + = oe 


De (Xs y; z) =—0° U(z) | | U(¢) aE Ox x Ve &P + y2 (z ay 


m/2 [oe] 

U(z 7) K? K* +k 

| U() Ft dé dt [ eos 040 | xe Rea? [—K*@ +2h+2)] 
Ss 0 0 


: 7/2 
x cos [y K sin ] sin [(x —-£) K cos 6] dK + = TO. [] U(C) FE ag dt [ cs @ dO 
0 


Ss 
(oo) 


x f K* {cosh [K*(z+h)] —(h/K*) sinh [K*(z-+h)}} {cosh [K*(C-+-h)] — (k/K*) sinh [K*(C-+h)]} 


cosh [K* h] 
0 
K* + (Ky + k cos? 0) sek? © 


* K*_(K, th cos*) sek? © gh [K*h] 


exp [— K*h] cos|y Ksin 0] sin [(x—&) K cos @) dK 


m/2 
aa U(C) SE dé de [ 08 @ exp LK" 
ch 


“Ss 


y. K* {cosh [K*(z+h)] — (k/K*) sinh [K*(z--A)}} {cosh [K*(E-+-h)] — (k/K*) sinh [K*(€--h)]} | 
K cosh [K* h] 
K* + (Ky + k cos? @) sek? © 


X T—(Kp + k cos? O)h sek? O sek h? [K*h] 


cos [y K sin 0] cos [(x —&) K cos 0] dO. (62) . | 


Hierbei ist Oy aus (52) zu entnehmen. Aus (62) erhalt man wegen (17) die Gleichung der freien | 
Beene Oe athe C(x, y) dadurch, da8 man in (62) z = 0 setzt und mit dem Faktor l/o g multi- | 
pliziert. 


Hier interessiert uns insbesondere noch ein Ausdruck fiir den Wellenwiderstand des Michelischen 
Schiffes. Nach Lunde kann man den Wellenwiderstand mittels der Formel 
oe , 7, On 7 7 . a] 
R=2 | | p00) de de’, mit y= n€',0) (63) 
Ss 


1 0. D. Kellog, Foundations of potential theory, New York 1953. 
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berechnen. Fiir p, ist der Ausdruck aus (62) einzusetzen. Macht man dies, so sieht man, daB die 
Beitrage zu R, die von den ersten drei Summanden von (62) herriihren, bei Vertauschungen der 
Integrationsvariablen &’,¢’ und &,¢ in die negativen Werte iibergehen. Die Beitrige dieser drei 
-Summanden zu R miissen also verschwinden, und der Wellenwiderstand ergibt sich zu 


n/2 
Ro 40 i ee 5 ee K*{K* + (Ky + k cos? 0) sek? ©} exp [— K* h] 
Ie 


K cosh [K* h] {1 — (K, + k cos? @) h sek? © sekh? [K* h]} dO 


%» 

| | U(c’) a {cosh [K*(¢’ + h)] — (k/K*) sinh [K*(£’ + h)]} exp [—i K cos 0 &’] dé’ dé’ 
Ss 

x | | U() a {cosh [K*(¢ + h)] — (h/K*) sinh [K*(¢ + A)]} exp [i K cos 9 £] dé d£. (64) 
Ss 


Setzen wir noch 


A+iB= i | U(C) FE {cosh [K*(é + h)] —(k/K*) sinh [K*(¢ +-h)]} exp [—i K cos OE] dE d, (65) 
Ss 


so wird 
7/2 ‘ é 
Eker 3 9, cos O K* exp [K* h] K* + (Ky + k cos? ©) sek? © 
R= 1 | (A _B ) K cosh [K* h] 1 —(K, +k cos? 0) h gek2 © sekh? [K* hl] dO. (66) 
> 


Unter Beachtung der Relation 
K* = (K, + k cos? @) sek? @ tgh [K* h] (67) 
1aBt sich (66) noch umformen zu 
m/2 


Ed, 2 ie 2 sek @ K, + k cos? © 
ser ca | a re ) K cosh? [K*h] 1 (Ky -++ k cos? @) h sek? © sekh? [K* h] dO (68) 
% 


oder 
n/2 


== sek © K**tgh [K* hl] (Ky + k cos? 6) 
ee i ae ~ K {sinh [2 K* h] — 2 K* hy dO. (69) 
9% 


Andere Darstellungen fiir R lassen sich noch angeben, wenn man mittels (67) anstelle von @ die 
Integrationsveranderliche K einfiihrt. Hierauf kommen wir spater zuriick. 


7. Die Lésung fiir das Hognersche Drucksystem. Wir betrachten hier ein auf den Bereich B 


der freien Wasseroberflache wirkendes Drucksystem p(x, y). Wir wollen hierbei annehmen, daB 
der vorgegebene Druck als Fourierintegral dargestellt werden kann: 


Posy) =F | [ PE) i dy | 40 [op [—i K {(x—£) cos O + (y —n) sin O}] KdK. (70) 
B —2 0 


- Die in den Abschnitten 2 und 3 hergeleiteten Differentialgleichungen und Randbedingungen 
gelten natiirlich auch hier. Fiir den Druck p(x, y, 2) werden wir damit auf den Ansatz 


poy) =P | [ pln) dé dy [46 f exp [—i K{(e—8) o08 6 + (y 1) sin BF 
B —2 0 


x {F,(0, K) cosh [K* (z + h)] + F,(O, K) sinh [K* (2 + h)]}} dK (71) 
gefiihrt. 


Hierbei ist wieder 
K*=/K?+2 und Uz) =Ve*, 
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Die Randbedingung am Tankboden 
ép 
Oz 


= 0 


z=—h 


liefert nun 


FO, K) =— = F(9,K), (72) 


und mit F,(0, K) = F(0, K) haben wir damit den Ansatz 


4 70? 


p(x, y, 2) =2@ 29 | | P(E.) dé dy [ dO [oo [—i K {(x —£) cos O + (y —7) sin OF] 


x F(O, K) {cosh [K*(z + h)] — a sinh [K*(e + h)}} dK. (73) 
Nach (12) haben wir als Randbedingung an der freien Oberflache 


kee a2 {OS ieee ae fir z=0. 


Mit (70) und (73) ergibt sich daraus 


K* K 
F(O,.K) = — San [K* bh] K* 22 eee one R* hy 


und somit 
4 


Pe ys2) = 89 | | PEs) a dr | 
E 


Fe, 


16 | exp [—i K {(x — &) cos 0 + (y — 7m) sin OF} 


K K* __ cosh [K*(z + h)] —(&/K*) sinh [K*(z + h)] yy 
cosh [K*h] | K* — (K, + k cos? ©) sek? © tgh [K* h] é 


x 


Der Ubergang zur reellen Darstellung liefert 


n/2 foe) 
Se UGyeiit K K* {cosh [K*(z + h)] — (k/K*) sinh [K*(z + h)}} 
P(% J. 2) = —S—— al P(E, 1) dé dy | dO | cosh [K* h] 
B 0 0 


cos [(y —”) K sin 6] cos [(x —£) K cos 6] 
‘ tS (Ky + k cos? ©) sek? © - tgh [K* hy dK. (74) 


Zur Umformung des Integrationsweges beziiglich K geht man wie in Abschnitt 6 von der Dar- 
stellung 


7/2 foe) 
_ _ eU@) ke K K* {cosh [K* (z + h)] — (k/K*) sinh [K* (z + h)}} 
BUS YS) ia | | P(S.1) dé dy | dO | cosh [K* h] {K* —(K, +k cos? ©) sek? 6 igh [K* h]} 
B 0 0 


x + fexp [—i K {(x —) cos © + (y —7) sin OF] 
+ exp [—1 K {(x — &) cos 0 — (y — 7) sin O}]} dK. 


aus. Andert man wie in Abschnitt 6 den Integrationsweg, so 1aBt sich genau wie dort zeigen, 
daB fiir x ++ oo 


m/2 
wail UO fiat K*® {cosh [K* (z + h)] — (k/K*) sinh [K* (z + h)} 
1 AEE eats th | | PIS, 0) dE dy |e [K*h] {1 — (Ky + k cos? 6) h sek? © sekh®[K* h]} 
B 0 


0 


_ X cos [(y —7) K sin 0] sin [(x — &) K cos O] dO + o(1) 


wird. Damit die Wellen weit vor dem Drucksystem aufgehoben werden, hat man also zu P(x, ¥, 2) 
nach (74) den Druck 


me] 2 
_eU@ ff. K* {cosh [K* (z + h)] — (kK*) sinh [K* (z + h)] 
Pity 2) 7 i) | P(E, 4) dE dn | oe [K* h] (1 — (Ky + k cos? O)h sek? © sekh? [K* 5 
(0) 


° 


x cos [(y — 1) K sin O] sin [(x — £) K cos 0] dO (75 
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zu addieren, und erhalt somit als Lésung unseres Problems 


: eyroe Eee U(z) a r oK K* {cosh [K*(z + h)] — (k/K*) sinh [K*(z-++ A] 
Pele») ne ig: Beal) a i ee cosh [K* h] {K* — (K, + k cos? ©) sek? © ee IK A 


x cos [(y 7) K sin 6] cos eee K cos 0] dK +25 | [Be n) dé dn 
B 


[2 
x K™ {cosh [K*(z + h)] — (k/K*) sinh [K*(z + h)} 
cosh [K* h] (1 — (K, + k cos? @) h sek? © sekh? [K* h]} 


Oo 
x cos [(y —7) K sin 0] sin [(x — &) K cos 0] dO. (76) 


Die Gleichung der freien Wasseroberflache ist nach (15) zu berechnen. Der Wellenwiderstand 
ist hier nach der Formel 


Re | i PEs 1!) ge AE! dy (77) 
B 


zu bestimmen, wo 0C/0x sich aus (76) mittels 


Bisel «[Op- op 
fae + a 


berechnet. Mithin wird 


: 1 — (k/K*) teh [K* h] 
= = ay =| | 7 Ntedy om 040 0K —(K, + k cos? ©) sek? © tgh [K* hi] 


7/2 
X cos [(y—n) K sin O] sin [(x—é) K cos 6] aK +z, | [ Pe y) dé dy [ cesOKK 
B 
1 — (k/K*) teh [K*h é 
XG area Oo rere yo [(y—n) K sin @] cos [(x—&) K cos 9] dO 


7/2 


Sag] | PE a a | cos 00 


x | K? cos [(y — 7) K sin 0] sin [(« — £) K cos O] dK. (78) 


Setzt man dies in (77) ein, so sieht man, da die Beitraige zu R, die vom ersten und dritten Sum- 
manden von (78) herriithren, bei Vertauschung von é’, 7’ und &, 7 in die entsprechenden negativen 
Werte iibergehen. Die Beitrage dieser beiden Summanden zu R miissen also verschwinden, und fiir 
den Wellenwiderstand erhalten wir: 


n/2 
_V [cos @ K K* (1 — (b/K*) teh [K* h}) 
R=, [UK Fh coe Gy h sok O sek? [K* h] dof |e nce 7) 


9 


xX exp [—i K (&’ cos 9 +7 sin @)] d& dr’ 7 P(E, 1) exp [i K ( cos O+-n sin O)] dE dy. (79) 
B 


Fihrt man zur Abkiirzung 


P= | { peésn) dé dy. C+iD=>{ | Pn) exp[—iK (Geos O +1 sin 6] dd (80) 
B z B : 


ein, so wird 


cos © K K** (1 —(k/K*) tgh [K* h])_ ag 


ee. j (Oe Le cen —(K, Ek cos? ©) h sek? © sekh? [K* h] (81) 


a. +) 


Ingenieur-Archiv 
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Beachtet man noch 
K* =(K, + k cos? @) sek? 0 tgh [K* h] , (82) 
so kann man hierfiir schreiben 
n/2 


Ee 2 sek @ KK* tgh [K* h] 

= 2 2 te ek O KIN La ee 

or t D) T=, Ecos 6) sek? Osckb [K* i] 0. (83) 
Gy 


Bericksichtigt man schlieBlich, daB wegen (82) 
(Ky +k cos? @) h sek? @ sekh? [K* h] =" * * _ 
nd ~ ‘sinh [2 K* h] 


gilt, so erhalt man 


m/2 
eR ge i . 5, sek @ K K* sinh? [K* h] 
R= | (C+ DY ge a Keg 19 (84) 
OQ 


Durch Einfiihrung der Integrationsveranderlichen K an Stelle von @ kénnen noch andere Dar- 


stellungen fir R angegeben werden. 


8. Lésung fiir den ebenen Fall des Hognerschen Drucksystems. Hier entfallt die Abhangigkeit 
von y und die Gleichungen (5) bis (8) reduzieren sich auf 


dU 

(ha # ‘dz op 
i code Sie (85) 

aU x 

1 dz 
o =varialy alee ZA. (86) 
me Wage | 

© = — G5 | ihe (87) 


In (85) ist die rechte Seite Null, da im Strémungsbereich — h < z < 0 keine Quellen liegen sollen. 
Die Randbedingung (12) gilt unverandert. 


Geht man aus der Darstellung des auf B vorgegebenen Druckes in der Form 
2. O Viele a J 
Po) = Sz | PE) a [ exp [—i K wg) aK, (88) | 
B f 
so wird man fiir p(x, z) auf den Ansatz 


U(z = - 
pee 2) = | pie) ag | exp [—i K (x —8)] {F,(0, K) cosh [K*(z + h)] 
B 0 
+-F,(@, K) sinh [K*(z +h)]} dK (89) II 
gefiihrt, wobei 2 ] 
KARE und.) Ue) a 
ist. Die Randbedingung am Tankboden 


gee 
ORNs pe 


ergibt 
F(0, K) =— 2; F(@, K). 
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Mit F,(0, K) = F(@, K) hat man also fiir p(x, 2) den Ansatz 


p(x, 2) =e) x | Pe) P(E) dé ij exp [—i K (x —8)] F(O, K) {cosh [K*(z + h) 


B 
— (k/K*) sinh [K*(z + h)]} dK, (90) 
und aus der Randbedingung an der freien Oberflaiche 
op gop 2 dp = 
ax? EZ V2 @z dx?” AONE 
ergibt sich mit (88) und (90) 
es 
F = — = Ss 
(0, K) cosh [K* h] {K* — (Ky + b) tgh [K* hp} ’ Ko = Tr 
und somit 
ee eG) K+ {cosh [K* (z+ h)] —(k/K*) sinh [K*(¢+h)]} 
P(*, 2) 2 ) ds i couh{K* AJ (K* —(K, £h) eh [KPa eh EE) dK. 91) 


Der in (91) auftretende Integrand hat nun Polstellen, wenn 


= (Ky + k) tgh [K* hl] 


‘ ist. Schreibt man diese Gleichung in der Form 


K*h 
Gre es 
kh 


so sieht man, daB der Integrand von (91) fiir (K, +k) h < the bl 


K-Achse liegende Polstelle hat, wogegen fiir 
(Kj +h) h> 


keine auf der positiven reellen 


kh 
tgh [k h] 


stets eine und nur eine positive reelle Polstelle des Integranden existiert. Ist also 


kh 


und andert man den in (91) auftretenden Integrationsweg wie in Abschnitt 6, so lat sich zeigen, 


daB fir x +--+ co 


P(x z) 2m ae ert [© K* {cosh [K*(z ate h)] at (k/.K*) sinh [K*(z ale h)}} sin [(x wat é) K] dé oe o(1) 


K cosh [K* h] {1 — (Ky + hy h sekh? [K* h]} 


ist. Damit die Wellen weit vor dem Drucksystem aufgehoben werden, haben wir also zu p(x, 2) 


_ nach (91) den Druck 


ue K*® {cosh [K*(z + h)] — (k/K*) sinh [K*(@ +h) [apn 
Pils 2) =" K cosh [K* h] {1 — (Ky + ht) h sekh? [K* h]} Bee eam) ede. (22) 


B 


gu addieren. Somit lautet hier die Lésung unseres Problems in reeller Darstellung: 


ist zu berechnen nach der Formel 


A ies 
af 


se K* {cosh [K*(z + h)] — (k/K*) sinh [K*(z + A)]} 
Pelt 2) =—"5 f? (6) a6 | Cae? APU =2(Ky + eh [Keay A TR 


0 U(2) K* {cosh [K* (z + h)] — (&/K*) sinh [K* ( + bY} (are ee 
58 K cosh [K* h] {1 — (Ky - ky) h sekh? [K* h}} p(é) sin [(x —¢) K] dE. (93) 


Hieraus ergibt sich die Gleichung fiir die freie Wasseroberflache mittels (13). Der Wellenwiderstand 


i [Pear (94) 
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wobei sich dé/dx wegen (88) und (93) aus 
dé 1 


ergibt zu 


d ees 7 KK*(—(&/K*)tgh[K*h} . 
- Ing | PS) 4. | Ke eek Ke es) Elis 
B 0 
V K*(1—(k/K*) tgh [K* A} ( — S 
Y 2g 1—(Ky +h) h sekh? [K* h] P(E) cos [(% —£).K] ds 
B 


2g 
B 


ang | PO) ag { Ksin [(x—é) K] dK. 
9 | 

LaBt man wie friiher die Bestandteile fort, die bei zweifacher Integration tiber die Drucklinie 
den Beitrag Null ergeben, so erhalt man 
V K*(1—(k/K*) tgh[K*h} ( —/¢, me SKE) dé. 95) 

a Fel = (Ko k) h sekh? [K* h] p(é’) exp [—# &'] dé p(&) exp [ ] ( 

B B | 

Mit den Abkiirzungen 


P= | pede, C+iD=> [ Peep l-iK ede (96) 


B 


ergibt sich aus (95) 
VP? K*. {1 — (k/K*) tgh[K*h]} || 
R= oy (CO + DY T= (K+ he sokh® [K* hl] ” a | 


falls 
kh 
(Ky sb k) h > tgh [k h] 
ist. Beachtet man noch 
K* =(K, + &) tgh [K* A], (98) 


so kénnen wir (97) in der Form 
K* tgh [K* h] 
(Ko Lk) h sekh? [K* h] 


R=s, (C+ DY) — (99) 


schreiben. Wird schlieBlich beriicksichtigt, daB wegen (98) 


2K*h 
sinh [2 K* h] 


(Ky +k) hsekh? [K* h] = 


gilt, so erhalt man 
K* sinh? [K* h] 
sinh [2 K* h] —2 K*h * 


R=* (C+ D2) (100) | 


Damit sind fiir das dreidimensionale und zweidimensionale Problem die Formeln fiir den Wel- 
lenwiderstand abgeleitet. Die Berechnung des Wellenwiderstandes nach diesen Formeln und ein 
Vergleich mit Versuchsergebnissen, wie auch eine Untersuchung des Wellenfeldes in groBem Ab- | 
stand vor und hinter dem Schiff, wird Gegenstand einer spateren Arbeit sein. 


(Eingegangen am 6. April 1960.) 


Anschrift des Verfassers: Dr. F. Kolberg, Laurensberg b. Aachen, Roermonderstr. 2 f. 
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Ein einfaches Rechenschema zur Aufstellung der Frequenzgleichung 
eines an den Enden freien Drehschwingers 


Von R. Priebs* 


1. Einleitung. Im folgenden wird gezeigt, wie die Frequenzgleichung eines an den Enden freien 
Drehschwingers als Polynom des Frequenzquadrates z = w? sehr iibersichtlich und schnell aufge- 
stellt werden kann. Das Wesentliche besteht dabei in der Herleitung von Rekursionsformeln zwi- 
schen den Koeffizienten der Frequenzgleichungen eines aus n Drehmassen und eines aus (n — 1) 
Drehmassen bestehenden Schwingers. Diese Rekursionsbeziehungen erlauben die Aufstellung eines 
fur die numerische Berechnung der Frequenzgleichung sehr praktischen Rechenschemas. Die 
Lésung der Frequenzgleichung kann nach dem Horner-Newtonschen Verfahren erfolgen (Abschn. 4). 

Das Vorgehen ist eng verwandt mit dem Verfahren von W. A. Tuplin, bietet jedoch rechen- 
technische und zeitliche Vorteile gegeniiber jenem (Abschn. 6). 

Hinsichtlich des Zeit- bzw. Rechenaufwandes diirfte das hier beschriebene Vorgehen insbeson- 
dere dann dem Verfahren von Holzer-Tolle iiberlegen sein, wenn keine Naherungswerte fiir die 

_Kigenfrequenzen vorliegen und mehr als eine Eigenfrequenz gesucht wird. Zur allgemeinen Orien- 
tierung tiber alle wesentlichen hierher gehérenden Verfahren sei auf den umfassenden Bericht von 
K. Klotter ,,Analyse der verschiedenen Verfahren zur Berechnung der Torsionseigenschwingungen 
von Maschinenwellen“+ verwiesen. Hinweise auf diese Arbeit werden im folgenden durch ,,Bericht 
Klotter, . . .“* gegeben. = 


2. Grundlagen. Unseren Betrachtungen liegt ein Drehschwinger nach Abb. 1 zugrunde. Er 
besteht aus n + 1 Scheiben mit den Massentragheitsmomenten 0, (A = 0,1, 2,...,7), die ela- 
stisch verbunden sind durch n glatte, zylindrische, tragheitslos angenommene Wellenstiicke mit den 
Torsionssteifigkeiten c,(A = 1, 2,...,). Der Winkel (A = 0,1, 2,...,) beschreibt die Ver- 
drehung der Masse @,. In dem Wellenstiick mit der Torsionssteifigkeit c, ist das Drillungsmoment 
- D, wirksam. Bei der Berechnung wird jegliche Dampfung vernachlassigt. 


G Gn-1 
Oy. 
b au By 
Re Cy ae Cn 
O, 
G, Q,, % 


G, 


Abb. 1. Gebrauchliche Darstellung eines an den Enden freien Drehschwingers. 


O,-1 


Die Bewegungsgleichungen des Drehschwingers sind gegeben durch die dynamischen Glei- 
_ chungen 


0,3, = D,—Djr41 (1) 
und die elastischen Gleichungen 
D, = & (,—1 — 94). (2) 
Mit dem harmonischen Ansatz fiir den Verdrehungswinkel 
0, =u, cosmt (3.1) 
und das Drillungsmoment 
; D, = X%, cos wt (3.2) 
erhalt man aus (1) und (2), wenn 
2 CDs (4) 


i i j i - Darm- 
* Verfasser ist den Herren Prof. Dr. phil. P. Matthieu (TU Berlin), Prof. Dr. Ing. K. Klotter (TH 
Peau Diplelng, H. Pfiitzner (TU Berlin) fiir wertvolle Hinweise sehr zu Dank verpflichtet. 
1 Ing.-Arch. 17 (1949) S. 1. 
| 10 


142 R. Priebs: Ein einfaches Rechenschema zur Aufstellung der Frequenzgleichung _Ingenieur-Archiv 


gesetzt wird, die Gleichungen 
XA 
41 =%, + O,u,z, u,= sy eae Fs . (5) 


Die x, bezeichnen gemaf (3.2) die Amplituden der in den Wellenstiicken wirksamen Drillungs- 


momente, die u, gema8 (3.1) die Amplituden der Ausschlage der Drehmassen. 
Da der Drehschwinger an den Enden frei ist, ergeben sich die Bedingungen 


0 (6.1) 


Xn+1 Ss 0 . (6.2) 


und 


Ferner.ist ¢) = ¢, 11 = 0- 
Aus (5) erhalt man somit ausgeschrieben fir 7 = 0,1,2,...,n die folgenden Gleichungen: 


Xy — 0, Up Zz » 

aT 
uy, =— Ug =r teay! 
A =a Opie, 

Xp 
Ug => Uy e 

(7) 

= 1+ Opa 23 

xy 
uj; gy a 


Xn4+1 = X ay (oF Uy, Zz - 

Die sukzessive Durchrechnung der Gleichungen (7) fiir einen festen Zahlenwert von z = w? 
liefert ein ,,Restmoment“ x,,,;. Berechnet man sich mehrere solcher Restmomente und tragt diese 
iiber z auf, so ergeben sich die Eigenfrequenzen aus den Schnittpunkten dieser Kurve mit der 
z-Achse; denn nur fiir diese Werte ist (6.2) erfiillt. Dieses Vorgehen ist als Verfahren von Holzer- 
Tolle bekannt. Fir die Durchfihrung dieser impliziten Methode sind einige fiir die numerische 
Rechnung zweckmaBige Rechenschemata entwickelt worden. Ein groBer Teil der Rechnung kann 
auch durch zeichnerisches Vorgehen ersetzt werden (Verfahren von Giimbel-Geiger)'. 

W. A. Tuplin benutzt die Gleichungen (5) zur expliziten Aufstellung der F: requenzgleichung, 
der gleiche Gedanke liegt auch dem folgenden zugrunde. 


3. Das Rechenschema. Ersetzt man in den Gleichungen (7), mit A = 1 beginnend, nacheinander 
die x,,1 und die u, durch die vorangegangenen x, und u,_j, so erhalt man Polynome in 2, deren 
Koeffizienten nur noch uy und die @, und c, enthalten. Fiir die Polynome, die sich fiir die Drill- 
momentenamplituden ergeben, kann geschrieben werden 

x, = Ul Ay, 2) » 
%_ = Up( Ap,  — Ago 2’), 
Xg = Ug (Ag, z — Agg 2” + Ags 2°) , 


(8) i] 


My ==tu( Ag sz — Ayn 2* + ---—(—1)* Aya), 
Sint == Ug Anii,1 2 —Anri.2 ¥ +. — (— It A ie) 


Die Polynome, die sich fiir die Amplituden der Ausschlaige ergeben, sind darstellbar durch 
Uy = uo(l — Bg; z), 
Uz = uo(l — Bs, 2 + Bsp 2), 
Us = Uo(l — By, z + Byo? — By 2), | 
(9) |] 
Ua = U(l — Byyi1 2 + Bayo,22®—--- + (—1)* By, 24), 


u, = Up(1 —— Buy, 12% + Bi 41,2 a oe + (—1)" Bene Z") ° 
1 Vgl. Bericht Klotter, S. 12—16. 
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Die Indizes der Koeffizienten A bzw. B wurden so festgelegt, daB der erste Index auf die Anzahl 
der Drehmassen hinweist und der zweite auf die zugehérige Potenz von z. 

Nun benutzen wir (8) und (9) bei der sukzessiven Durchrechnung der Gleichungen (7), wobei 
wir den Parameter z mitfiihren. Es ergeben sich dabei ganz einfache Rekursionsbeziehungen 
zwischen den Koeffizienten, die durch eine besondere Umrahmung gekennzeichnet sind: 


Ng Or, e = Uo(Ay 2), 
also 
Ay, = % |; (10.1) 
Uy = Uy — = Uy (1 — Bz), 
A 
also 
By = Sa ; (10.2) 
z | 
X= % + O,u,z = Ug (Ag, = — Age 2”) 5 
= Uy (A), 2 + O; (1 — By, 2) 2), 
= Up (An te 0,) re 0, By, 2?) ’ 
also 
(10.3) 
21% —— (Ag, 2 — Ap, *)) 2 
Co y 
also 
(10.4) 
As, = Ag; a O, 2 
Az, = Ayo + O Bs, | ; (10.5) 
“A — 0, Bo 
(10.6) 
Any = As, - (OR 9 
Aga = Age + O3 Bar, | | (10.7) 


Ax, — Ags + O; By, 
Aas = Os Bas 


10* 
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Das Bildungsgesetz fiir weitere Koeffizienten ist sofort erkennbar. Es kann leicht durch vollstandige 
Induktion bewiesen werden. 

Die aus den Rekursionsformeln (10.1) bis (10.7) ersichtlichen GesetzmaBigkeiten kann man 
durch das Koeffizientenschema Tabelle 1 sehr iibersichtlich darstellen. In der untersten Zeile 
dieses Schemas sind die zu den Koeffizienten einer Spalte gehérenden Potenzen von z und das je- 
weilige Vorzeichen zur Erleichterung noch angegeben. Auferdem soll noch bemerkt werden, daB 
die Addition der in einem Kastchen stehenden beiden Gréfen die unter dem Kastchen stehende 
Gre ergibt. 


Tabelle 1. Das Koeffizientenschema. 


Das Koeffizientenschema kann ohne weiteres fiir die numerische Berechnung der Koeffizienten 
herangezogen werden. Es ist jedoch vorteilhaft, wenn die dabei erforderlichen Operationen mit 
dimensionslosen GréBen durchgefihrt werden kénnen. Zu diesem Zweck schreiben wir die ein- 
zelnen Massentragheitsmomente und Torsionssteifigkeiten als Vielfache eines Bezugs-Massentrag- 
heitsmomentes @* baw. einer Bezugs-Torsionssteifigkeit c*: 


0, = 0, O*, =, ¢*. (11.1) 
Es ist also @, das OQ;-fache von O* und c, das c)-fache von c*. Durch Wahl der Bezugsgréfen O* 
und c* kénnen wir dann geeignete dimensionslose GréBen 


eo = & 


a = ee baw. 4 = = (11.2) 


c* 
bilden. 


Denkt man sich nun das Koeffizientenschema mit den Werten 0, @* und cj c* durchgerechnet, 
so erkennt man, daf jedes Glied der ersten Spalte den Faktor @*, jedes Glied der zweiten Spalte 
den Faktor @*/c*, jedes Glied der dritten, vierten, fiinften, ... Spalte jeweils den Faktor 9*2/c*, 
O*/c*2, O*3/c™, ... enthalt. Diese gemeinsamen Faktoren jedes Gliedes einer Spalte kann man 
herausziehen und mit den jeweiligen Potenzen von z (vgl. Koeffizientenschema) zusammenfassen. 
Mit den Abkiirzungen 


caf: (12.1) 
und 
p=c*l=O*:% (12.2) 


gelangt man dann zu dem Rechenschema von Tabelle 2. 


Man erkennt, da die Durchrechnung des Rechenschemas mit den Zahlenwerten 0; und cj die 
x, und uy jetzt wie folgt liefert: 


%4(C) = Up p (Agi — Ajo +-+-—(— 1 Aj, -}), 
u,(2) = uy (1 — Biyiig +---+(—1) Biyial). 


Bei der numerischen Berechnung der 4;, und B;, nach den durch das Rechenschema gegebenen 
Vorschriften erweist es sich als vorteilhaft, daB haufig Multiplikationen bzw. Divisionen mit der 
gleichen GréBe durchzufiihren sind. 


(13) 
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In praktischen Fallen gentigt es, ein leeres Formblatt etwa in der Art von Tabelle 3 zu benutzen. 


Beim Ausfiillen des Formblattes wird es nur zu Anfang nétig sein, das Rechenschema daneben 
VAD legen, um sicherer rechnen zu kénnen. Ubersichtlichkeit und mechanischer Ablauf der Rechen- 
operationen pragen sich schnell ein. 


Tabelle 2. Das Rechenschema. 


3 


AL 
/ Uf SKS er ed) 
Ays Fis = on Ajy=O, 845 


« | | ! 
“1 ) | | Use A \ | | 
| l ! 


vipa es Sse eet ol a5 gee fs 


Tabelle 3. Formblatt fiir die 
Anwendung des Rechenschemas. 


4, Aufstellung und Lésung der Frequenzgleichung. Die Frequenzgleichung ergibt sich (vgl. 
Abschn. 2) aus der Bedingung x, ,; = 0. Mit Hilfe des Rechenschemas erhalt man sie in-der Form 


F,41(¢) = “netts) = Any1,1— Angie 6 +++» —(— 1)! Ansan ot = 0. (14) 
Hier ist bereits der die triviale Lésung €, = 0 enthaltende gemeinsame Faktor p = c* € aller 
Koeffizienten abgespalten. Aus ihren nichttrivialen Lésungen €, (i = 1, 2,..., 7) ergeben sich die 
Eigenfrequenzen w; gemaB (4) und (12.1) zu 


wo, =~% = eee (15) 


Die Ermittlung der Nullstellen ¢; wird man im allgemeinen am schnellsten durch die Kombination 
von Hornerschema und Newtonscher Wurzelverbesserung! erreichen. Als Naherungswert fiir die 
(it + 1)-te Wurzel kann dabei die vorangegangene i-te Wurzel benutzt werden, wobei man mit der 
bekannten trivialen Lésung z = 0 baw. €, = 0 beginnt und natiirlich jedesmal mit der so gefun- 
denen Wurzel reduziert, sobald diese in der erforderlichen Genauigkeit vorliegt. Die Reduktion ist 
bei Benutzung des Hornerschemas sehr einfach! Durch dieses systematische Vorgehen ist es nicht, 
wie z. B. bei Holzer- Tolle, notwendig, durch Probieren die Nullstellen aufzusuchen, sondern die gute 
Konvergenz des Verfahrens fiihrt, ausgehend von der trivialen Lésung, rasch zur Grundfrequenz, 
von dieser wiederum zur Higenfrequenz der ersten Oberschwingung usw. 


Die Aufstellung der Frequenzgleichung ist besonders einfach, wenn viele der Drehmassen 0, 
und Federzahlen ¢, den gleichen Wert O baw. c haben. Wird namlich O* = O und c* = ¢ gesetzt, 
dann folgt @’ =c’ =1, s0 daf alle im Rechenschema mit diesen Werten auszufiihrenden Divi- 


1 Vel. R. Zurmiihl, Praktische Mathematik fiir Ingenieure und Physiker, S. 42ff., 2. Aufl., Berlin/Gottingen/ 
Heidelberg 1957. 
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sionen und Multiplikationen mit dem Divisor baw. Faktor | durchzufiihren sind. Die Durchrech- 

nung der vollkommen ,,homogenen™ Maschine (O, = O und c, = ¢) fibrt selbstverstandlich auf 
° 7 (74 

die von R. Grammel eingefiihrten ,,reduzierten Frequenzfunktionen“?. 


5. Ein Zahlenbeispiel. Es sind die Eigenfrequenzen des durch folgende Daten gegebenen Dreh- 


schwingers” zu bestimmen: 


OQ, = 0,36 cm kg sek’, ¢, =— ‘10 emkg, 


@, = 0,15 cm kg sek’, C = —:10° emkg, 


0, = 0,21 cm kg sek’, é==-—~2 102 cuikg., 
0, = 0,09 cm kg sek” . 


Wird [vgl. (11.1) und (11.2)] 
O* —1lcmkgsek® und c* = 10° cmkg 


gewahlt, so ergibt sich das Rechenschema von Tabelle 4 (vgl. Tabellen 2 und 3). 


Tabelle 4. Ausgefiilltes Rechenschema 
zum Zahlenbeispiel. 


462 0.3402 
N32 | 41647 


12,9417 45049 


2,5326 


Man erhilt also die Frequenzgleichung 


F6) = 2D = Ay — Aig + Aig O— A = 9, 


— 0,81 — 2,5326 £ + 1,5049 C2 — 0,21433 3 — 0 


oder 
C3 — 7,0214 ¢2 + 11,816 € — 3,7792 = 0. (16) 


Die Lésung nach dem Horner-Newtonschen Verfahren (vgl. Abschn. 4) liefert €, = 0,4171. 
Nach der Reduktion von (16) mit diesem Wert verbleibt die quadratische Gleichung 


6? — 6,6043 € + 9,061 — 0 
mit den Wurzeln ¢, = 1,944 und ¢, = 4,660. Die Eigenfrequenzen ergeben sich gem48 (15) zu 


TET eee oe ee es 
oO; = \é ar = 6," 10° sek, 
@ = 4,171 +104 sek = 204,2 sek2, — w, = /19,44- 104 sek“! — 440,9 sek, 
Ws = 46,60 - 104 sek? = 682,6 sek. 


6. Weitere Anwendungsméglichkeiten des Rechenschemas; Vergleich mit dem Verfahren yon 
W. A. Tuplin. Wir sahen [vgl. (13)], daB dem Rechenschema alle x, und u, unmittelbar entnommen 
werden kénnen. Fiir den im Zahlenbeispiel behandelten Schwinger entnimmt man dem zugehérigen 
Rechenschema (Tabelle 4) z. B. die Ausschlagsamplitude der Drehmasse O, zu 


Up = Up (1 — Bs, 6 + Bs. 0), 
Uy = Up (1 — 5,13 ¢ + 1,62 22). 


ae C. B. Biezeno u. R. Grammel, Technische Dynamik, Bd. 2, S. 366, 2. Aufl., Berlin/Gottingen/Heidelberg 


2 Vel. Bericht Klotter, S. 13. 
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Mittels Hornerschema erhalt man bei der ersten Oberschwingung wegen ¢, = 1,944 die Amplitude 
u(C,) = — uy: 2,6930 . 
Auf diesem Wege kann man sich leicht alle (auf u, bezogenen) Ausschlagbilder verschaffen. 


Wenn die Bestimmung der Eigenfrequenzen und damit die Bestimmung der Frequenzgleichung 


auch die wichtigste Aufgabe ist, so sollte doch gezeigt werden, daB die Ermittlung der Ausschlag- 
bilder mit Hilfe des Rechenschemas sehr einfach ist. 


Noch wertvoller erweist sich das leichte Auffinden aller x,(¢) und u,(C) bei der Behandlung 
verzweigter Drehschwingungssysteme!, insbesondere dann, wenn diese Systeme viele » inhomogene“ 
Teile enthalten. Darauf soll hier jedoch nicht naher eingegangen werden. 


Kin Vergleich des Verfahrens von Tuplin mit dem hier beschriebenen Vorgehen an Hand ein- 
facher Zahlenbeispiele® zeigt, da hinsichtlich der Anzahl der Rechenoperationen keine nennens- 
werten Unterschiede bestehen. Die Durchfiihrung der Rechenoperationen ist bei Tuplin jedoch 
etwas umstandlicher und zeitraubender: Aus den Daten des Drehschwingers miissen mehrmals fiir 
die Rechnung geeignete Konstanten gebildet werden, wahrend hier sofort mit den Daten gerechnet 
werden kann, und auferdem jeder einzelne Wert nur einmal fiir eine Reihe aufeinander folgender 
Divisionen oder Multiplikationen gebraucht wird. Auch ist es fiir das numerische Rechnen vielfach 
stérend, wenn drei Zahlen mit verschiedenen Vorzeichen zu addieren sind, wie dies bei Tuplin der 
Fall ist, nicht aber bei Anwendung des hier beschriebenen Rechenschemas. Wollte man als Ver- 
gleichsmaBstab die Anzahl der anzuschreibenden Zahlen heranziehen, so wiirde das Verfahren von 
Tuplin ebenfalls schlechter abschneiden. 


Als sehr vorteilhaft gegeniiber dem Verfahren von Tuplin erweist sich das hier beschriebene 
Vorgehen in den zuvor angedeuteten Fallen, in denen die Ausschlagsamplituden oder die Drill- 
momentenamplituden in ihrer expliziten Form u,(z) baw. x,(z) benétigt werden: Sie kénnen dem 
Rechenschema sofort entnommen werden, wahrend fiir ihre Ermittlung nach Tuplin eine Reihe 
von Umrechnungen erforderlich ware (s. auch Bericht Klotter, S. 33—34). 


SchlieBlich sei noch bemerkt, daB sich fiir die Behandlung beidseitig eingespannter Dreh- 
schwinger (analog zu Abschn. 3) ebenfalls ein einfaches Rechenschema aufstellen lat. An die Stelle 
der Randbedingungen x) = 0 und x,,, = 0 treten dann die Bedingungen u, = 0 und u, = 0. Die 
Frequenzgleichung eines einseitig eingespannten Drehschwingers (x) = 0, u, = 0) kann mit dem 
hier abgeleiteten Rechenschema aus der Bedingung u, = 0 ermittelt werden. | 


7. AbschlieBende Bemerkungen. Zusammenfassend kann festgestellt werden, daB das be- 
schriebene Vorgehen besonders fiir die Berechnung inhomogener Drehschwinger geeignet ist. Dies 
liegt in erster Linie daran, da8 man die Frequenzgleichung mit geringem Aufwand explizit berechnen 
und auch beliebig genau lésen kann (s. Abschn. 4). 


Einen interessanten Hinweis bekam der Verfasser von Herrn Prof. Dr.-Ing. L. Cremer (TU 
Berlin): Die umgekehrte Problemstellung der vorliegenden Arbeit hat man in der Regelungstechnik, 
namlich die Bestimmung der zu einer vorgegebenen Frequenzgleichung gehorigen Schwingungskette 
mit ihren einzelnen Gliedern. Die Lésung dieses Problems ist mit dem bekannten Routhschen 
Schema’ méglich. Das hier abgeleitete Schema 1aBt sich also als eine Umkehrung des Routhschen 


Schemas deuten. 


Mitteilung aus dem Heinrich-Hertz-Institut fiir Schwingungsforschung Berlin-Charlottenburg. 


(Eingegangen am 25. April 1960.) 


Anschrift des Verfassers: cand. ing. Ralf Priebs, Heinrich-Hertz-Institut, Abteilung Mechanik, 
Berlin-Charlottenburg 2, JebensstraBe 1. 


1 Biezeno-Grammel, Technische Dynamik, Bd. 2, S. 385, 2. Aufl. ; 

2 Miele Bojesielin Bericht Klotter, S. 20 und Beispiel in Abschn. 5 der vorliegenden Arbeit. A 

3 Siehe hierzu: E. J. Routh, Die Dynamik der Systeme starrer Korper, Bad. 2, S. 222ff., Leipzig 1898; 
W. Cauer, Theorie der linearen Wechselstromschaltungen, Bd. 1, S. 196ff., Leipzig 1941; L. Cremer, Z. Rege- 


lungstechn., 1 (1953) S. 38. 


148 F. Manna: Wirkungslinie des Massenwiderstandes eben bewegter Scheiben Ingenieur-Archiv 


Zur Ermittlung der Wirkungslinie des Massenwiderstandes eben bewegter 
Scheiben und Getriebeglieder 


Von F. Manna 


1, Einleitung. Zur Ermittlung der Wirkungslinie der resultierenden Tragheitskraft P = — m bg 
eines eben bewegten starren Systems — mit der Masse m und dem Tragheitsmoment J fiir den 
Schwerpunkt —, dessen Beschleunigungszustand durch die Beschleunigungen by eines beliebigen 
Punktes A und b, des Schwerpunktes S bestimmt wird, sind bis jetzt vier verschiedene graphische 
Verfahren vorgeschlagen worden. 


Das erste! setzt die Kenntnis des Beschleunigungspoles P* voraus, da die gesuchte Gerade t 
durch den Schnittpunkt der Verbindungslinie P*S mit der Antipolare von P* beziiglich des um S 


geschlagenen Kreises vom Halbmesser 9 = /I/m hindurchgeht. Dieser ist, nebenbei bemerkt, 
der antireziproke oder antikonjugierte Punkt von P* (von Federhofer? als Antipol von P* be- 
zeichnet). 


Das zweite Verfahren liefert ¢ als Antipolare in bezug auf den gleichen Kreis eines nur mit 
Hilfe des Beschleunigungsplanes bestimmbaren und deshalb gleichfalls mit der Stellung der Scheiben 
veranderlichen Punktes?. 


Die iibrigen beiden Verfahren beruhen entweder auf der Méglichkeit, die Scheibenbewegung 
in eine durch b, gekennzeichnete Translation und eine Drehung um A zu zerlegen® oder auf einfachen 
geometrischen Eigenschaften*. Beide verlangen die vorlaufige Bestimmung des Tragheitsmittel- 
punktes in bezug auf A, d.h. des antireziproken Punktes L, von A in der Antipolaritat beziiglich 
des wiederholten Kreises, und sind weitaus einfacher als die vorhergehenden, da die Stellung von Ly 
unveranderlich ist. 


Die Konstruktion, die wir unter Benutzung der Ahnlichkeit einiger Dreiecke finden werden, 
beruht gleichfalls auf der Eigenschaft des Tragheitsmittelpunktes (Schwingungsmittelpunktes). 
Bei der Betrachtung des in der Getriebepraxis ziemlich haufig vorkommenden allgemeinen Falles, 
in dem die Beschleunigungen zweier beliebiger mit S auf einer Geraden liegender Gelenkpunkte 
A und B gegeben sind, werden wir zeigen, daB unsere Konstruktion sich noch einfacher als die von 
O. Tolle erweist; denn mit L, und Lg vermeidet sie sowohl das fiir die Anwendung der ersten 
Konstruktion dieses Verfassers notwendige Aufzeichnen des Beschleunigungsplanes, als auch die 
Bestimmung von bg, die zum Voraus von allen tibrigen Verfahren gefordert wird. 


Aus dem Folgenden wird sich sowohl in unserer wie auch in den Tolleschen Konstruktionen der 
Vorteil zeigen, den einzelnen baw. die beiden Tragheitsmittelpunkte méglichst durch ,,direkte“ 
Versuche festzustellen, nimlich durch Versuche, welche jede Konstruktion von antireziproken 
Punkten eribrigen. 


2. Beschleunigung cines Gelenkpunkts gegeben. Wir erinnern uns zuerst daran, daB by ohne 
weiteres durch die augenblickliche Winkelgeschwindigkeit und -beschleunigung der Scheibe um 
den Drehpol bestimmt werden kann. Ist by gegeben oder bestimmt worden, so kann bg bekanntlich 
als resultierender Vektor von b, und der Relativbeschleunigung des Punktes S um A angesehen 
werden, welche sich ihrerseits (Abb.1) aus der Normalbeschleunigung CA = b§ 4 und der Tan- 


ss 
gentialbeschleunigung DC = b§ , zusammensetzt. 


Ist @ die Winkelbeschleunigung der betreffenden Bewegung und stellt also — I @ das Moment 
von P beziiglich S dar, so wird der lings der Verbindungslinie 4S gemessene Abstand s der Gera- 
den t von S zu I @/(m bg sin «) wobei a der von bg mit der Verlingerung von AS gebildete Winkel ist. 

1 H. Alt, Z. angew. Math. Mech. 6 (1926) S. 58. 

2 K, Federhofer, Z. VDI 74 (1930) S. 234. 


3 Q. Tolle, Z. VDI 74 (1930) S. 1453. 
4 0. Tolle, Ing.-Arch. 19 (1951) S. 355; M. T. Z., 15 (1954) S. 47. 
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Wird I=m- AS. SL, und w = bi 4/AS eingesetzt, so ergibt sich 
Si bs4 
i Ly by sina 9 (1) 


ee die s darstellende Strecke wird sich in Abb. 1 links von S finden, wenn ® Gegenzeigerrichtung 
at. 


P=-mb, 


Abb, 1. 


Abb, 2. 


Dies vorausgesetzt, verbinden wir mit S die Endpunkte D und A, des dem Vektor bg gleich- 


wertigen Vektors DA, dessen Verlangerung AS in E schneidet. Wir ziehen durch L, und A, 
die Parallelen zu DS, die sich mit 4,S in Fy baw. mit der Verlingerung von AS in G schneiden. 
Ist X der Schnittpunkt von AS und der Parallelen durch F'4 zu DA,, so ergibt sich aus der Ahn- 
lichkeit des Dreieckspaares (XF'4S, EA,S) und (SF4L4, SA,G) 


= = ES 
Ne () 


Wird mit S, der Endpunkt von bg bezeichnet, so ergibt sich aus der Ahnlichkeit der anderen 
Dreiecke EDS und SS,G die Beziehung ES/SG = DC/S,H = bt, ,/(bs sin 4), wodurch die rechten 
Seiten von (1) und (2) identisch werden. Hieraus folgt XS = s, und die gesuchte Gerade ¢ fallt 
mit der zu bg konstruktionsmaBig parallel gezogenzn Verbindungslinie XF, zusammen, die man 
durch die Kenntnis des Schnittpunktes F', von A,S mit der Parallelen durch L4 zu DS bzw. zu 
A,S, leicht feststellen kann. Die Aufgabe lést sich also durch Aufzeichnen der drei Geraden DS, 


A,S und L4F'4 (Abb. 2a) baw. A,S,, A,S und L4F', (Abb. 2b), je nachdem bg in der obengenannten 
Weise unter Benutzung des Beschleunigungsplanes gefunden oder als in S angreifender Vektor 
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gegeben wird. In diesem letzten Falle darf bg nicht willkiirlich, sondern mu8 bekanntlich so an- 
gegeben werden, da8 die Verbindungslinie der Endpunkte der ,,reduzierten’* Beschleunigungen 


bs/w? und b,/w? senkrecht zu AS ausfallt}. 


3. Beschleunigung zweier Gelenkpunkte gegeben. Nehmen wir nun ganz allgemein an, daB 
diese Bedingung fiir die Beschleunigungen b, und bg zweier Gelenkpunkte A und B, die mit dem 
leicht zu ermittelnden Schwerpunkt S ein Dreieck bilden (Abb. 3), sich als befriedigend erweist. 
In diesem Zustand mu man zuerst bg finden, indem man iiber A. das zu ABS ahnliche Dreieck 
A,B,S, zeichnet und dann von einem der beiden Tragheitsmittelpunkte aus, z. B. vom Ly, die 


Parallele zu A,S, zichen, welche A,S im Punkt Fy schneidet, durch den t hindurchgeht. 


Liegt jedoch S auf AB (Abb. 4), dann lohnt es sich folgendermafen zu verfahren: Man zeichne 
auf AB die Punkte L,, S und Lg und verbinde A, mit B,. Man ziehe durch L, und Lz 


zwei Parallelen zu A,B, die sich mit 4,S und B,S baw. in F i i 

awe 7 -in fy und Fz schneiden, deren Verbind - 
linie schon die gesuchte Gerade t darstellt. So eriibrigt sich jede rately Sit Baie ibotheon 
das letzte Verfahren fahig ist, die graphische Lésung aller Probleme bedeutend zu beschleunigen 
wie z. B. die Bestimmung des Tragheitspols?, bei welchen die Kenntnis der GréBe und Richtaae 


1 M. Manarini, Ricerche di Ingegneria, 1938, H. 3, S. 70. 
* 0. Tolle, Feinwerktechnik, 61 (1957) H. 9. 
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von P entbehrlich ist. Schneidet man tibrigens t mit einer weiteren, zu 4,B, durch S gezogenen 
Parallelen sowie mit A,B, selbst, so erhalt man die Endpunkte M und M, des mit bg gleichwertigen 


und daher fiir die Ermittlung der noch unbekannten Elemente geeigneten Vektors MM,. 


4. Bemerkungen iiber die Tragheitsmittelpunkts-Bestimmung. In allen vorhergehenden Kon- 
struktionen ist nun zu beachten, da8 die Strecke AL, die Lange des mit dem in A aufgehingten 
isochronen physikalischen Pendels darstellt, welches um die normal zur Bewegungsebene durch A 
gehende Achse kleine Schwingungen ausfiihrt, und daB daher 


AL, = cs Te, (3) 


ist, wobei mit T', die entsprechende Schwingungszeit bezeichnet wird. 

Liegt der betreffende Kérper in Ausfiihrung vor, und ist wegen der Unregelmafigkeit seiner 
Form die experimentelle Bestimmung von T, notwendig, so ist es vorteilhafter, statt diesen Wert 
zu benutzen, um zuerst J und dann 9g zu erhalten, den gemessenen Wert von T', gleich in (3) ein- 
zufiihren und sich ebenfalls jede Konstruktion von antireziproken Punkten zu ersparen. Nimmt 
man noch an, dafi AL, das Verhaltnis zwischen einem Tragheitsmoment und einem statischen 
Moment ist und daher unabhangig von der Dichte des schwingenden Materials, dann versteht man 
auch, daf der obengenannte T'4-Wert an einer einfachen Schablone aus Blech oder Pappe gemessen 
werden kann, wofiir Gleichung (3) auch dann anwendbar bleibt, wenn von dem Kérper nur die 
Zeichnungen zur Verfiigung stehen und man sich die bekannten zeitraubenden graphischen Methoden 
ersparen will. Handelt es sich also um ein homogene Scheibe gleicher Dicke, so hat diese Schablone 
den gleichen Umri® wie die Scheibe; im Falle irgendeines anderen gleichartigen Kérpers muB sie 
aber so profiliert werden, daB die auf Kérper und Schablone — mittels Zylinderflachen gleiches 
Halbmessers, deren Achsen die Aufhingeachse bzw. die durch A gehende Achse sind — durch- 
geschnittenen Flachen gleich sind. 

Bleibt der Beschleunigungszustand durch b, und bg bestimmt, so ware es lohnend, auch T, 


zu messén und BL» durch die auf die Schwingungen um B beziigliche Formel zu erhalten. Auf 
diese Weise kann man S bestimmen, entweder als Schnittpunkt der beiden Strecken AL, und BL 
oder, wenn diese tibereinstimmen, als Endpunkt der Senkrechten zu AB durch den Punkt N, 
in welchem sich die beiden Halbkreise, deren Durchmesser AL, und BL, sind, schneiden. 


Diese zuletzt in Abb. 5 angegebene Konstruktion 1aBt sich dadurch rechtfertigen, daB sich als 
Folge der Gleichungen 


0? = SL,- AS = SL: BS 

NS = 0 ergibt. Man kann hier auch der analytischen Ermittlung den Vorzug geben, indem man 
die Verhaltnisse €, = AL 4| AB und €, = BL,/AB herstellt und auf AB vom A ausgehend die 
Strecke AB (1 — &,)/[2—(€4 +&g)] abschneidet, die, wie man leicht feststellen kann, AS darstellt. 

Auf jeden Fall ist sowohl die Konstruktion nach Abb. 5 als auch die besagte zuphylisehe Be- 
stimmung bei weitem allen anderen Methoden, um S aufzufinden, vorzuziehen, da man die cae 
gungszeiten Ty und Tg mit geniigend groBer Genauigkeit erhalten kann, indem die Messung au 
eine Reihe zahlreicher Schwingungen ausgedehnt wird. 


(Eingegangen am 2. Mai 1960.) 
Anschrift des Verfassers: Professor Ing. Filippo Manna, Casalnuovo di Napoli (Italien), Corso Umberto, 389. 
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Berichtigung 


zu meinen Aufsatzen ,,Die Lenkstabilitat des luftbereiften Kraftwagens gegeniber 
kleinen Stérungen“ in Band 27, S. 88, und ,,Der EinfluB der Kraftwagenfederung auf 
die Lenkstabilitat“ in Band 29, S. 100 des Ingenieur-Archivs. 


Von G. Mitterlehner 


Zur erstgenannten Arbeit ist zu erginzen, daB bei Bericksichtigung der Tragheitskopplung 
zwischen Lenk- und Gierschwingung in der Gleichung (lc), in der dritten Gleichung (14) und 
in der dritten Gleichung (25) y durch (7 — ) zu ersetzen ist. Dadurch werden die Ergebnisse der 
Abschnitte 7 bis 9 nur bei hohen Frequenzen der Lenkschwingung brauchbar. In Abb. 7 fallt der 
ganze erste Quadrant der (n, v)-Ebene in den stabilen Bereich, und in Abb. 8 wird die Stabilitats- 
grenze eine Gerade v = konst. 


In der zweiten Arbeit verandert der durch die dritte Gleichung (4) eingefiihrte Fehler je nach 
Fahrzustand und Konstruktion die Koeffizienten der charakteristischen Gleichung nur geringfigig, 
da in diesen Beziehungen die von den Kreiselmomenten der Rader stammenden Glieder tiberwiegen. 


(Eingegangen am 21. Dezember 1960.) 


Anschrift des Verfassers: Dr. G. Mitterlehner, Wien II, Halmgasse 1/20. 
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ZUR INFORMATION 
Moderne Turbomaschinen-Anlagen stellen oft schwierige regeltechnische Probleme, deren Lésung eine 
sehr genaue Kenntnis des Regelverhaltens voraussetzt. Dieses wird daher im vorliegenden Buch ein- 
gehend behandelt. Weiter sind die Festigkeitsprobleme des Turbomaschinenbaues vor allem im Zusam- 
menhang mit den hohen Temperaturen und Driicken sehr komplex geworden und erfordern eine ent- 
sprechend genaue Erfassung. Endlich zwingen die in der Praxis immer noch gelegentlich auftretenden Lauf- 
stdrungen und Schaufelbriiche dazu, die dynamischen Berechnungen viel weiter zu treiben als dies frither 


iiblich war. Alle diese Fragenkomplexe erfahren daher in diesem zweiten Band eine griindliche Darstellung. 
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